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ABSTRACT 

Mass transfer studies with a single sphere absorber were used to 
determine the diffusivity of hydrogen in liquid aminomethane over a temp- 
erature range of -20°C to -30°C. The diffusion coefficients were found 
to lie in the range of 90 x 107° to, 110 x 10° chic /sec. A modified 
Ferrell and Himmelblau correlation has been derived from the diffusion 
data for hydrogen in ammonia and in aminomethane. It is recommended 
for use in estimating the diffusion coefficients for hydrogen, hydrogen 
deuteride and deuterium in aliphatic amines where experimental data is 
not available. The accuracy of the diffusion data obtained is estimated 
to be within £5%. 

Pseudo-first order reaction rate constants for the liquid phase 
hydrogen deuteride aminomethane exchange were determined over the tempera- 
ture range of -10°C to -30°C. This was achieved by measuring the rates 
of absorption of hydrogen deuteride into a thin film of liquid amino= 
methane flowing over a sphere absorber operating in the fast reaction 
regime. Potassium methyl amide in solution with aminomethane was used to 
catalyze the exchange process. Rate constants were found to be in the 
range of 230 to 5/70 Bet at catalyst concentrations ranging up to 15.8 
(gms K/kg Amine). The rate constant appears to vary linearly with 
catalyst concentration up to a value of 5.0 (gms K/kg Amine). The activa- 
tion energy of the reaction was estimated to be 3.1 (kcal/gm mole). An 
expression was derived which represents the data in the range of investi- 
gation to £4% and is recommended for predicting the rate constants at 
conditions of temperature and catalyst concentration where data are not 
available. 


It is estimated that the accuracy of the rate constants obtained is 


within 410%. fg 
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CHAPTER | 


INTRODUCTION 

The isotopic exchange of deuterium between gaseous hydrogen and 
liquid aliphatic amines has considerable potential as a basis for a 
commercial process to produce heavy water. The economics of this 
route indicate that there is promise to produce heavy water more 
cheaply than processes based on the exchange between hydrogen-ammonia 
or hydrogen sulphide-water. Amongst the aliphatic amines, amino- 
methane has been reported to be the most attractive [1]. 

In order to design industrial contactors for systems such as 
hydrogen-aminomethane which undergo a simultaneous absorption and 
chemical reaction, it is desirable that sufficient information for 
estimating mass transfer coefficients be available. The diffusion 
coefficient for physical absorption and the kinetic rate constant are 
two parameters which are useful for estimating transfer coefficients 
when chemical reaction is involved. 

Physical gas absorption into a liquid flowing over a sphere has 
been successfully used by a number of workers [2,3,4,5,6] to obtain 
estimates of the diffusion coefficients in gas-liquid systems. The 
single sphere contactor minimizes the end effects common to other gas 
liquid contactors. It is also reported [3,7] that relatively high 
liquid flow rates, with absence of surface rippling, may be employed 
with this absorber. Compared to the laminar jet, the single sphere 
absorber provides longer contact times (0.1 sec to 1.0 sec) and a 


higher surface area per unit volume of flowing liquid which makes 
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the sphere absorber more suitable for studies involving sparingly soluble 
gases. 

Preliminary measurements of the diffusion coefficient of hydrogen 
in aminomethane were made by Moore [8] but the results obtained indi- 
cated inconsistencies which were attributed to experimental difficulties. 
Subsequent refinements of the experimental apparatus and additional 
operating experience have allowed reliable results to be obtained. 


These results are reported as part of this investigation. 
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CHAPTER 2 


LITERATURE SURVEY 
A fairly comprehensive literature review of the various experi-~ 
mental techniques employed in the measurement of diffusion coefficients 
has been undertaken by Moore [8]. A repetition of that material is 
deemed unnecessary for this investigation. A review of the operating 
experience pertinent to the sphere absorber and details of the math- 


ematical techniques used to analyze the data obtained will be presented. 


2.1 Past Experience with the Single Sphere Absorber 


A single sphere absorber consists of a sphere over which a liquid 
flows in a laminar film. The liquid is introduced through an orifice 
onto the top pole of the sphere and is collected in a constant liquid 
level take-off tube situated directly below the sphere. The sphere 
is enclosed in a gas tight chamber which is in turn located in a constant 
temperature environment. The molecular diffusivity can be calculated 
from the measured rate of absorption of a gas into the liquid film on 
the sphere. 

Past experience with the design and operation of an absorption 
system using the single sphere absorber has indicated the following 
important observations: 

(i) A large volume of the gas absorption chamber results in 
excessive errors in the measured gas absorption rate caused by minor 
temperature fluctuations of the gas phase [5]. 


(ii) A more stable surface film of liquid over the sphere may 
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be obtained by using a laminar jet as the liquid feeding device [7,9] 


in preference to having the support rod extend through the liquid feed 


nozzle [3]. 


(iii) The surface film is very sensitive to vibrations of the 


apparatus. Use of vibration arrestors is recommended. 


(iv) The sphere size should be chosen in a manner to eliminate 


wall effects of the absorption chamber and at the same time prevent 


the effects indicated in (i). 


(v) Depending on the system being studied, there is: 


(a) 


(b) 


(vi) The use of the single sphere absorber has associated with it 


an upper and a lower limit to the jet length that 
should be used; and 

a range of permissible distances between the sphere 
and the take-off level for which the stagnant layer 
reported by previous workers [2,3,4,7], may be main- 
tained on the support rod. Shorter distances cause 
the stagnant layer to move on to the sphere absorber 
thus reducing the effective interfacial area and 
resulting in a marked decrease in the measured 


absorption rates. 


some end effects which may be subdivided as follows: 


(a) 


the fresh liquid issuing from it will absorb some solute gas. 


Entrance Effects 


When a laminar jet is used as the liquid feeding device, 


of gas absorbed on the jet is reported to be significant [9]. This 


amount may, however, be calculated using the equation for physical 


absorption by an ideal jet [9,10,11]: 


The amount 


ee _ _ > = 
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5. 
Gracey vh LD (Cc. ~ Cc’) (1-1) 


Assuming complete mixing of the solute gas and the liquid at the 
sphere-jet junction [9], the bulk concentration of the liquid flowing 


on to the sphere may be evaluated as 
G = we (1-2) 


Goettler [9] investigated the cases of complete mixing and no mixing 
at the sphere-jet junction while studying carbon dioxide absorption in 
water. His results indicated that complete mixing at the sphere-jet 
junction was a more likely representation of the fluid conditions at 
the top pole of the sphere. 

(b) Exit Effects 

If the distance between the bottom of the sphere and 
the take-off level is properly adjusted the absorption of the solute 
gas on the exposed support rod is reported to be negligible due to the 
presence of a stagnant layer [2,3,4,9]. This stagnant layer is attributed 
to the presence of trace impurities of surface active agents in the out- 
going liquid [4,9]. As a consequence of this layer, it is suggested 
that the velocity profile of the liquid over the support rod might be 
fully parabolic. This would reduce the gas absorption on the rod 
sufficiently to be neglected [4,9]. 

Experience with the sphere absorber as reported above is limited 
to studies involving water as the absorbing liquid wherein the experi- 
ments were conducted at temperatures generally within 15 to 30°C. Moore's 
investigation [8] has indicated that use of the single sphere absorber 
for low temperature studies involving the hydrogen-aminomethane system 


requires that emphasis also be placed on the following aspects of its 


operation: 


2 
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6. 

(i) Minor temperature fluctuations of the absorption system or 
the room where the experiment is conducted can cause large errors in 
the measurement of the gas absorption rate. This requires that: 

(a) The experiments be conducted in a controlled 
temperature laboratory. 

(b) The incoming liquid feed be pre-cooled to a temperature 
which is close to the constant temperature environment 
surrounding the absorption system. 

(c) The liquid leaving the absorber must be maintained at 
a temperature equal to, or lower than, the temperature 
of the absorption system to prevent inverse pressure 
gradients from impeding flow. 

(ii) The take-off level control device (lute vessel) used by 
Davidson et al. [3] may not be used for this system. A fine metering 
valve on the outgoing liquid line is a suitable substitution. 

(iii) The choice of the size of the sphere absorber should take 
into consideration the saturation levels of the outgoing liquid. Very 
high saturations can cause the errors due to temperature fluctuations 


and measured gas rates to be magnified considerably. 


2.2 Theoretical Analysis 
2.2.1 Hydrodynamics of the Spherical Film. The hydrodynamic 


behaviour of the spherical film of liquid flowing over the sphere 
absorber may be described by using the solution of Nusselt [12] for 
laminar flow down an inclined plane. This approach was first attempted 


by Kramers et al. [2] and was subsequently verified and used by other 
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investigators also [3,9,13,14]. 

It may be assumed that the film thickness of the liquid f lowing 
at any latitude on the sphere is the same as it would be for the same 
flow rate per unit length on an inclined plane making the same angle 
with the vertical. Thus, for a liquid flowing over the surface of 


a sphere (Figure 1) the liquid film thickness is given by 





1/3 =2/3 
a Sie : = 
A | aaRa | sin ) (133) 
and the velocity distribution is represented as 
- 2 
Vy = e Gite Ore ) (1-4) 


where te is the velocity at the surface of the film (xe="0) and is 


given as 
ZA dat 
Pree. ee Sil Td < 
us 2v (1-5) 


Combining equations (1-3) and (1-5) one gets 


9 1/352/3 5-2/3 -1/3 
V o«= L R sin ) (1-6) 
re) 2 

32vT 


For an element of liquid flowing over the sphere ® varies from 0 to 7; 
the diffusion or contact time is therefore given as 


T 
7 Rd@ _ 9g 
t = — = 
V 2 
0 re) 32a 
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FIGURE 1 


Cross Section of Sphere and Liquid Film 
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sin 6 dé 


1¢273)me (1/298, 
mah = 2558 (1-8) 
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2.2.2 Physical Absorption of a Single Gas ina Spherical 


Liquid Film. The equation representing the physical 
absorption process of a gas in a spherical film may be derived by a 
simple steady state mass balance around an elemental section ABCD of 
the liquid film (Figure 1). 

In making this balance it is assumed that, (i) transport by 
molecular diffusion in the direction of bulk flow, i.e., across the 
radial lines AB and CD is negligible, and (ii) transport of mass by 
convection, i.e., bulk flow across the streamlines AC and BD is 
negligible. 

The net transport of material across the radial lines AB and CD 
by convection is: 

oC 


27R sin 6 Vy a] dx dé (1-9) 


The net material transferred by diffusion across the streamlines 


AC and BD is: 


; . 
ray tin eel ms - 2R a dx dé (1-10) 
aX 9 9 


Under steady state conditions equations (1-9) and (1-10) may be equated 


to give: 
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Since the liquid film around the sphere is usually very thin, it may be 


assumed that: 


2 

9“C 2 | 9c 

s ea aca (1-12) 
C6 Q 


2 V 
eC oC 
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The above equation has been solved by a number of authors [2.351547 

The analysis of Davidson and Cullen [3] has been used for this investi- 
gation. The above authors considered the effect of the stretching of 
the liquid film and introduced a dimensionless length y = ~ into 


equation (1-13) to obtain: 
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2 Vee 8 
ceo aoe oc (1-14) 
wa RD | 36 


Substituting for y from equation (1-4) one obtains: 
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Define a stream function w(y) as 
x 


v(y) = | V. 2n(R+A-x') sin 6 oe (1-16) 


o o.! 


By substituting for Yi from equation (1-4) and neglecting A and 


x' in comparison to R equation (1-16) may be represented as 
v= SL (y - y3/3) (1-17) 


Since ~ is a function of y only one may say that 


(3), -(38, 0 


and defining 





R 
do 2 B, 
Ea Ven (1-19) 


One obtains 


2 
aC 22 dC 
$ 


which is the final form of the differential equation derived by 


Davidson and Cullen [3]. 
The above authors considered two sets of boundary conditions to 


cover cases where the depth of penetration of the solute in the liquid 


film may be low or high. 


(a) Low Depth of Penetration Model 


For this case y << 1 and equation (1-20) becomes 
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The solution of equation (1-21) with the B.C. set (1-22) is 





y 
(c - C) = (C. - c) erfc a (1-23) 


The total rate of gas absorption may thus be obtained from 


9 
 Neeaa te yi) le a 
c= Zu | [rlcet (1-24) 
0 Fe . 
1/3 7/3 
For a sphere $, = 5 x 1.687 E 73 (1-25) 
u 
Thus 
172 a 1/ bye 7/6 1/372 
Gate (12 x 1.68) i R L D (C.-C) (1-26) 


(b) High Depth of Penetration Model 
In this case a may not be neglected in comparison 


to 1 and the B.C. are: 


i} 
QO 


(1) Jo. = 0; y 2 0; C 
(2) > 0; y = 0; Cece Cie 7) 


(3) $20; yal; Se =o 


B.C. (3) merely states that no solute crosses the solid sphere boundary. 
Davidson and Cullen [3] adapted the solution of Pigford [16] and 


Vyazovov [17] for gas absorption in a liquid film flowing down a 
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wetted wall column. The solutions for the single sphere and the 
wetted wall column are equivalent so long as the value of bo is 
the same for the two equations. 


The total gas absorption rate may be expressed as a series 


given by: 
n 
BaeeeiCeee ) al lhe & B; exp(-y,o) ] (1-28) 
or |G. L(C.-C.) F (a) (1-28a) 
where a = 3Do, (1-29) 


The coefficients B. and y,; are reported in Table 1. 

Davidson and Cullen [3] recommend use of equation (1-28) with 
coefficients as in Table 1 for cases where the outlet liquid satura- 
tion exceeds 40%. 

Olbrich and Wild [15] extended the solution of Davidson and 
Cullen so that it may be used for both low and high depth of 
penetration cases. Their solution corresponds to equation (1-28) 
with ten terms being used in the series whereas Davidson and Cullen 
used only four. For high depths of penetration the calculated gas 
absorption rate using the coefficients of Davidson and Cullen or 
Olbrich and Wild is almost identical. 

Equations (1-26) and (1-28) consider the surface area of the 
liquid in contact with the gas to be equivalent to the surface area 
of a dry sphere. Goettler [9] and Ratcliff and Reid [4] pointed 
out that the actual surface area of the liquid is larger than the 


area of a dry sphere owing to the thickness of the liquid film 
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Coefficients for Series Solution for G 
(Davidson and Cullen [3]) 
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TS. 
flowing over it. Goettler used an average correction which was 
based on the average liquid depth over the complete sphere. The 


correction is of the form 





+ =i) Eo ess = (1-30) 
where 

An = corrected interfacial area 

A. = surface area of the dry sphere 

An = depth of the liquid film at the sphere equator given as 

pee i re cg 


Equation (1-26) was used to estimate a starting value of 
diffusivity and this first estimate was used in equation (1-28) to 
generate an iterative scheme which was considered to have converged 
to a final value when two successive estimates of diffusivity did 
not vary by more than 0.5%. 

Calculated diffusivities are very dependent on solubility data, 
and accurate values can only be obtained if very accurate solubility 
data are available. This is evident from equations (1-26) and (1-28) 
where C. may be taken to be a saturation concentration of the gas in 
the liquid at the gas-liquid interface. The solubility of hydrogen 
in aminomethane at pressures up to 20 atm. and over the temperature 
range -62°C to 23°C has been determined by Moore [8]. The Henry's 


Law constant is given by 


In H, = 4.1403 + 1.66128 - 0.11608" (1-32) 
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B = TUK) (1-32a) 


and “P. = x.H (1-32b) 


The viscosity of aminomethane has been determined as a function 
of temperature by several authors [18,19,20]. This data has been 
assimilated and a mean curve of viscosity vs. temperature has been 


reported by AECL [21]. 


fies Empirical Methods for Prediction of Diffusion Coefficients 


Moore [8] has reviewed a number of correlations employed for 
the purpose of predicting diffusion coefficients of sparingly 
soluble gases in non-aqueous solvents. Only those correlations 
that have been used for comparison with data obtained during this 
investigation will be given here. 


(a) Wilke-Chang Correlation [22] 


The empirical form of this correlation is 


1/2 
Du aM 
1“ 20 Me ( ee li) 
Tami rhea ae) yore (1-33) 
2. 


where ‘a' is an "association parameter'' which defines the effective 
molecular weight of the solvent. The "'association parameter'' must 
be determined from experimental diffusivity data for the solvent 
in question. Wilke and Chang [22] report a value of 2.6 for water 


and 1.0 for non-associated solvents. 
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(b) Lusis and Ratcliff Correlation [23] 


The Lusis and Ratcliff correlation is of the form: 
ae Ps Walia a4 
Tae 8.52 x 110 (V,) 1.40 = + |j— (1-34) 


This equation is reported to be applicable for the entire range of 
solvent and solute molar volumes but its applicability at tempera- 
tures substantially below 0°C is questionable owing to the lack 
of available data at low temperatures. 

(c) Ferrell and Himmelblau Correlation [24] 

This correlation was developed for the solution of 

sparingly soluble gases in water and as such one cannot expect it 
to predict accurate data for a non-aqueous solvent. In spite of 
this, it may be used to provide reasonable estimates of the 


diffusivity of hydrogen in ammonia. The correlation is of the form: 








a 
Duy 5} Ee oe 
= 4.8 x 10 (1-35) 
it V 
2 
=D 
where a, was reported by the authors to have an equality with EE 
U 
if cr is computed from the reJation 
es eal 
a iz - 
Vv \" (1-36) 
we 
N 


o in equation (1-36) is the collision diameter of the solute molecules. 
A’ is the quantum parameter of Hildebrand et al. [25] which 
accounts for the presence of quantum effects when gases such as 


helium or hydrogen diffuse through liquids with high internal pressures. 





[ES] noiseistye9 Vitl>seA bas 2taud a) 


“mor sl) to 2] poltsisrias THIER bas 2iavl SAT 7 7) | 






i vy) EM ey) ] EN us 
iets Pas | 4 R-~ i ' 
(dg-1) tc a eitg CW) | Of *« 82.6 a 


\ { 
S Spne1 syiine aft 10% sideailqgs. sd of badvoge? 27 noijeipe elaT 
“syeqms) Js yiitidesifaqe 211} Jud eamilov islom stufoe baa Jnevios 


A265! of12 of pniwo asidenoidzsup 21 3°O woled yilsisnesedie gov 


[AS] noiz6fsyx09 veldismmit} bie Iistzey (9) 


St ae 


tO neiluioe Shs 70% bsgolevsb zaw noltslaiias 2iAT 


; 
.2aiuisisqmes wol gs 6ieb sidetisvae to 


Ji j9eqgxe tonnes sno dope 26 bas isd6w Al 2526p alduloe yipntisqe 


tO siige ni .jnevioe evosups-non 6 Toh. eteb siswiaos Ioibeiq oF 


aij to eetemites sldenotss1 sbhivorg of baeu sd yam 3] «zing 





2mI07? Si Jo zi noijs!s1709 4dT .sinomme ni asporbyd Fo yJivizutiib 
D 
(2e-1) ie ee i So f. 
" fi x 5.8 = 
ee | wer 
a - ; - 
| ¢ 
5, Atiw ydileups ne ove oF etoliue sila yd tetiogey e6w ;? 9 1Sriw 
g: 
nois6is7 sit moit bszvemes 2l he ii 
ar Nae: 7 





. Fi > S 7 

OF ton im _ 
_ 7 : : ia : 
iresvany 23. da 


Wes 
Ferrell and Himmelblau [24] suggested the use of the theory 
of absolute reaction rates to provide a molecular model which 
could serve as a basis for correlation of diffusivities. 
The expression of Eyring [26] for the molecular diffusion 


coefficient is written as: 


E 





D = A,T exp (- ar ) (1-37) 
E 
D 
d Deals ‘ 
we [Ct ja ies 


A similar treatment of the corresponding expression for viscosity 


yields 
d ale 
Sn Glan a) A aoe (39) 
RET? 
d D 
, a8 [i a) ip eB 
OS tn W)] Fy 
= 
o fine 2 | - 595 = [ln (u)] (1-40) 


Integrating equation (1-40) with respect to T one obtains: 


in (2) = —" In(u) + In (A,) (1-41) 


where In (A,) is a constant of integration. A, is reported to be 
independent of temperature but it is a function of the properties of 


the solute and the solvent. 
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Accepting the equality of a with a equation (1-41) may be 
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D= ee (1-42) 


Once again AS is a function of the properties of the solute-solvent 
pair and u is the viscosity of the solution. For sparingly soluble 
gases wu may be taken to be the viscosity of the solvent also. 


The quantum parameter A* of Hildebrand et al. [25] may be 


defined as: 


Ax = a (1-43) 


o(me) 


Ferrell and Himmelblau postulated that for sparingly soluble 


gases in water A, would be represented by 


Anes a aay goer (1-44) 


Values of 'a' and 'b' were determined by linear regression and the 
resulting expression for Ay when substituted in equation (1-42) 
yielded the final form of the correlation represented as equation 
(1-35). The authors have reported a much better agreement of 


equation (1-35) to their data than the correlation of Wilke-Chang 
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CHAPTER 3 


EXPERIMENTAL 


3.1 Apparatus 
Moore [8] designed and built the basic equipment for the measurement 
of molecular diffusion coefficients of hydrogen in aminomethane. 
As mentioned earlier, preliminary measurements made by him 
indicated experimental difficulties. The apparatus was thus 
modified to improve its performance and more reliable data was 
obtained with it for this investigation. Detailed description 
of the equipment used may be obtained from Moore's thesis [8] 
and only a short review of the set-up is included here. Details 
of supporting equipment used such as: 
(a) The refrigeration system 
(b) The hydrogen purification system 
(c) The aminomethane purification unit 
(d) The vacuum system 

are given in Moore's thesis. 

A schematic diagram of the single sphere absorber is given in 
Figure 2. 

The design of the single sphere absorption chamber was similar 
in concept to the absorber of Goettler [7,9] in that a laminar jet 
was used for introducing aminomethane onto the sphere. A number 
of modifications on the equipment used by the above authors were 
necessary due to the properties of aminomethane. A diagram of the 


absorption chamber is given in Figure 3. 


— RATIAKD 


JAVMIMIASIXS 


| ZOeeISaED) ‘*at 


Shamstvesem oft vO? Josmyiups >)zed ens jf ivd bd banoieeb [8] s100M 
shedismeains ni nspotbyA to 2insisitisos soleutilh yefgssian Fo 
win ea abam z2zinsmsyuebem yiarimbisig vit inés benolinam 24 
euns 264 2udsieqqe sy -@84a ipoittib. lesnem)vadxs bsgeot bat 
tev Steb sidsifsy siom bos sonamottsq ef) SV¥OTGM) oF BeTTigGm 
moliqitsesb bslistsd .niclispiszaval gids vot 2) ddiw Genieids 
18] elezsni 2's xe0% mor} bsnisido sd yen bsev Idamgiwps efds To 
élisse@. .sisr bsbulsai <i du-tsa sit to weivsx 216ne & vind One 
:26 lave bseu Indmalups eplsveqque to 
majzye noljsispliiat sat (es) 
msjJay2 aolisoitiawe meporbyd ont. ta) 
jiu noitsciti wd snsdeenonine ohh. “fag 
natzve muuosvesdT tb) ty 
-2izud), 2! StooM ni nevig a6 
: pp fevig zi iedtoeds evarnge slont2 dy To merpeib Sigemedse Ain 





~ 4S s1ught 
| 
=i 
7 . 
| | 
) 












VACUUM 
SYSTEM 


= 
S) 
G 
= 





PURIFICATION 
SYSTEM 


Ze 
uw 
O 
O 
a 
O 
> 
als 


P=A0psig 


oO 
uw 
re 
a 
=) 
ae 





E. 
DG x) 4 


Spee 





$=) (8) 


Ne Tho eee UA 
SENG SSS NSO OSE 


Q) 
J 
dX 


FIGURE 2 


Schematic Diagram of Single Sphere Absorber 
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Symbols for Figure 2 
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Liquid Amine Feed Tank 

Relief Valve 

Glass Wool Filter 

Liquid Feed Chiller 

Liquid Feéd Rate Control Valve 
Liquid Feed Chiller 

Liquid Rotameter 

Liquid Distributor 


Single Sphere Absorption Chamber 


. Take-Off Level Control Valve 

. Liquid Outlet Chiller 

. Spent Liquid Amine Dump Tank 

. Hydrogen Low Pressure Regulator 
. Soap Film Meter 

. Hydrogen Vent Valve 

SO dp | fap 

. Hydrogen Circulator 

. Hydrogen Saturator 

. Solenoid Valve Sy 

. Solenoid Valve So 

. Solenoid Valve (Hydrogen Saturator) 
. Thermocouples 


. Mercury Manometer 
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FIGURE 3 


Schematic Diagram of Absorption Chamber 
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The laminar jet was designed in a similar manner to those of 
Goettler [9], Scriven [27,28] and Hatch [29]. The nozzle design 
was modified to incorporate the use of ruby watch jewels for the 
cylindrical throat sections. Their use was considered desirable 
because they could be accurately bored and press fitted into the 
nozzle caps. The nozzle diameter used for this work was 0.046 cm. 

The spheres used in this study were 1] inch and 3/4 inch 
diameter ball bearings which had been softened by heat treating. 
They were suitably drilled so that they may be mounted on a 1/8 
inch diameter shaft or support rod and had a 1/16 inch diameter 
mark machined on the top pole to facilitate aligning the laminar 
jet. Precise alignment of the laminar jet over the top pole of the 
sphere was achieved by an eccentric mechanism mounted in the top 
flange of the absorption chamber. Acetone was passed through the 
jet during the alignment procedure. 

Liquid flowing over the sphere and the support rod passed 
through a glass collector cap mounted on a take off tube. This assembly 
passed through the absorber base flange before being connected to the 
liquid return line. The distance between the bottom of the sphere 
and the level of the outgoing liquid in the glass collector cap 
(take-off height) could be varied by moving the steel support rod 
up or down. 

Details of the top flange housing, the eccentric mechanism and 
the bottom base flange may be obtained from Moore's thesis [8]. 

The support rod was roughened to enhance wetting and was 


centered in the liquid take-off tube by two teflon guides. The rod 
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24. 
passed through the base of the take-off tube and was attached to the 
frame supporting the entire absorption apparatus by a threaded 
adjustment nut which was used to adjust the sphere height relative 
to the floor of the absorber. 

A 1/2 inch Swagelok connector was used to seal the glass 
connector cap to the top of the 1/2 inch stainless steel take-off 
tube. The take-off tube in turn was sealed to the base flange 
by an O-ring seal. 

All the equipment, as shown in Figure 2, was fabricated out 
of 316 stainless steel except for those parts where visual 
observation was necessary. The latter parts were made of QVF 
glass sections and Pyrex glass. Neoprene was used exclusively 
for all gasketing material and O-rings. Zytel (Nylon) ferrules 
were used for glass to stainless steel Swagelok fittings. 

The absorption chamber and the hydrogen saturator were enclosed 
in a refrigerated air bath to enable operation in the temperature 


range of -30°C to -20°C. 


3.2 Operation and Procedure 

Fresh aminomethane was purified by contacting with lithium 
metal to remove water followed by a simple flash vacuum distillation. 
Aminomethane has a normal boiling point of -6.33°C and thus it could 
be fed to the absorber under its own vapour pressure by maintaining 
the feed tanks close to room temperature. The liquid feed was 
precooled before entering the refrigerated air bath to minimize 
thermal gradients, and after leaving the absorber it was collected 


in a dump tank immersed in liquid nitrogen. The amine was subse- 
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25% 
quently degassed and reused. In order to obtain an unimpeded flow 
through the exit liquid line to the dump tank, it was necessary 
to chill the liquid line through the use of a concentric tube cooler. 
When this was not done, a back pressure was set up in 
the liquid line and the flow of aminomethane was impeded. The 
liquid feed and liquid outlet control valves were both ‘Whitey 
21RS4-316'' micrometering valves and their use allowed excellent 
flow control. 

Hydrogen was introduced through a ''Matheson Model 70A'! low 
pressure regulator. The absorption chamber was pressurized to a 
preselected system pressure which was recorded with a mercury 
manometer and the gas phase was saturated with aminomethane prior 
to a run. The hydrogen absorption rate was measured with a soap 
film meter constructed from a 50 c.c. burette with 0.10 c.c. 
graduations. 

The operating variables that required consideration were 
liquid flow rate, absorber temperature, absorber pressure, sphere 
size, distance between the sphere and the feed jet (i.e. jet length), 
and the distance between the base of the sphere and the liquid level 
in the take-off tube (i.e. take-off height). 

The liquid flow rate should be such that laminar flow is 
maintained on the sphere surface. Although the flow of liquids 
in a thin film over a solid surface has facets a large amount 
of attention [30,31], the criterion for conditions that will eliminate 
convective mixing is not well defined. Previous work with gas- 


water systems [3,4,9] based on a Reynolds number criterion evaluated 
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268 
at the sphere equator suggests that a film of aminomethane should 
be laminar at flow rates upto 1.0 cc/sec at temperatures from 
-15°C to -30°C. However, experimental data which agreed with the 
absorption model could only be obtained at liquid flow rates upto 
approximately 0.45 cc/sec or at Reynolds numbers less than about 
50. The occurrence of surface activity was much more apparent on 
the support rod than on the sphere and it was generally observed 
that when ripples were present on the rod, they were also present 
on the sphere. 

The Reynolds number criterion, however, does not incorporate 
the effects of surface tension of the liquid which appear to affect 
the rippling phenomena. The surface tension of aminomethane is 
about one third that of water and this significant difference 
could well be responsible for reducing the upper limit of the 
Reynolds number for stable flow over the sphere absorber. The 
effect of surface tension needs to be studied in more detail than 
has been done at present before conclusive evidence in this direction 
is available. 

The laminar jet was operated within its upper and lower limits 
of stability. These limits have been defined and observed by 
Goettler [9] and reconfirmed by Moore [8]. A jet length of 2.0 mm 
and a take-off height of 1.25 cm was used throughout this investigation. 

The startup and operating procedure for the single sphere absorber 
is essentially the same as that reported by Moore [8] with some 
modifications. Details of this procedure are included in Appendix Bl. 

The manner in which the data obtained from a run was used to 


determine the diffusivity of hydrogen in aminomethane using equations 
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(1-28) and (1-29) has been described in detail by Moore [8]. From 
the measured values of the different variables indicated earlier 
like liquid flow rate, temperature and hydrogen absorption rate, 
the molecular diffusion coefficient D was taken to be the value 


which minimized the function 


E= ) | [06.), - ((c,-c) F (a)),] | (1-45) 


where N = number of determinations. 


The experimental data is reported in Appendix Al. 
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CHAPTER 4 


RESULTS 

Experimental measurements of the rate of physical absorption 
of hydrogen in aminomethane were taken on the 3/4 inch and the 
1 inch spheres. These are reported in Tables 2 and 3 along with 
the calculated values of the diffusion coefficient. 

Figure 4 is a plot of the hydrogen absorption rate corrected 
to 1 atmosphere hydrogen partial pressure versus the liquid flow 
rate over the sphere at -29°C. Average values of the gas flow 
meter temperature and system temperature were used to generate 
the plots shown. The gas flow meter temperature is required to 
estimate the appropriate gas density for use in expressing mass 
flow rates. The maximum deviation of both these temperatures 
about the mean over a period of approximately four hours, did not 
exceed +0.2°C. The solid curves shown represent the simulated 
variation of the hydrogen absorption rate versus the liquid flow 
rate, using the high depth of penetration model, for two assumed 
values of the diffusion coefficient. This was done using a computer 
program written by Moore [8]. The plot depicts the degree of agree- 
ment of the absorption rates with those predicted by the model. 

Figure 5 shows the variation of the absorption rate with the 
length of the exposed support rod. This plot indicates that a 
length of 1.25 cms. corresponds to the case where the stagnant 
layer lies at the junction Brethe sphere and the support rod and 
shorter lengths would cause it to move on to the sphere. 


The logarithm of the diffusion coefficient as determined with 
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TABLE 2 


Experimental Results for the Diffusivity of Hydrogen 
in Liquid Aminomethane Using the 3/4 inch Sphere 


6 
L G(1 atm) x10? 2 i 1p 

T(E) (cc/sec) (gms/sec) % Saturation (cm”/sec) 
= 229 0.349 Q0...15,7 75.6 106.3 
= 25 42 0.350 Olio] 76.8 hiO.2 
-25.8 02303 0.146 80.7 105.4 
-26.6 0.351 0.156 74.9 105.0 
-27.8 0.336 0.149 75.6 OZ ea 

TABLE 3 


Experimental Results for the Diffusivity of Hydrogen 
in Liquid Aminomethane Using the | inch Sphere 


L G(1 atm) x10° eas 10° 
TCG) (cc/sec) (gms/sec) % Saturation (cm’/sec) 
=e) 0.400 OmrgZ 83.9 91.4 
S20 ee 0.378 0.188 86.5 94.1 
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-28.9 0.3/2 0.188 O63 94.8 
-26.7 0.402 0.207 86.9 LOZ .2 
-26.7 0.405 0205 84.9 S eee 
-2h .0 0.418 0.218 85.8 isles 
Pace i) 0.418 OVZ19 86.1 OZ ah 
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=205.0 0.426 e235 87.0 LOW 


220.0 0.421 O5235 87.9 11050 


S: Sd8AT 


nspeibyH to yilviabttid srg aot peas i rprend 
ejodge dani A\E add phied spstsmonine | don 





) 


of : Soix (mis Rye y _ 
(se2\" nia.) anitewise ¥ (oae Nem) (ge2\a9) | RIE a 7. 
—— — tne een ne ee : y 


£.301 d.2f Feb. 0 eve .0 2% : 





L) 
Av 
e.or! 5. at fai<0 02g, S.2S- 
A, 201 r08 att ..0 e0f.0 ~ Bae : 
a 
0.201 2, AY 320.0 f2e.0 2. a8- 
0 5.2 CF T.0 d££.0 8.7S- 
E 3eAt 
y neporbyH Fo yiivievt710, sia 107 27ldesh lésqemlyeees: 
eisrigé toni | sig eniel snsnismontimA bitipis at » 
ae “Oix (mas 10 4 
(9s2\ ins) wisJeiuiEe F (sa2\eng) (pa2\o5) (3°9T . 





Pete e088 set.0 004.0 0. 0€- 7 
a 2,98 sot.0 8te.0 8.85 aS 

Pee | e, a8 se 1.0 POe.0 0.88- 
8 Xe £.48 881.0 SVE.0 2,.8g- 
$.S01 ee: GOs.0 £0F.0 1.85, 
1,ee e886 Ce 10H > §.a8- _ 


e 
a a 





a 

eee ae 
“ro 

ot 


“he 


CORRECTED HYDROGEN ABSORPTION RATE (gm/sec)x10°> 









D =100«107° cm2/sec 
D=94x10°° 






cm2/sec 











JL | cole be slate aa 
Cc Geile ele 
Bed SI el 1 


0.2 0. 
LIQUID ey RATE a 


FIGURE 4 





Rate of Absorption of Hydrogen in Aminomethane at -29°C on a 1 inch Sphere 








r ' -a ' : - - 7 
ie = = 7 A ' 7 > 7 
u 1 33 ; ‘ ‘ 


- *Ofx(se2\me) 3TAI UOITIROZBA UIOORGYH G3TDIRROD 





— 








; 
I 
| 
+ cua 





' 
a 





| 
——s 
' 


ae @ 


ra 
} 
a CAGE Seeeneent meee comer tae 
A: 


} 
t 
‘ 

i 





2°08 5 
™ 
7 


nl ee 
ene 
> 7 
ai 
ef 


~  k0 


SAUTARISM3T ADARAVA 


. 
3 





of 
a 


CORRECTED HYDROGEN ABSORPTION 


RATE 


(gms /sec) x10 























0.230 
0.220 
0.210 
0.200 
0.190 
L 
( cc/sec ) 

0.180 4 0.418 

: @ 0.392 
0.170 ; 


0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 
TAKE - OFF HEIGHT (cms) 


FIGURE => 
Effect of Take-off Length on 


Hydrogen Absorption Rate 
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the I inch sphere is plotted versus the reciprocal absolute tempera- 
ture in Figure 6. Experimental data for hydrogen ammonia [32] is 


also shown for comparison. 
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FIGURE 6 


Comparison of Predicted and Experimental Diffusivities for 


Hydrogen in Ammonia and Aminomethane 
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CHAPTER 5 


DISCUSSION 
5.1 Experimental Data 

It is apparent from Figure 4 that for liquid flow rates where 
the absorption data agree with the theoretical model (<0.45 cc/sec), 
a small change in the measured absorption rate will cause a large 
change in the calculated diffusion coefficient. This also means 
that, at flow rates of about 0.4 cc/sec a 2% error in the measure- 
ment of the hydrogen absorption rate will cause a 6% error in the 
value of D, and emphasizes the importance of having very good control 
of the experimental conditions. 

Use of the 3/4 inch diameter sphere was prompted by the fact 
that the resulting lower saturations of the outgoing liquid might 
provide better agreement with the absorption model by reducing 
possibilities of high errors occurring at very high saturations. 

As may be seen from Tables 2 and 3 the values of diffusivity obtained 
were always higher than those obtained using the | inch sphere. 

The reason for this is uncertain but may have been due to an end 
effect on the take-off rod. 

Deviations of the measured absorption rates from the theoretical 
model for liquid flow rates in excess of 0.45 cc/sec are probably 
due to convective mixing in the liquid film around the sphere. 
Although, at these flow rates, considerable activity was observed 
on the take-off level, the liquid film around the sphere appeared 
smooth. The condition of the take-off level as observed visually, 


was therefore taken to be a more reliable criteria for detecting 
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35. 
convective mixing than a similar observation of the liquid film 
around the sphere. 

The reported values of the take-off length used in previous 
investigations [3,4,7,9,33] with gas-water systems cover a fairly 
wide range. Lynn et al. [33] indicate that the allowable take- 
off length varies inversely as the liquid flow rate and report 
values ranging from 0.1 cm for high to 2.0 cms for low flow rates. 
Davidson and Cullen [3] took all their data at 1.5 cms though 
their data indicates that 1.0 cms would also have been adequate. 
Goettler [7,9] and Ratcliff and Reid [4] used values of 2.0 and 
4.0 cms respectively. Their data indicates no discernable 
dependence of the take-off length on liquid flow rate. It is 
apparent that each investigator repeated the original experiment 
of Lynn et al. to determine the take-off height best suited to 
their system. It is probable therefore that apart from the hydro- 
dynamics, system properties and general surface characteristics 
may also play a role in deciding the best take-off height for any 
one investigation. The value best suited for this investigation 
was 1.25 cms. 

When the stagnant layer is on the support rod the velocity 
profile on the rod may be assumed to be fully parabolic. This 
implies that the rate of gas absorption on the rod is reduced 
considerably. A plot of the gas absorption rate versus take-off 
height will thus show a decrease in slope when the stagnant layer 
reaches the junction between the bottom of the sphere and the 
support rod. If in addition, the surface area of the rod exposed 


for gas absorption is very much smaller than that of the sphere, 
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36. 
the slope of the plot will tend towards a limiting value of zero. 
In such an event using a value of take-off height somewhat larger 
than that for which the stagnant layer is at the bottom of the 
sphere would cause a negligible error in the measured gas 
absorption rate. The data of Goettler [7,9], Davidson and Cullen 
[3] and Ratcliff and Reid [4] apparently fall in this category. 
If, however, the surface area on the exposed rod is such that the 
gas absorbed on it could form a considerable part of the total, 
then the take-off level should be maintained such that the stagnant 
layer is at the bottom of the sphere. This condition would minimize 
errors in the measured gas absorption rate and thus in the calculated 
value of the diffusivity. 

Ratcliff and Reid [4] have reported an equation for the rate 
of gas absorption in a liquid flowing over a sphere assuming the 
velocity profile in the liquid film to be fully parabolic. The 
equation is of the form: 

— (3x1 68) aie arm (¢.-c.) (1-46) 
They assumed the low depth of penetration model of Davidson and 
Cullen [3]. The ratio n of the absorption rates given by equations 
(1-46) and (1-26) thus represents the reduction of mass transfer 
due to the stagnant layer. Since the equation for gas absorption 
over a sphere was derived [3] using the approach of Pigford [16] 
and Vyazovov [17] for gas absorption over a cylinder, the above 
ratio is identical for a cylinder of an equivalent radius. It 


may also be assumed that the ratio is the same for both the low 
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3/7 
and the high depth of penetration models. Table 4 gives the cal- 
culated values of the diffusivity assuming the exposed rod to be 
(i) completely inactive (n=0), (ii) completely active (n=1.0) 
and (iii) partially active (n=n). It is seen that under condition 
(iii) about 1% of the total measured absorption occurs on the 
support rod and this yields diffusivity values about 3.5% lower 
than those under condition (i). 

It was suggested earlier that the diffusivity values obtained 
using the 3/4 inch sphere were always high possibly because of an 
end effect on the support rod. Table 5 gives the calculated values 
of the diffusivity in a manner similar to that in Table 4 for the 
1 inch sphere. It is apparent that the end effect can be much 
more significant for the 3/4 inch sphere than for the 1 inch sphere. 
The diffusivities obtained under condition (iii) are within 3% of 
the corresponding values obtained using the | inch sphere. In the 
range of liquid flow rates investigated the average values of n 
obtained with the 1 inch and 3/4 inch spheres were 0.3024 and 0.3157 
respectively. The | inch sphere data reported in Figure 6 corresponds 
to the case where the support rod is considered inactive (i.e. n=0). 
The reliability of the data being within +5% it was considered 
unnecessary to incorporate condition CaL)e 

Moore [8] has indicated that the model of Davidson and Cullen 
[3] or Olbrich and Wild [15] is sensitive when the liquid saturation 
levels get high resulting in an increasing possibility of errors in 
the diffusion coefficient. A detailed derivation of the equations 
involved in developing the error analysis of the data obtained 


during this investigation is given in Appendix Cl. The final form 
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of the error equation is 
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where On» o and a. are the standard errors of the 


L7eG--Ce?. FC 
iE axe: S 
diffusivity, liquid flowrate, concentration driving force, gas 


absorption rate and viscosity respectively. Z is defined as 


- dig 79 


f(a) and a are as defined by equations (1-28a) and (1-29) and B(a) is 


given as 


(gp 3 Sa (1-476) 


It may be noticed that as the saturation level approaches 100%, 
F(a) + 1, Bla) + 0 and Z + © which implies that aD +> 6 Jat high 
saturations. Examining these equations it is also apparent that 
a eee as a" > 1 and econsequent yeas: 2 > Oe eInis simpli ecminat 


the saturation level of 100% is a limiting value at a limiting liquid 
flowrate of zero. Figure 7 is a plot of a vs. L over a liquid flow 


range of 0 to 0.65 cc/sec. The saturation levels have also been 


indicated relative to the liquid flow rate. It is apparent that 
= is very large at low values of L where the errors in diffusivity 
would indeed be large, but a decreases sharply with increasing 


liquid flow rate. In the range of flow rates corresponding to this 
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42. 
investigation (i.e. saturation levels of 80 to 85%) the value of 
o. is not very different from that at flow rates of 0.65 cc/sec. 
The reproducibility of the model in this region is therefore 


thought to be adequate. 


5.2 Experimental vs. Estimated Values 

Several semi-empirical correlations are available which are 
suitable for estimating the diffusivities of sparingly soluble gases 
in liquids. Figure 6 compares the measured diffusivities with values 
predicted using the Wilke-Chang [22] and Ferrell and Himmelblau [24] 
correlations. A comparison is also made with the experimental data 
of Haul and Puttbach [32] for hydrogen in ammonia. 

Moore [8] indicated that the correlation of Lusis and Ratcliff 
[23] was best suited for the data obtained with the hyd rogen-amino- 
methane system. The diffusivity data obtained by Moore was, however, 
preliminary and further work during this investigation indicated 
that the correlation of Ferrell and Himmelblau [24] is better suited 
for both the hydrogen aminomethane and hydrogen-ammonia systems. 
Despite the fact that the equation of Ferrell and Himmelblau was 
developed for water, the estimated values of the diffusivity of hydrogen 
in aminomethane using this equation were found to be within 12% of 
the experimental values obtained during this investigation. In 
addition, the deviation of the estimated values from the experimental 
data was approximately constant over the entire range of temperature 
investigated. Although, the average deviation obtained using the 
Lusis and Ratcliff correlation was 11%, the errors at higher tempera- 


tures were larger than those at lower temperatures. In other words, 
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43, 
a plot of diffusivity vs. temperature as shown in Figure 6 is not 
parallel to the mean line through the experimental data. For the 
hydrogen-ammonia system the equation of Ferrell and Himmelblau 
predicts values of diffusivity that are low, on an average, by 
about 9%, the deviation being about 12% at -10°C and about 4.5% 
at -40°C. In comparison, the equation of Lusis and Ratcliff 
predicts values that are low by about 19% through the entire range of 
-10°C to -40°C. The experimental data of Haul and Puttbach [32] was 
used for comparison. 

It is evident from Figure 6 that the correlation of Wilke- 
Chang [22] with an assumed association parameter of 1.0 is quite 
unsuitable for estimating the diffusivity of hydrogen in amino- 
methane. A similar situation is encountered for the hydrogen- 
ammonia system. 

The Wilke-Chang correlation employs an association parameter 
which accounts for the effective molecular weight of the solvent. 
In the absence of more reliable data, a value of 1.0 was used for 
the curve presented in Figure 6. In general, predictions obtained 
from the Wilke-Chang equation are poor when the solute gas is 
hydrogen. The experimental data for the hydrogen-aminomethane, 
hydrogen-water [24] and hydrogen-ammonia [32] systems were used 
to evaluate average apparent association parameters that may be 
used in the Wilke-Chang equation for predicting data for these 
and other similar systems. The average apparent association 
parameters are given in Table 6. Molar volumes for this correla- 
tion and the correlation of Lusis and Ratcliff were estimated by 


the method of Le Bas [34]. 
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TABLE 6 


Average Apparent Association Parameters for 
the Wilke-Chang Equation 


Average Apparent Reference of 
Association Parameter Experimental Data Used 
ae Ate) 


Ferrell & Himmelblau [24] 


This Work 


Haul and Puttbach [32] 


Haul and Puttbach [32] 
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It was pointed out earlier that Ferrell and Himmelblau [24] used 
the theory of absolute reaction rates to arrive at a generalised 
correlation that expressed the change of diffusivity with temperature. 


The generalised form of the equation was given as 








D = A i 
oF) a 
I 
ml 
C1 
Du, 
or : Sane (1-48) 
where A, is a constant for a particular solute-solvent pair. This is in 


contrast to the Stokes-Einstein [35] or the Eyring [26] models which 
Du 
suggest that — iS a constant for a particular size of the solute 


molecule and for a particular solvent. Ferrell and Himmelblau indi- 
Du 





cated that oh tends to decrease with an increase in temperature for 
the hydrogen-water and helium-water systems. The same fact is bee 
out a the hydrogen-aminomethane system. Average values of zd 

and at are given for different temperatures in Table 7. A word of 
caution is, however, in order with regards to the use of Table 7. The 


experimental data have been obtained between the temperature range of 
-20°C to -30°C but Table 7 includes values at -10°C and -40°C which were 
obtained by extrapolating the mean line of the D vs. 1/T curve shown in 
Figure 6. 

In order to take account of the presence of quantum effects in the 


diffusion of lighter gases dissolved in liquids, the quantum parameter 
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A* of Hildebrand et al. [25] was used to give the following equation: 





Sy) See A) temalal (1-44) 


Although Ferrell and Himmelblau [24] used this equation for predic- 
ting diffusivities of sparingly soluble gases in water, in principle 


there is a set of constants a and b for each solvent. 
Du l 


= were determined from the data of Haul and 





Average values of 
Puttbach [32] for hydrogen-ammonia and deuterium-ammonia systems along 


with the values for hydrogen-aminomethane obtained from the best fit line 

Loe 

———— | on Figure 8. The best 
Vo 

values of a and b determined from this plot are given as 


of Figure 6. These are plotted vs. | 


a 


as ore 


i} 


and b 


0.53 


The general equation (1-44) thus becomes 


DOS dh gee! - bee et (1-49) 


The molar volumes of the solute species used in this equation are as 
reported by Scott [36]. 

Equation (1-49) was used to predict the diffusivities of hydrogen, 
hydrogen deuteride and deuterium in aminomethane and ammonia. Represent- 
ative values of the diffusion coefficient using this equation are compared 


with available experimental data in Table 8. The predicted data are 
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FIGURE 8 


Estimation of Best Parameters for Equation (1-44) 
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50. 
within 5% of the experimental values. It is suggested that this equation 
be used to predict diffusivities of these systems in absence of experi- 


mental data. 


Bas) Hydrodynamics 


In order to be able to use the models of Davidson and Cullen [3] for 
evaluating the diffusivity of hydrogen in aminomethane, it is necessary 
that the flow over the sphere absorber be laminar. This is evident from 
the fact that both the low depth and the high depth of penetration models 
assume the film thickness around the ball to be represented by Nusselt's 
eduatlon OhutlOw.ovegea Llateplates| 12). 

Based on a Reynolds number criteria the data of Davidson and Cullen 
[3] and that of Goettler and Pigford [9] indicate that the condition of 
ripple free flow was maintained for values of Re in excess of 110. In 
comparison with this, the experimental data of a number of other authors 
[37-42] have indicated that wave inception begins on the liquid at 
extremely low flow rates well within the viscous flow region. Kapitza 
[43] has deduced that the Reynolds number at wave inception is given by 


1/11 
Re = 2.4 Seely (1-50) 


and has shown that in the pseudo laminar region for Reynolds numbers 
Re > Re, the velocity profile in a liquid film flowing over a vertical 
plate remains parabolic (as derived by Nusselt), but the mean film thick- 


ness is given by 


2.4 vr 173 
Vv 


A = is gow (1-51) 
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The work of Tailby and Portalski [44-48] has contributed to a better 
understanding of the flow of liquid films. Portalski defines the follow- 
ing flow regimes with increasing liquid flow rate [47]: 

(i) Steady Laminar (no waves or ripples) 

(ii) Pseudo Laminar (wavy flow) 

(iii) Transition 

(iv) Pseudo Turbulent 

(v) Turbulent 
(i) and (ii) constitute what has been known as the viscous or laminar flow 
regime, and (iv) and (v) constitute the turbulent regime. Portalski [47] 
suggests that Nusselt's equation applies only in the steady laminar region 
but not in the pseudo laminar regime where Kapitza's equation represents 
the experimental data better. Beyond the pseudo laminar region even 
Kapitza's equation deviates from experimental data considerably. He also 
indicates that the physical properties of the system, namely specific 
gravity, viscosity and surface tension affect the film thickness and 
suggests that in the pseudo laminar region as the surface tension becomes 
smaller the deviation between the values predicted by Nusselt's and 
- Kapitza's equation becomes smaller. This implies that the use of 
Nusselt's equation in the pseudo laminar region for a system such as 
aminomethane could well be justified. 

In order to establish the above regimes firmly, however, it is nec” 
essary that more work be done with systems where effects of viscosity, 
specific gravity and surface tension can be independently investigated 
through actual measurements of the liquid film thickness. 


For the flow of liquid aminomethane over the sphere absorber it was 
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52. 
found that the maximum Reynolds number that could be attained without 


any rippling as seen through the naked eye corresponded to about 50. 
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CHAPTER 6 


CONCLUSIONS 

The single sphere absorber has been used to obtain reasonable values 
of the molecular diffusion coefficient for hydrogen in aminomethane. The 
accuracy of the reported diffusivities is dependent upon the careful 
control of experimental conditions, on the accuracy of solubility and 
viscosity data, and on the applicability of the hydrodynamic model. It 
is estimated to be within the +5% range reported by others who have used 
the sphere absorber. 

Agreement with the theoretical model of Davidson and Cullen [3] was 
obtained only for liquid flow rates where the depth of penetration of the 
solute gas into the liquid film was high. Deviations from the model are 
thought to be due to convective mixing. Surface tension may be an important 
factor in determining allowable flow rates, but additional work is 
required to establish criteria useful for predicting the onset of con- 
vective mixing in films over spherical surfaces. 

The method suggested by Ferrell and Himmelbiad [24] has been used to 
correlate experimental data for the diffusivities of hydrogen in amino- 
methane and ammonia. Equation (1-49) is recommended for use in esti- 
mating the diffusion coefficients for hydrogen, hydrogen deuteride and 


deuterium in aliphatic amines. 
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NOMENCLATURE 
A interfacial area aaah 
A, corrected interfacial area lo 
A. surface area of a dry sphere fom>) 
Ay a constant for a given solute-solvent pair defined by 
equations (1-42) and (1-44) | eat | 
sec K 
A constant in Erying's equation (1-37) 
A, constant of integration in equation (1-41) 
a association parameter in equation (1-33) 
a constant in equation (1-44) 
b exponent in equation (1-44) 
C concentration of solute (SS) 
Cc. concentration of solute Mette interface a) 
Co concentration of solute in the bulk liquid (MS) 
Ch concentration of solute in the bulk liquid are the nozzle ac 
D diffusion coefficient of the dissolved solute — a 
d diameter of take-off tube (cm) 
d; diameter of the jet (cm) 
E given as the function defined in equation (1-45) 
ED energy of activation for diffusion (852) 
E energy of activation for viscosity (ee) 
G rate of absorption of the solute gas in the liquid eS 


G(1 atm) rate of absorption of solute gas corrected to 1 atm solute 


‘ gms 
partial pressure ome 
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rate of absorption of solute gas in the liquid jet (=) 


rate of absorption of solute gas in liquid over the sphere oad 


rate of absorption of solute gas in liquid over the sphere 


after area correction due to surface film Ge 


rate of absorption of solute gas due to a temperature 


ms 
gis) 


fluctuation ( 
sec 


rate of absorption of solute gas due to a take-off level 


fluctuation (=—— 
Sec 





acceleration due to gravity | et | 
\ =e atm 
Henry's Law coefficient SE ESERIGA. 
Planck's constant (erg-sec) 
height or length of jet (cm) 
take-off height (cm) 
pets cm 
liquid flow rate aa) 
molecular weight of solvent (—_) 
gm-mole 


molecular mass (gms) 


atoms 
Se 


Avogadro's number ( ayia 


number of moles 

total number of terms in equation (1-28) 
vapour pressure of the solvent (atm) 
partial pressure of solute gas (atm) 
total system pressure (atm) 


3 


volumetric flow rate of solute gas () 


radius of sphere (cm) 


aR - 
(52) 401 bivptt ong oi zap s2uloe Yo noisqvoeds Toveaes 
(FRB) sverige eds isvo biupi!l ni eap Stuloe to no! tqueeds To 6287 
sisrige sz 1sve bivpil nl eep stufoe to noizquoeds % S967 

(see) mti¥ sostiue of aub noliger0s 6816. 199% 
siute1sq9met 6 62 sub esp stuloae to no} 3qvezds to sie, 

(288) noizeus oul? 
{eval tio-sds3 6 09 eub 269 siulo2? to noltqveeés To oga7 


(=) noi teusoul? 








e | Viiveie oF sub nolisielsacs 


m6 


(—saeatreton insivititsos wed e*yiAek 
noljos’ 

(ose-py9) jinssenoo 2@*Hons!4 

(m>) te{ to dtens! yo tdgted 


(ma) -thpied tto~ades 





t 
<=) ste wor? biupit 
(See insyfoz Jo selistane 16{uoalom 
(emp) 2esm 16luoslom 
(Sere ) sedmun 2! o vbspevA 
2g1Gm, to, 19dmun 
(88-{) noissups ni 2mis) to isdmun 16303 
ute ssi lye 2702 peer 





elites 





56. 


cal 


Rx gas constant rcnoice 
Re Reynolds number 
Re, Reynolds number at wave inception 
1, temperature (°K) 
Ty absorber temperature (°C) 
Ty hydrogen soap-meter temperature (°C) 
t contact time (sec) 
V volume of absorption system (cm?) 
vA volume of absorption system excluding take-off volume (cm?) 
3 

V molar volume of solvent (—“~——) 

] gm-mole 
V, molar volume of solute (—““—) 

gm-mole 
f par: 4 cm 

Vy velocity of liquid along a stream line ee 
Vo velocity of liquid at the surface of the film a 
x distance measured perpendicular to the stream lines 


from the surface of the film inwards (cm) 


Xo mole fraction of the solute gas in the liquid 
y dimensionless distance => 
Z function defined in equation (1-47a) 


Greek Symbols 


A thickness of the liquid film around the sphere (cm) 


Re thickness of the liquid film at the equator of the 


sphere (cm) 


: ree cm 
v kinematic viscosity of aminomethane liquid aa) 
6 the angle between the vertical and a radial line (radians) 


wv a stream function as defined by equation (1-16) 
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p a transformation variable defined by equation (1-19) 
5 the limit of ¢ for a sphere; defined by equation (1-25) 
oO a transformation variable defined by equation (1-29) 
Oy a property of the solute defined by equation (1-36) 
Oy standard error in diffusivity 
by Standard error in liquid flow rate 
Ae o¢ standard error in driving force 
io 
Oe standard error in gas absorption rate 
S 
at, standard error in kinematic viscosity 
B a transformation variable defined by equation (1-32a) (eK!) 
B pre-exponential coefficients of equation (1-28) 
Y; exponential coefficients of equation (1-28) 
u viscosity (centipoise) 
My viscosity of solvent (centipoise) 
Ax quantum parameter defined by equation (1-43) 
Oo collision diameter of the solute molecules (A) 
o; surface tension of the solvent (dynes/cm) 
2 
€ collision energy of a solute molecule —s- 
sec 
n ratio of absorption rates as defined by equations (1-46) 
and (1-26) 
Ms cm 
6 kinematic surface tension of the solvent = — Sei 
p 
sec 
0 density of the solvent (2) 
emp cm? 
Ne liquid volumetric rate per unit length eae 
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APPENDIX Al 
Experimental Data 


The experimental data on the physical rate of absorption of hydrogen 
in aminomethane are summarized in Tables Al-1 and Al-2. 
Only those runs where the experimental control of the variables was 


satisfactory have been included. 


The significance of the variables listed in the tables is indicated 


in the nomenclature. 
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APPENDIX Bl 


The following start-up and operating procedure may be used for the 


measurement of diffusivity of hydrogen in liquid aminomethane using the 


Single sphere absorber. 


(1) 


(2) 


(3) 


(4) 


Charge the system with hydrogen and make a few soap bubbles in the 
Soap-meter to check for possible gas leaks in the system. 

When the soap bubbles are stationary for about 30 minutes, charge 
the hydrogen saturator with pure aminomethane. 

Open the solenoid valve Sy between the hydrogen saturator and the 
gas absorption chamber. 

Start the peristaltic (hydrogen circulation) pump and circulate 
hydrogen for about one hour. The system temperature will rise 
during this time. 

Stop the peristaltic pump, close solenoid valve Sy and allow system 
temperature to stabilize. 

Immerse the aminomethane dump tank in liquid nitrogen. 

Open the solenoid valve Sy and the solenoid valve Sy between the 
liquid distributor and the absorption chamber. 

Introduce aminomethane to the system by opening the micrometer 
control feed valve. Maintain the flow rate at a level higher 
than the discharge rate of the nozzle. This will cause amino- 
methane to accumulate in the distributor. 

Step 8 will cause the system pressure to increase substantially 
due to aminomethane vapour pressure. Vent the system periodi- 


cally to ensure that the pressure remains above the setting of 
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65. 
the gas regulator. 

When the aminomethane level appears in the take-off tube, open the 
valves in the liquid return line and adjust the micrometer 
control liquid return valve so as to maintain the liquid level 
in the take-off tube. 

When the level of aminomethane in the liquid distributor reaches 
the glass tube above the distributor, close the solenoid valves 
Sy and Sy and adjust the feed control valve to maintain the 
desired flow rate as seen on the liquid rotameter. 

Bring the pressure of the system to the value corresponding to the 
regulator setting by venting the excess pressure from the system. 

Open solenoid valves Sy and Sy momentarily to balance system 
pressures and then close them again. 

Wet the soap-meter wall with soap solution by allowing some hydrogen 
gas to bubble through it and go through the soap-meter. 

Carry out the final adjustments on the aminomethane feed rate and 
liquid level in the take-off tube. 

Open the hydrogen feed valve and allow the system to stabilize. 

Record the readings of the four system thermocouples during this 
transient period. A steady temperature record indicates that 
the system is reaching stability. 

After reaching a steady state with regard to the temperatures, form 
two soap films in the soap-meter and record the time taken by the 
bottom bubble to sweep a given volume (5.0 cc). 

Control the liquid level in the take-off tube at the desired level 


for the entire duration of the run. A cathetometer should be 
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used for this purpose. 


The liquid flow rate or take-off height may be changed during a 


run, but at least 20 minutes should be allowed for the system 


to stabilize after any change is made in the absorber conditions. 


Upon completion of a run, close the hydrogen feed valve and drain 


the liquid aminomethane from the absorber flow lines and the 
liquid saturator into the dump tank. The remaining aminomethane 
in the feed tanks should also be drawn into the dump tank. 
Isolate the dump tank from the absorber and evacuate the dump tank 
to a pressure of about 5 microns. Ascertain that sufficient 
liquid nitrogen is present around the dump tank. 
Isolate the dump tank from the vacuum system and allow it to warm 


up to room temperature. 


Evacuate the gas-absorption system so as to clear it of any 


residual aminomethane and hydrogen. 


On the day following the run, immerse the feed tanks in liquid nitro- 


gen and evacuate them. lsolate the vacuum system and feed 
aminomethane, under its own vapour pressure, from the dump tank 


to the feed tanks. 


(26) When all the aminomethane is transferred, isolate the dump tank and 


evacuate the feed tanks so as to degas the aminomethane further. 
A pressure of about 1 micron should be maintained in the feed tanks 
for about four or five hours. This requires that the liquid 


nitrogen level around the feed tanks be kept high during this time. 


(27) Isolate the feed tanks and allow to warm up to room temperature. A 


period of about forty hours should be allowed for this purpose. 
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APPENDIX Cl 
The equations used for the analysis of the errors in the data 
obtained for the determination of the diffusivity of hydrogen in amino- 
methane may be developed as follows. 
The error of a number of measured quantities may be expressed in 


terms of a standard error defined as 


k= ail = 
a . (C1-1) 
where 
Oo = standard deviation of the measured quantity 
and N = the total number of observations. 
For any function 
Gee=e ft (ASB, Gad) 
The standard error Oe in G is given by 
Z Z 2 2. 
2a of 2 of 2 of 2 of. Zz # 
Seo Meme Ge Gye) ie er Ie Ue rries Sy (C152) 
where G is the mean of G and nr ps Aes Ay are the standard errors in the 


‘means of A, B, C and D. 

The rate of gas absorption on the sphere absorber is represented by 
the high depth of penetration model of Davidson and Cullen [3].  Neg- 
lecting the last two exponential terms in f(a) the equation may be written 


as 


@ 
i] 


L(c.-C.) [1 - 0.7857 exp (-3.414a) - 0.1001 exp(-26.21a)]  (C1-3) 


ea es} 
Wierendme=12 xl 60m [Sy “73 D (1-29) 
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68. 
The contribution of the two exponential terms neglected in equation (C1-3) 
is negligible at the saturation levels involved in this study. Equation 


(C1-3) may be expressed in general as 


@ 
Il 


f,[L, (C.-C), D, »! (c1-4) 


3G 2 9 
+oie——— | o (C1-5) 
av Vv 
Representing 


ia) = 42-70. /85/ exp (-3.4140) - 0.1001 exp (-26.210) (C125) 


3G dG dG dG 
the derivatives — 2 — and — may be represented 
aby a(C.-C.) > aD av 


as follows: 


G 
a = (C.-C\) F(a) + Ec ec.) [-0.7857(-3.414) exp (-3.414a) 
ele} 
- 0.1001 (-26.21) exp(-26.21a)] a0 
8G 
er = (lease) HEN L(c,-C.) [2.6824 exp (-3.414a) + 2.6234 
ele} 
exp(-26.21a)] <> (C1-7) 
Representing 
B(a) = 2.6824 exp(-3.414a) + 2.6234 exp (-26.21a) (c1-8) 


equation (C1-7) reduces to 
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The remaining derivatives are obtained in a similar manner: 


sb L(C.-C) B (a) 


L f(a) 


L(C.-C.) B (a) 


The partial derivatives of ‘'a' 
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may be obtained from equation 
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(C1-9) 


(Ci-10) 
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(G1e12) 


(1-29) as 


(C1-13) 


(C1-14) 


(C1-15) 


Substituting from equations (C1-9) to (C1-15) in equation (C1-5), Og may 


be written as 
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Considering the fractional errors to be all additive (thus representing 


the worst situation), a, may be expressed as 
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The total measured hydrogen absorption rate in gms/sec into the 


absorber may be written as: 
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ep hydrogen soap-meter temperature (sec) 
Qu, = volumetric flow rate of hydrogen (cc/sec) 


Following a similar pattern of analysis as shown for Oe» the following 
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equation may be written: 
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The total measured hydrogen absorption rate may be considered to 


consist of the following: 
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where 
Gy = total measured gas absorption rate a 
G. = actual absorption rate on sphere (ae) 
A sec 
Gye das absorbed on jet ) 
ii sec 
Gr = hydrogen absorption rate due to a change in absorber 
temperature during a run (a=) 
sec 
Gre = hydrogen absorption rate due to a change in the take-off 


level during a run (2) 
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G. and Gr are assumed to have a mean of zero so that only the expected 
error due to their presence is considered. From equation (C1-25) 


the fractional error may be expressed as 
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The effect of a change in the absorber temperature and the liquid 
take-off height may be determined from the equation describing the 


hydrogen partial pressure in the absorber. 


Po PL cs P, (C1-28) 
where Po = hydrogen partial pressure 
PL = total system pressure 
P, = aminomethane vapour pressure 
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Assuming the gas law to hold 
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Since the effect of Gt and Ge is to change the number of moles at 


constant total pressure. 


From equation (C1-30) 
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The change in absorber volume due to a change in the liquid take-off 


height may be evaluated as follows: 
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where d = diameter of the take-off tube 
hoy = take-off height 
V = volume of the absorption chamber excluding the take-off tube 
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The actual hydrogen absorption rate over the sphere G, is related to the 
rate G. through the area correction accounting for the film thickness over 


the sphere. 
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Using the definition of ACOR (equation 1-30) and equation (1-31), the 


following equation may also be written down: 
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The bulk concentration Co of hydrogen in aminomethane flowing over 


the sphere is given by 
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where My = molecular weight of aminomethane 


density of aminomethane 
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D partial pressure of hydrogen 
Following the procedure outlined earlier, the fractional error in the 


interfacial concentration is given by 
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The derivatives of Py» Py and H, with respect to temperature are obtained 


from the defining relations of Py> P, and H, asvant UnGt ono Ths The 


relations for Pi and H, are given as equations (C1-39) and (1-32) respec- 
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tively. The variation of density with temperature is given as 
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Using the relationships derived for the different variables used during 
an experiment (equations (C1-20) to (C1-50)) the error in the diffusion 


coefficient may be explicitly evaluated using equation (C1-19). 
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CHAPTER | 


INTRODUCTION 

The kinetics of the hydrogen-aminomethane exchange reaction have 
received limited attention in the past [1,2] and the data available do 
not conform to the operating conditions of interest to industry. 
Industrial contactors are expected to operate with a high concentration 
of the homogeneous potassium methyl amide catalyst and within a 
temperature range of -40°C to 70°C. The data of Bar-Eli and Klein 
[1] was restricted to very low lithium catalyst concentrations of 
less than 0.3 gms Li/kg solution and temperatures of less than -22°C, 
while Rochard [2] obtained her data at potassium methyl amide concen- 
trations of less than 5 gms K/kg solution and temperatures not 
exceeding -62°C. The above mentioned investigations were both conducted 
in stirred type reactors and difficulties may have been encountered 
in eliminating mass transfer resistance at higher temperatures and 
catalyst concentrations. 

The work reported here uses the single sphere absorber which allows 
the measurement of the rate constant at conditions of greater commercial 
significance. 

The exchange of isotopic hydrogen with aminomethane in the region 
of low deuterium concentrations can be described by the following 


reversible reaction: 


k 
| 
HD + CH.NH +H - 
D+ CHaNH) © CH,NHD + H, (2-1) 
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which takes place in the liquid phase and is catalyzed by potassium 
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methyl amide. Bar-Eli and Klein [1] and Rochard [2] investigated the 
kinetics of the above reaction and concluded that the exchange process 
is of first order with respect to HD and CHNHD concentrations. They 
indicated that the reaction rate is substantially higher than for the 
analogous hydrogen-ammonia system. Since the latter reaction is known 
to be fast [3] it is reasonable to suggest that the hydrogen amine 
system also involves a fast reaction. 


Assuming that the concentrations of CH and H, are large and 


Be 2 
substantially invariant, the exchange process may be represented by 
a pseudo-first order reaction given as 
HD : CH,NHD (2-2) 


3 
Ko 


The equilibrium constant for the reaction (2-2) may be written down as 


; [CHNHD], ky 


gr es = (2-3) 
[HD], k, 
where [CH,NHD] . and [HD] . represent the concentrations of the deuterated 


species at equilibrium. The local rate of reaction 'r' may be expressed as 
r= k, [HD] = k, [CHNHD] (2-4) 


which on substitution from (2-3) gives 


sy 
r= k, [HD] -— [CH NHD] 
[CHNHD] 
= k,({HD] - —+— ) (2-5) 
K 


Equation (2-5) could also be written as 


r= k, (Cc - c.) (2-5a) 
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where 'C' is synonymous with [HD] and C. is a concentration of HD that 
would be in equilibrium with [CHNHD] . If the bulk volume of the 
liquid aminomethane used is such that the molar ratio of the liquid 


to hydrogen gas is large then the product concentration [CH.,NHD] 


B 
is essentially constant and C. in equation (2-5a) would also be a 
constant value. 

The rate constant, Ky» of the exchange reaction can be simply 
measured by knowing the total liquid holdup and the absorption data, 
provided the resistance for transfer of deuterium in the liquid phase 
can be sufficiently reduced so that the reaction becomes the rate 
controlling process. This condition has been defined as the kinetic 
regime by Astarita [4]. When, however, the reaction is a fast one it 
becomes exceedingly difficult to reduce the mass transfer resistance 
sufficiently to operate in the kinetic regime. The data of Bar-Eli 
and Klein [1] and Rochard [2] was obtained in stirred type reactors which 
were provided with mechanical vibrators and stirrers so that mass transfer 
resistance may be reduced. In spite of this, their work was restricted 
to low catalyst concentrations and low temperatures because of diffi- 
culty in operating in the kinetic regime. 

In contrast to operation in the kinetic regime there are a number 
of contactors which make it possible to conduct the exchange reaction 
in the fast reaction regime [4]. These include the wetted wall column, 
the laminar jet, the multiple sphere absorber, the single sphere 
absorber, the disc column absorber and the rotating drum absorber. 

The fast reaction regime is typified by high diffusion times as com- 


pared to the reaction times and a narrow reaction zone close to the gas- 
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liquid interface. For the pseudo-first order reaction (2-2) operating 


in this regime, the overall chemical absorption coefficient K, is given 


i 
by 


Kee vice D (2-6) 


where Ky is the pseudo-first order rate constant and D is the molecular 
diffusion coefficient of the solute in the liquid phase. The rate of 


absorption of the solute specie is represented as 
5 say) (2-7) 


where C, and C. are the interfacial and chemical equilibrium concentra- 
tions of the solute in the liquid phase. If K is known through experi- 
mental determination, then using the values of D estimated in the first 
part of this thesis, it is possible to evaluate the kinetic constant 
k, from equation (2-6). 

In preference to the other contactors mentioned above, the single 
sphere absorber was used for this study. Its merits over the other 


absorbers may be summarized as follows: 


(i) [lt minimizes the end effects common to other gas liquid 
contactors. 
bi) lt provides longer contact times (0.1 sec to 1.0 sec) which 


is desirable for operation in the fast reaction regime. 

(iii) Relatively high liquid flow rates without surface rippling 
may be employed with this absorber. 

(iv) Its surface area exposed for absorption is large which makes 


it suitable for absorption of sparingly soluble gases such as hydrogen. 
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CHAPTER 2 


LITERATURE SURVEY 
A study of the kinetics of the exchange of hydrogen with ammonia 
Or aminomethane involves an investigation of the manner in which the 
rate of exchange varies with the concentrations of hydrogen and of 
the catalyst in solution. The work of previous investigators [5,6,7,8] 


in this area is outlined below. 


2.1 Hydrogen-Ammonia System 

Aminomethane is a simple substituted derivative of ammonia and as 
such, exchange studies of the hydrogen-ammonia system with alkali metal 
amides as catalysts are pertinent and useful towards an assessment of 
the behavior of the hydrogen-amfnomethane system. 

Wilmarth and Dayton [5] used a stirred reactor to obtain exchange 
data for the hydrogen-ammonia system at -53°C using potassium amide 
as a catalyst in solution. Their work was restricted to very low 
catalyst concentrations ranging from 0 to 10 millimoles/litre (mm/1) . 
The rate of exchange was reported to be proportional to the concentra- 
tions of hydrogen and the amide ion. A rate constant Ky = thas} nine 
was obtained at a potassium amide concentration of 2.7 mm/]. Bar-Eli 
and Klein [6] also used a stirred reaction cell and extended the 
hydrogen-ammonia study to higher catalyst concentrations ranging from 
5-200 mm/1. They covered a cempenertre range of -80°C to 0°C 
but their data of k vs. temperature appears to be inconsistent 


when compared with the data reported for k vs. catalyst concentra- 


tion. Based on the data of k vs. temperature, however, they 
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indicate that an activation energy of 7.5 kcal/mole was obtained for 
the bvdtcaen annonis exchange reaction. Dirian et al. [7] investigated 
a fairly wide range of potassium amide catalyst concentrations extending 
from 0.4 to 1400 mm/1. They indicate that the mass transfer resistance 
in the liquid phase could not be eliminated at conditions of high 
catalyst concentrations and high temperatures inspite of stirrer 

speeds of about 9000 rpm. At -50°C the mass transfer resistance could 
be eliminated only for catalyst concentrations upto 60 mm/1. From a 
plot of the kinetic rate constant vs. the catalyst concentration at 
-50°C reported by them, it appears that kK = 1.5 Mn at ai aur 2.5 mm/] 
and kK = 2.5 ion at about 5.0 mm/1l. This would indicate concurrence 
with the work of Wilmarth and Dayton [5]. An activation energy of 8.0 
kcal/mole was reported for the reaction. Delmas et al. [8] also 
conducted a study of the exchange of the hydrogen-ammonia system using 

a stirred reactor. Their work was directed towards the evaluation of 
the dissociation coefficients of the potassium amide catalyst in solution 
and their conclusions have helped to identify the catalytically active 
specie. The study covered an investigation of amides of potassium, 
rubidium and cesium within a concentration range of 1 to 200 mm/1 at 
-45.3°C, -59.3°C and -70°C. Typical values of the rate PG at 
-59.3°C and at potassium amide concentrations of about 5 mm/1 and 150 
mm/1 are approximately 1.0 ian and 20.0 muni, respectively. A value 
of 5.5 kcal/mole is reported for the activation energy of the reaction 
and the discrepancy with the value of 7.5 kcal/mole obtained by Bar-Eli 
and Klein [6] is attributed to different values used for the solubility 


of hydrogen. The work of Bourke and Lee [3] appears to be the only one 
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where the:exchange reaction was conducted with contactors which allowed 
accounting for the diffusional resistance in the liquid phase. They 
used a wetted rod column and a disc column to obtain the dependence 
of the rate constant with temperature and catalyst concentration. 
They reaffirmed the first order dependence of the rate of reaction on 
the hydrogen concentration and indicated that the reaction zone in the 
liquid was very narrow, with the bulk of the reaction occurring within 
a depth of about 0.001 in. below the gas-liquid interface. Rate constants 
were found to lie in the range 10-500 oe at potassium amide catalyst 
concentrations of 4-40 gms KNH ,/1 (73-730 mm/1). The temperature 
range investigated was -40°C to 20°C, and an activation energy of 8.0 
kcal/mole was reported. The above authors, however, used diffusivity 
values for hydrogen-deuteride in ammonia as predicted by the Wilke- 
Chang equation using an association parameter of unity. This equation, 
as indicated in the first section of this investigation, underestimates 
the diffusion coefficients for the hydrogen-ammonia system by as much 
as 100%. The data of Bourke and Lee [3] when corrected for diffusivity 
using equation (1-49) gave values of the rate constant in the range 
10-160 ee for the same temperatures and catalyst concentrations. The 
activation energy corresponding to the corrected data was found to be 
5.8 kcal/mole. A plot of the corrected values of the rate constant Kk, 
vs catalyst concentration at 20°C indicated that kK, varied linearly 
with catalyst concentration upto about 365 mm/1. 

Values of the rate constant at two temperatures and catalyst 
concentrations as obtained by some of the authors mentioned above [3,6, 
7,8] are compared in Table 1. In order to obtain a common basis for 


comparison the data had to be extrapolated linearly both for tempera- 
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ture and catalyst concentration wherever necessary. It may be noted 
from Table | that there is a reasonable agreement between the values 
obtained from the data of all the investigators at 163.7 mm/1 and 
-59.3°C whereas at 655 mm/1 and -23°C the value obtained from the 

data of Delmas [8] is significantly lower than that obtained from the 
data of the other authors [3,6,7]. This appears to be due to the 

fact that a linear extrapolation in catalyst concentration upto 655 
mm/1 underestimates the rate constant obtained from the data of Delmas 
[8] and this cannot be compensated for through a temperature extrapola- 
tion using an activation energy of 5.5 kcal/mole. The higher activation 
energies reported by Bar-Eli and Klein [6] and Dirian [7] seem to com- 
pensate for the underestimation and result in a better agreement with 
the measured value reported by Bourke and Lee [3]. 

201 AMPEfTectwor Catalyst Concentration? In order to understand 
the manner in which the catalyst concentration affects the rate constant 
K)> it is desirable to have some information regarding the catalytically 
active specie. Wilmarth and Dayton [5] reported that the rate of 
exchange of the hydrogen-ammonia reaction is first order with respect 
to the amide ion concentration [NH] which implies that the amide ion 
is the catalytically active specie. The rate constant could thus be 


expressed as 


k, = k [NH] (2-8) 


The authors suggested that the dissociation of potassium amide in solution 


occurs according to the equation 


> - a ‘n 
KNH, 2K’ + NH, (2-9) 
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and that the reaction proceeds through the formation of a hydride ion 
ine. 


a ue = 
NH. EHO eaees NH.D +H 


(2-10) 
H+ NH, > NH, + Ho 
Bar-Eli and Klein [6], however, suggest that in addition to the amide 
ion, the undissociated potassium amide also catalyzes the reaction 
although its activity is 20 times lower than that of the amide ion. 
Their data indicates that the reaction rate is linear with catalyst 


concentration through the entire range of concentrations investigated 


(5-200 mm/1). The rate constant k, was expressed as 
= i} 
KS k [NH] wake [KNH.] (2-11) 


Dirian and his coworkers [7] reported the catalytic effect of the alkaline 
cation and suggested the following equation for the rate constant: 


kK, = Kol NH + Kol °NHS at (2-12) 


where NHS and ay,+ are the activities of the anion and the cation 
species respectively. The rate of exchange was found to depend linearly 
on catalyst concentration only up to a value of 20 mm/1 after which the 
first order dependence was found to be invalid. 

The apparent uncertainty about the catalytically active specie 
was due to the lack of information regarding the dissociation coefficients 
of amides and the activity coefficients of ions in liquid ammonia. 
Delmas et al. [8] conducted a kinetic study of the hydrogen-ammonia 
exchange reaction and carefully calculated the dissociation coefficients 


for amides of potassium, rubidium, and cesium. They conclude that the 
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TABLE | 


Comparison of the Reaction Rate 
Constant for the Hydrogen-Ammonia System 


Catalyst Activation 
Reference Concentration Temperature Rate Constant Energy 
mm/ | } kcal/mole 


Bar-Eli & 
Klein 


[6] 


Dirian 
et al. 


[7] 


Delmas 
et al. 


[8] 


Bourke 
& Lee 


[3] 





* Indicates that the variable was extrapolated outside the reported range. 
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/ 


free amide ion is the only catalytically active specie and the rate of 
exchange is linear with the concentration of the amide ion of potassium 
upto a value of 5 mm/1. This corresponds to a potassium amide overall 
concentration of about 65 mm/1. They suggest that the deviation from 
the first order law above these concentrations is due to the presence 


of triple ions like (K,NH) 


2.2 Hydrogen-Aminomethane System 
(a) System Properties 


The investigation of the hydrogen-aminomethane exchange 
system is of recent origin and as such the extent of data available on 
fundamental properties of the system is very limited. The properties 
of particular interest from the point of view of this investigation 
include: 

(i) The hydrogen solubility. 

(ii) The diffusivity of hydrogen-deuteride in aminomethane. 

(iii) The equilibrium constant of the exchange reaction. 

Rochard [2] has reported data for the solubility of hydrogen in 
aminomethane at three specific temperature namely -62.2°C, -77.2°C and 
-90°C. These were the temperatures at which her exchange studies were 
conducted. Moore [9] undertook a study to measure the solubilities of 
hydrogen in amines over an extensive temperature range. His work for 
aminomethane covers a temperature range of -62°C to +23°C. Assuming 
Henry's Law to apply the temperature dependence of the Henry's Law 


coefficient is reported as 


In H, = 4.1403 + 1.66128 - 0.116087 (2-13) 


oe Ws ! — 
a 
| ee 


. 88 





to ote, adj bas sisaqe svitoés Vileoigyletso Vine BM 2) fol Shims e577 
mizz26309 Yo noi sbime say to noiyenznesaes os) hei, veeni! 21 spraioxe ; ‘ 


ifersvo sbime mileeetog 6 of sbnayes too 2/AT !\n 2 Yo sulev. 6 oxqu 


mot? noltel¥eb ond feds tesoeu2 yonT . t\n 20 Segs ae nos 1813 I8QNOD | } 
iis 
saneesig st 01 sub 2i his | Senhinniics Ged svods wel rebto y2317 orld 
| YH) St anos StQIit Ye" 
Meteyé Shatgenon ipa $.§ 
#2 |d7aqoa8 syeNs (s) 
spréisks snediomonime=aspolbye afd to noltepiteovns ot 
no sidelisve stab to Inagxs ond roue es bie mibise Gneset Yo 21 mazaye 
28is2eqoig sit .bszimil yvav el meseye en to esisysqow | sinemsbdauT 
noiIgce ein cidj} 3a wely to tnidq sd} mov T2etednl rsluotd18q to 
sebuloni 
vil lidulee neporbyd ed? (i) 
siedismionine nl sbivasush-naposbyd to yiivieuttib et (71) 
Jnoljoses sonciioxe sf3.'70 ané2znoo mutrd! fidpe, sat (iii) a a 
i} ngpoibyA Yo yxi!l iduloe erly 107 eieb betnoqe 26d ($] bretoon 
bans 2°S.5f- .2°S.Sd- yisman saujeieons? oi tiseqe soind tq snerizemonime 
g1aw 29ibuse ogiétione 19d piity 36 -eonuaeisqnes Sree aot + 3°0R- ; 


ee Ftdain pi eer He! 





89. 


where B = TORK) (2-14) 
and P, = oH. (2-15) 


Po is the partial pressure of hydrogen in atmospheres and Xo is the mole 
fraction of hydrogen in aminomethane. Equation (2-13) is used for esti- 
mation of the hydrogen deuteride solubility for this investigation. 

The data obtained in the first part of this study is perhaps the 
only data available for the diffusivity of hydrogen in aminomethane. 
The diffusivity of hydrogen deuteride is estimated using the modified 
Ferrell and Himmelblau equation reported as equation (1-49) in this 
thesis. 

The equilibrium constant for the hydrogen-aminomethane exchange 


reaction (2-1) may be written as 


[CHNHD] [H.] 


K = (2-16) 
eq [CH NHI THD] 


where [CHNHD] , [H.] etc. represent the concentrations of the species at 
chemical equilibrium. The separation factor relating the deuterium 


_ concentrations in the amine and hydrogen at equilibrium is expressed as 


DN per 
a. = D/H (2-17) 
hydrogen 


Rearranging equation (2-16) one obtains 


[CHNHD] ™ (2) 
[CHNH, ] HUN e 
Keg . ‘ 2 ds : amine (2-18) 
HD 
[Ho] mL 


hydrogen 
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or K = ~ 0, (2-19) 


where m* and n* are the number of exchangeable hydrogen atoms per mole- 
cule of the amine and hydrogen respectively. Since aminomethane and 
hydrogen both have two exchangeable atoms per molecule m* is equal to 


n* for the hydrogen-aminomethane system. This implies that 


os3 = Obe (2-20) 


AECL [10] have conducted some studies to measure the separation 
factor for the hydrogen aminomethane system. Their present conclusions 
indicate that the separation factor for this system may be taken to be 
equal to that for the hydrogen-ammonia system. The expression for the 
change of the separation factor with temperature for the hydrogen- 


ammonia system is given as 


te eee 237059 . 
1319 aes 0.2422 + aaa (2-21) 
where T is the temperature in °K. In the absence of indications to the 


contrary oer is taken equal to oer for this investigation also. 
(b) Exchange Studies 

Kenyon and Pepper [11] noticed in 1961 that small 
additions of certain amines to the hydrogen-ammonia system enhanced the 
rate of reaction. Bar-Eli and Klein [1] undertook to measure the rate 
of exchange of hydrogen-deuteride, tritium hydride and deuterium with 
aliphatic amines in the presence of lithium and potassium alkyl amide 
catalysts. Their study was conducted in a stirred reactor in a manner 
very similar to the study of the hydrogen-ammonia system [6]. The 


data for HD-CH NH, in presence of LiNHCH, indicates the rate constant 
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to be about 6.5 mtn: at -29°C and about 3.5 ine at -46°C. The 
maximum lithium amide concentration studied appears to be about 30 

mm/1. A comparison of the rates of exchange of deuterium in amino- 
methane and in ammonia indicated that the rate in the amine was up to 

35 times higher at -66°C. Based on the kinetic data obtained for the 
HD-CHNH, system an estimated value of the activation energy was reported 
to be 3.6 kcal/mole. 

The work of Rochard [2] also conducted in a stirred reactor inves- 
tigated the exchange of the hydrogen-aminomethane system between -62°C 
am -90°C and a potassium methyl amide catalyst concentration of 0 to 
70 mm/1. For a catalyst concentration of about 9 mm/1 a rate constant 
or’55 mint was observed at -62.2°C. In contrast to the value of the 
activation energy reported by Bar-Eli and Klein [1], Rochard [2] estimates 
the activation energy to be 6.8 kcal/mole and attributes the discrepancy 
to a different value of hydrogen solubility used by Bar-Eli and Klein [1]. 

2.2.1 Effect of Catalyst Concentration. The works of Bar-Eli and 
Klein [1] and Rochard [2] are the only ones dealing with the effect of 
catalyst concentration m the exchange of the hydrogen-aminomethane 
system. 

Bar-Eli and Klein [1] noted that at -22°C the rate of the reaction 
became invariant to increases in catalyst concentration beyond a value of 
10 mm/1. This effect was attributed to the saturation of the amino- 
methane with lithium methyl amide at 10 mm/Il. A linear variation of the 
rate constant Ky with catalyst concentration was observed ap to the 
suggested saturation value. Rochard [2] also observed a linear 
dependence of ky with catalyst concentration upto a value of about 20 


mm/1 at -90°C and suggested that the subsequent levelling off might be 
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due to the variation of the degree of dissociation with catalyst con- 
centration. By analogy with the work of Delmas [8] for hydrogen- 
ammonia, Rochard suggests that the free alkylamide ion [CHNH ] may 

be the only catalytically active specie. Additional work is, however, 
required to determine the degree of dissociation of the alkyl amide 
catalyst in the hope of establishing whether the alkylamide ion is 


in fact the only catalytically active specie. 


2.3 Catalyst Solubility 

Since catalyst concentration is an important variable for the 
kinetic study of the hydrogen-aminomethane exchange reaction it is 
desirable that the limits of solubility of the alkyl amide catalyst 
in the liquid be known. 

Hayashitani [12] conducted a study to determine the solubility of 
potassium alkyl amide in amines over a temperature range of -78°C to 
+40°C. His work indicates that the solubility of potassium methyl 
amide in aminomethane decreases with increasing temperature. Typical 
values of the solubility reported are 40 gms K/kg amine at -78°C and 


25 gms K/kg amine at 25°C. 


2.4. Reaction Mechanism 

The reaction mechanism for the hydrogen-aminomethane exchange is 
derived from the mechanism of the analogous hydrogen-ammonia reaction. 
Wilmarth and Dayton [5] suggested that hydride ions were formed in liquid 
ammonia solution as a reaction intermediate but this was contested by 


Bar-Eli and Klein [6] on the basis that lithium hydride failed to 
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catalyze the reaction. They believe that the free electron pair in the 
catalyst jon or ion pair initiates the exchange. This, they suggest, 
is done by stretching the D-H bond and partially forming a N-D bond. 
The solvent molecule then supplies a hydrogen atom and the exchange 


is completed by the transfer of charge as shown below. 


— = Ris : af : ed r ee : oo 3 
NT Csoly) HD = H N D:H Fada t N3 DH SH N H 
(solv) 
H H H 
Transition Intermediate 
State Comp 1 ex (2522) 


+ NH + H.ND 


a HOND + He He NE me =e Neoly) 9 


2 iz 


H 


Delmas et al. [8] suggest that the exchange occurs between a HD molecule 
and the solvated NH, ion. A N-H bond of the solvated ammonia molecule 
is polarized by NH, and then weakened. A four center mechanism is 


involved, i.e. 


‘4 = 
H - N|----- H-N-H Me te 
f / 
i x re Seu ae 
D-H : ; 
pi (2-23) 
H H(-) 
[ l 
> H-N+H+ |N -H 
le zai 
Dy 


Rochard [2] suggests that the above mechanism proposed by Delmas [8] is 


the most plausible one for the hydrogen-aminomethane exchange. 
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CHAPTER 3 


ANALYSIS OF MASS TRANSFER AND CHEMICAL REACTION OVER A SPHERE 


In order to obtain a theoretical support to the experimental study 


of mass transfer with chemical reaction in a liquid flowing over a 


sphere it is necessary to solve the equation which describes the above 


phenomena. In essence this equation is the equation of continuity of 


the solute specie and is written as 


wr tee VC AD Velho t (2-24) 
pre a The accumulation term 
v.VC = The convective transfer term 
pvc = The molecular diffusion term 
r = The rate of reaction (disappearance of the specie) 
D = The diffusion coefficient of the solute in the liquid phase. 
Equation (2-24) may be written in spherical coordinates to represent 


a solute mass 


to the set of 


balance in a liquid flowing over a sphere. With reference 


coordinates given in Figure | and neglecting convective 


terms in the R and $ directions and the diffusion terms in the 86 and 


d directions, 


equation (2-24) may be simplified to give 


+ Rr + <= (2-25) 
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Under steady state conditions and assuming a first order or pseudo- 


first order reaction, equation (2-25), becomes 


eel oC ae (of) 


= vy. (se ; + Rk, (C-C.) (2-26) 


The boundary conditions associated with equation (2-26) may be stated 


as follows: 


C = c. x =00 ey sage’ (= 2/3) 
eG. io _ 
Tk Woe x= A 8 >O0oO (2-27b) 

Ca Co x >0 8@=0 (2-27c) 


where A is the liquid film thickness. 

The boundary condition given by (2-27b) is used when the depth of pene- 
tration of the solute in the film is high so that the presence of the 
sphere influences the concentration profile. When the depth of penetra- 
tion is small (2-27b) may be replaced by 


C= Cy xX = © 8 > 0 (2-27d) 


which considers the film to be infinitely deep as compared to the depth 
of penetration. The low and high depth of penetration cases were 
reported for the case of pure physical absorption in a liquid flowing 
over a sphere by Davidson and Cullen [13]. 

When the reaction rate is such that the concentration in the bulk 
of the liquid is at chemical equilibrium then C, may be replaced by C. 
in (2-274). A fast reaction, in fact, has the effect of reducing the 
concentration level to C. very close to the interface thus steepening 


the concentration gradient in the film. This also implies that when a 
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chemical reaction occurs simultaneously with gas absorption the depth 
of penetration of the solute species is lower than the case for pure 
physical absorption. For higher rates of reaction this effect is more 
Pronounced. The boundary condition (2-27d) would thus apply in most 
cases where a fast reaction was involved. 

Ratcliff and Holdcroft [14] presented a solution to a simplified 


form of equation (2-26) wherein R was assumed much larger than A and 


2 
oC e e e 
ax Was neglected in comparison to R a Their solution was asymp- 
OX 


totic to the equation for physical absorption of Davidson and Cullen 
[13] but failed to approach the equation describing the gas absorption 


rate at high values of the rate constant namely 


G = 4nR? VD (C,-C.) (2-28) 


where (4nR7) is the interfacial area for absorption on a sphere. Astarita 
[15] pointed out this anomaly and assuming no stretching effect of the 
liquid film obtained the following solution to the original differential 


equation of Ratcliff and Holdcroft [14]. 


5 : exp(-k, t) 
G = 2R~ VOn oF | (aq erf (vk t) + ) sinads (2-29) 
0 t 
R 1/3 
wheres) t= =~ | Sin 6 dé (2-30) - 
V 0 


and ve is the interfacial velocity at the sphere equator. The boundary 
conditions used for equation (2-29) correspond to conditions (2-27a), 


(2-27c) and (2-27d) with Cy and C. assumed equal to zero. 
k,R 


For a fast reaction when kt or = Secret (vk t) and 
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exp (-kt) 


Vt 
equation (2-28). It is important to note that equation (2-28) is 


approach | and 0 respectively, and equation (2-29) reduces to 


independent of liquid flow rate which implies that at high reaction 
rates the hydrodynamics over the sphere do not affect the rate of 

gas absorption. This, however, is only true as long as the surface 
area of contact can be estimated by the expression ae For values 


of k,t or x <<] equation (2-29) degenerates to 


V 
fo) 


ee 2(C.-C.)R° “ov */R 


m  Sinéede 
1/3 v2 


) 
5 Uf Sin ada ) (2921) 


For k,t = 0, which corresponds to the case of pure physical absorption, 
Cc. the equilibrium concentration in equation (2-31) may be replaced by 
Co the bulk liquid concentration to give an equation which is identical 
to the one given by Lynn et al. [16]. In developing equation (2-31), 
however, Astarita [15] followed the approach of Lynn et al. [16] who 
neglected the stretching effect of the spherical liquid film. Wild 

and Potter [17] indicate that when the stretching effect is neglected 
the solution given by (2-31) for low reaction rates is about 18% in 
error. For high reaction rates the effect is indeed negligible. They 
solved equation (2-26) numerically for low depths of solute penetration 
and reported different equations representing the absorption rate in 


different ranges of k,- Defining 


2 5/3 
dy = ode = OE Sin NEE (2-32) 
vi 3LA 


where Vo? the interfacial velocity, and A, the film thickness over the 
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Shek: 


sphere, are given by equations (1-5) and (1-3) respectively, in the 
Firstepartkot this thesis. AT is the film thickness measured at the 
sphere equator and is given by equation (1-31). The maximum value 


of y for the sphere is given as 
v= | dy (2333) 
Le 


and Se the maximum contact time is expressed as 
1/3 
Sas ha le | sin ee (2-34) 


The solutions obtained by Wild and Potter [17] may thus be written as: 


(i) ee = 0 (physical absorption) 
te 
ates OOS (2-35) 
L(C,-C,) Wb, 
Cia) kt Pe hee 
G 
ees es, OS NR Tee + 1.693 (2286) 
1 “max 
L(Ci-C.) Y¥5 
Gl eOpy ai tee 530 
ear GA Paice Wet + 0.689 (2-37) 
L(C-C Mo BONE 
(iv) 5.0 < cree < 25.0 
ah 2 eee EE (2-38) 
L(C.-C.) Yo iin 
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and 

(v) kitty > 25-0 
Seen SO Poel Seria s IETS (2-39) 
L(C iC), 


[t may be noted that for equations (2-36) to (2-39) it is assumed that 
the reaction rate is high enough to keep the bulk liquid concentration 
at the chemical equilibrium value Co: 


Wild and Potter [17] also indicated that their numerical solution 





could be closely approximated over the entire range of salleass, by the 
equation 
T - 
sesepleh remem 5 0.714 ree | (ert 2 aie ) sineds (2-40) 
L(C.-C_) Wp, 0 vam 
where 
kt p : 
n= EES sin’ 36 | Sin?/? ede (2-41) 
: 0 


Equation (2-39) is representative of high values of Kitmas and is 
in fact identical to equation (2-28). Using equations (2-35) and (2-40) 


one may express the enhancement factor » due to chemical reaction as 


T 
g = 0.421 EE | (ertv oa ) sineae (2-42) 
max 0 Van 


Danckwerts [18] obtained a solution for $ for a semi-infinite slab 
model (i.e. neglecting the stretching effect of the film) which is 


expressed as 
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VT ] exp(-k)t..,) 
$= — ae O ae SA, + ———_—_—— | (2-43) 
f emen ae 0 


Wild and Potter [17] compared equations (2-42) and (2-43) for various 


values of 1S hired and indicated that for low reaction rates neglecting 
the stretching effect predicted absorption rates about 18% lower than 


the values obtained by considering the sphericity of the liquid film. 
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CHAPTER 4 


EXPERIMENTAL 


4.1 Apparatus 


4.1.1 General. The mechanics of the hydrogen-ami nomethane 
exchange experiment involve contacting the aminomethane flowing over 
the ball absorber with hydrogen gas which is enriched in deuterium. 
The hydrogen-deuteride (HD) transfers from the gas phase to the 
liquid where it dissolves and reacts. The reaction transfers the 
deuterium atom (D) to the aminomethane in exchange for a hydrogen 
atom (H) which combines with the remaining H atom of the HD resulting 


in a slight supersaturation of the aminomethane liquid with H This 


9° 
excess hydrogen gas then diffuses back to the gas phase. The consequence 
of the exchange is therefore to deplete the gas phase of its deuterium 
content. 

The solubility of HD in liquid aminomethane is very low and 
consequently the driving force for mass transfer in the liquid phase 
is small resulting in low gas absorption rates. This implies that if 
the gas phase were contacted with the liquid on a single pass basis 
the depletion of the deuterium in the gas phase would be very small 
which in turn would cause difficulties in detection of the small change. 
In order to alleviate this problem the gas phase was contacted with the 
liquid on a recirculating basis making the depletion of deuterium 
cumulative with time. A gas sample drawn from the circulating gas 
phase at a suitable time would indicate a deuterium content which is 


sufficiently lower than that of the previous sample to make the analysis 
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reliable. 

The potassium methyl amide-aminomethane (PMA) solution is hazar- 
dous to work with and certain general precautions should always be 
observed to ensure adequate safety. These precautions may be briefly 
summarized as follows: 

(i) Potassium methyl amide reacts violently with water libera- 
ting hydrogen which ignites creating explosive conditions in an air 
atmosphere. Moisture must therefore be kept out of the system at all 
times and contact with air must always be avoided. 

(ii) | When handling PMA solution in a refrigerated condition care 
should be taken to allow for volume expansion in case a warm up should 
occur or become necessary. Cold PMA solution should therefore, not be 
left trapped between valves. 

(iii) Aminomethane is toxic and a powerful irritant. Continued 
exposure to it can cause eye irritation and respiratory problems. 
Adequate ventilation and purge fans should be available when working with 
this system. 

(iv) The equipment used to feed and to store PMA solution must 
be protected against possible pressure buildup by providing pressure 
relief valves. 

(v) Amyl alcohol and methyl alcohol can be used to destroy the 
PMA solution and must be available at hand when required. Adequate 
fire fighting equipment should also be available. 

(vi) Face protection shields should be worn when handling PMA 
solution to avoid severe burns of the skin and eye damage from possible 


splashes. 
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The materials of construction used in the design of the hydrogen- 
aminomethane exchange equipment were limited to 316 stainless steel 
and glass due to the corrosive nature of the PMA solution. The stain- 
less steel equipment was designed for 500 psig pressure and temperatures 
lower than -200°C. Operating pressures and temperatures in the steel 
equipment during an experiment were in the range of 40-50 psig and -40°C. 
The glass equipment was used where visual observation was desired. It 
was designed for a pressure of 50 psig although the most severe operating 
conditions were limited to 15 psig and -30°C. Neoprene rubber was used 
for most gaskets and O-rings. Polyethylene gaskets were used only where 
liquid PMA solution contact was anticipated and polypropylene glands 
were used for all valve packings. Polypropylene tubing was used in 
place of 316 ss tubing where flexible lines were considered desirable. 

Lys lax: Hydrogen-Aminomethane Exchange Apparatus. The basic purpose 
of the hydrogen-aminomethane exchange system was to contact a circulating 
Stream of enriched hydrogen with PMA solution flowing over the sphere 
absorber on a once through basis. An adequate volume of liquid feed was 
available to allow the liquid to flow on a once through basis and yet 
provide a sufficient time for the experiment. The planned temperature 
range of investigation was between -10°C and -30°C and to achieve this 
the entire exchange apparatus was enclosed in a refrigerated air bath 
whose temperature could be controlled between -5°C and -30°C with a 
sensitivity of about +1°C. The refrigeration was provided by a ''York'! 
unit model F62C-502E, six cylinder, two stage internally compounded 
refrigerant 502 compressor, with a capacity of 7920 BTU/hr. at -100°F 


saturated suction temperature. Methyl hydrate was used as the external 
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heat transfer medium. Auxiliary refrigeration units of smaller capacity 
were used to service the feed and catalyst wash sections of the apparatus. 
Details of the equipment that constitutes the exchange system are shown 
in Figure 2. 

Pre-cooled liquid feed (PMA solution) from the liquid feed 
section (II!) came into the absorption system through a 'Whitey 21RS4'! 
316ss feed control valve (6) provided for accurate flow control. Any 
solid particles that might have been carried over by the liquid stream 
from the feed section were retained in a 60u size 3l6ss filter (15). 
The feed was then passed through a temperature stabilizing coil (16) 
and its flow rate was measured with a rotameter (18). A distributor 
(19) was provided in line to dampen any fluctuations in liquid flow 
before it was admitted on to the sphere absorber (23). The PMA solution 
passed over the sphere absorber through a jet nozzle (20) and the 
return liquid was collected in a take-off tube (24) from where it was 
led out of the air bath. Care was taken to keep the return liquid 
outside the air bath cold by installing concentric pipe heat exchangers 
(7,8) in the return line. This was done to ensure that no inverse 
pressure gradients, created by a possible fluid warm up, would impede 
the flow of the liquid to the receiving tanks (IV). A flow control 
valve (5) similar to the one in the feed line (6) was provided for 
accurate regulation. 

Details of construction of the distributor, the nozzle, the bal] 
absorber and the absorption chamber (25) have been included in the 


section on the diffusivity of hydrogen in aminomethane. 
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FIGURE 2 
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Gas-Liquid Exchange Apparatus 
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Symbols for Figure 2 
Amyl Alcohol Container 
Methyl Alcohol Container 
Alcohol Pump 
Drain 
"Whitey'' Control Valve on PMA Solution Return Line 
"Whitey'' Control Valve on PMA Solution Feed Line 
Concentric Tube Heat Exchangers on Liquid Return Line 
Saturator Amine Retrieving Tank 
Saturator Amine Feed Tank 
Pressure Gauge 
Gas Phase Circulation Pump 
Gas Saturator 
Solenoid Valve 
316 SS Micron Filter 
In Line Cooling Coil 
Thermocouples 
Liquid Feed Rotameter 
Liquid Distributor 
Jet Nozzle and its Approach Tube 
Gas Reservoir 
Pressure Balance Solenoid Valve 
Ball Absorber 
Take-Off Tube 
Absorption Chamber 


Gas Rotameter 
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Vi. 


1065, 


Refrigerated Air Bath Enclosure 
Gas Vent Line 
To Vacuum System 


To Enriched Hydrogen Feed Tank 


. To Liquid Feed System 


To Liquid Return Tanks and Catalyst Wash Tank 
To Gas Sampling System 


To Mercury Manometer for System Pressure Record. 
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Enriched hydrogen was introduced into the exchange system at the 
beginning of an experiment and then was kept circulating through the 
absorption chamber with the aid of a diaphragm pump (12). The 
diaphragm was constructed from neoprene rubber and had to be 
replaced periodically due to wear. In order to prevent evaporation 
of aminomethane from the surface of the sphere absorber the circulating 
ditcueh was kept saturated with aminomethane. This was accomplished 
by passing the gas through a column of pure, dry, catalyst free liquid 
kept in the saturator (13). It was introduced into the saturator from 
a separate feed tank (10) through a solenoid valve (14). If necessary 
this aminomethane could be retrieved through the same line into a 
chamber (9) immersed in liquid nitrogen. From the saturator the 
circulating hydrogen passed into a tank volume (21) where flow fluctua- 
tions were dampended and the gas was led into the absorption chamber 
through a rotameter (26). The volume of the gas phase is a parameter 
to which the time taken to deplete a given fraction of the deuterated 
specie is proportional. The tank volume was thus suitably selected so 
that sufficient time for exchange might be available. Gas samples for 
analysis were taken from this tank through the gas sampling section (V). 

In order to achieve a smooth gravity flow of the liquid over the 
sphere absorber it was necessary to balance the pressures in the liquid 
distributor and the absorption chamber. This was accomplished through 
a solenoid valve (22) connecting the vapor spaces of the two chambers. 
The pressure in the absorption chamber was measured through this line 
by a mercury manometer (VI). 


At the end of an experimental run the lines of the gas-liquid 
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exchange system had some residual PMA solution that could not be 
retrieved into the receiving tanks. In order to avoid plugging of 

the jines and the jet nozzle with solid catalyst due to evaporation 

of the amine from the PMA solution, this trapped liquid had to be 
disposed of in a suitable manner. This was achieved through the 
catalyst wash and destruction systems. The wash system as shown in 
Figure 3 was used to dilute the trapped PMA solution in the lines with 
pure aminomethane until the straw yellow color of the solution could 
no longer be seen. 

The remaining catalyst in the diluted PMA solution was destroyed 
using the catalyst destruction system as shown in Figure 2. In order 
to destroy the catalyst, methyl hydrate was pumped through the lines 
and the reaction products were led into the drain or waste disposal 
area. 

Potassium methyl amide reacts with alcohols to give a non-reactive 


alkoxide and aminomethane. The reaction is given by 


CHNHK + ROH > ROK + CH..NH, (2-44) 


where R is an alkyl radical. With an excess of the alcohol the products 
of reaction can be taken out of the system without difficulty. The 
reaction (2-44) generates substantial amounts of heat the magnitude of 
which depends on the alcohol being used. Higher alcohols give a mild 
reaction but are disadvantageous from the standpoint of subsequent 
system evacuation because of their low vapour pressures. 

During an experiment, hydrogen gas samples were taken at suitable 


intervals from the tank volume of the exchange apparatus into a line 





801 


od jon blues sedt noituloe AMG Isubize1 smor ben megeye spnersks 


to enippulg biove o} ashe mn! 6 .ained enivisosa eng osm! bavalajzat 


noitsvoqevs o7 oul Jeyieisa biloe di2iw elsson Jol sd2 bas een)i ead 


ad oJ ben biupil beqge s zidz ,noituioe AMY ond mor? enims afd to 
aly deuoidd baveidos caw 2inT .180nem sidetive 6m) to beeogeib 
ni nwode 26 meseye Heew ofT .emeteye noldoursesb bas deew Jeyiedso 
fisiw eenif sft ni moljuloe AMI beqgqe 1s off osgulib.od beeu 2eow € saupit 
biluos aoiguloe sft to toloo wolley wevte end litnu enendesmonims sug 
-naer ed tepnal on 

bsyovs2sb caw noituloe AMS besullb edd ni seyistes paiaiemes sat 
tabio nl .S axupit ni awore 26 med2zye nolzouiteseb tdeyieges ed} onie 
eani! ads Hpvoils beqmuq esw syerbyd tyrison .Jeylaiso eng yous zab 93 
is2oqzib eteew 10 nierb eft ofn] bel ovsw etoubesg noldasen sit bas 
, 6976 

svijoset-non 5s svip o3 2tonools djiw esoee, ebims fydjam muleeeso4 


yd covip 2i nolijes) scT .snadsemoqime bas ebixolis 


oD + OR» HOR + nin AD 


~ 


(AA-s) git 


@iavberg efi forosls sit io 2asaxe née d3iW -leoiber tylls oe 2i A stedw 
eT .viiu2itlib tuorsiw meseyz st3 to tuo Aedpd 9d ns2,n0ijose1 to 


Hoh soar *o adteanem: Soisntcadie aaemaaaig iit) malas 





a Me ae 


109. 


wWazSAS ysemM 2SAL ee) BUL 


S sheliisdke 


WSS SSA SAN MSAANSSAS N 





OmOmMOme YY) 
OZ OTE ea a I A 


= 


v 
Ne 
No 
N 
El 
RE 
Re 
Np 


00! 









10 


ee, eS, ee ee ee ee 





MLD CR ee AM aT 


109a. 


symbols for Figure 3 


1, I-A Aminomethane Holding Tanks 

2, 2-A Mild Steel Drums to hold liquid nitrogen when required 
3, 3-A Steel Straps 

4. Cooling Bath with Ethylene Glycol-Water Mixture 
Dee atn@ Line Cool ing? Coll 

6. Intermediate Holding Tank 

7. Catalyst Wash Receiving Tank 

8. Glas Wool Filters 

9. Mercury Manometer 

10. Pressure Gauge 

fs To Vacuum System 

[1. To Liquid Feed-Receiving Tank System 

ltl. From Gas-Liquid Exchange Apparatus 

IV. To Catalyst Destruction System (Alcohol Lines) 
V. Aminomethane Wash Feed to Exchange Apparatus 


VI. Alcohol to Drain during Catalyst Destruction Cycle. 
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whose volume was about 5cc. The sample was then expanded into sample 
bombs. The aminomethane vapours in the gas sample were condensed in 
chilled traps. At the end of a run the samples were analyzed using 
mass spectrometry. 

In order to characterize the data completely it was necessary to 
take PMA solution samples and analyze them for concentration of 
potassium. This sampling was not done as a regular feature at the 
end of every experiment but was done frequently to ensure that the 
catalyst concentration was well established. The analysis was done on 
an atomic absorption spectrophotometer. 

Additional equipment used to carry out the various operations 
outlined above, for a hydrogen-aminomethane exchange experiment, is 
described below. 

4.1.3 Liquid Feed System. The liquid feed system as shown in 
Figure 4 was made up of two parallel trains which served as the feeding 
tanks (1-1A) and receiving tanks (2-2A) on an alternating basis. Each 
set of tanks was made from two 6 inch diameter, 316 SS chambers immersed 
in a bath of ethylene glycol-water mixture (11,12). The temperature 
of the bath around the feeding tanks was controlled at 25°C and the 
bath around the receiving tanks was kept refrigerated at -40°C. Methyl 
hydrate was used as the heat exchange medium. A constant vapour pressure 
of the PMA solution at 25°C was thus used to feed the liquid to the gas 
absorption chamber. Liquid was fed from the feeding tanks through two 
dip tubes and was cooled en route to the absorption chamber by concentric 
pipe heat exchangers (13, 14, 15 16) and an intermediate cooling bath 
(10) also containing ethylene glycol-water. Adequate residence time 


was provided in the cooling bath by introducing an in-line cooling coil 
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FIGURE 4 


The Liquid Feed System 
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Symbols for Figure 4 


1-1A. Liquid (PMA) Feed/Return Tanks (Ist Set) 
2-2A. Liquid (PMA) Feed/Return Tanks (2nd Set) 
3. Pressure Relief Valve 

4. Pressure Gauge 

5. In Line Cooling Coil 

6. Intermediate Holding Tank 

7. Glass Wool Filter 

8. Cooling Coil carrying Refrigerant 

Few tivanes 

10-12. Cooling Baths 

13-16. Concentric Tube Heat Exchangers 

Le To PMA Solution Storage Tank 

11. PMA Solution Feed Line from Storage Tank 
l11. To Argon Pressure Cylinder 

IV. To Alcohol Lines for Catalyst Destruction 
We To Gas-Liquid Exchange System 


Vi. To Vacuum System 
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(5) and an intermediate tank (6). Secondary filtering of the PMA 
solution was achieved through a glass wool filter (7). Each set of 
tanks was provided with vacuum, argon blanketing, pressure measuring 
and pressure relief services. 

4.1.4 Catalyst Wash System. Two stainless steel holding tanks 
(1,1A) (Fig.3) containing pure aminomethane were supported on the edge 
of two outer mild steel drums (2,2A) by steel straps (3,3A). The 
holding tanks were connected in parallel to ensure an adequate supply 
of aminomethane and to balance pressures. The aminomethane was pre- 
cooled by passing it through a cooling coil (5) and a tank volume (6) 
both immersed in a refrigerated bath (4) of ethylene glycol-water. 

The wash amine was then filtered through a glass wool filter (8) before 
being finally led to the gas-liquid exchange section. The diluted PMA 
solution was returned into a mild steel refrigerated holding tank (7). 
After three or four experiments the bulk of the amine in tank (7) was 
flash distilled into the original holding tanks (1,1A) for re-use. 

The flash distillation was achieved by filling the mild steel drums 
(2-2A) surrounding the holding tanks with liquid nitrogen and connecting 
the vapour spaces of the holding tanks (1,1A) with tank (7) which was 
brought to room temperature before initiating the distillation proce- 
dure. 

4.1.5 Catalyst Destruction System. Two glass bottles (1,2) (Fig.2) 
containing amy] alcohol and methyl alcohol respectively, were connected 
to the gas-liquid exchange system through a centrifugal pump (3). The 
amyl alcohol was kept stocked as a standby in case of an emergency 
situation. The methyl alcohol was used to destroy the remaining 


catalyst in diluted PMA solution left in the lines of the gas-liquid 
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exchange system after the wash operation was completed. The alcohol 
was pumped into the absorption chamber through the PMA feed line, 
return line and the drain line. These were the only lines that come 
in contact with the liquid. The same pump was then used to pump 

the liquid out of the chamber into a drain (4). The above procedure 
was followed only when some PMA solution spillage had occurred during 
an experiment. If such was not the case then the alcohol was simply 
pumped through the liquid feed line and out through the return line 
into the drain. 

4.1.6 Catalyst Make-up Section. Potassium methyl amide catalyst 
was made by contacting freshly cut potassium metal with dry aminomethane 
in a 10 inch diameter 316 SS storage tank (1) (Fig. 5) of about 35 
litres total capacity. 

The potassium metal was cut into small pieces, washed in n-pentane 
and introduced into the tank through a 1/2 inch bore under a blanket 
stream of argon. This was done to prevent contact of air and moisture 
with the highly reactive metal. The storage tank was then evacuated 
through the main vacuum system (IV) to eliminate traces of air, argon 
and n-pentane. When a vacuum of the order of 1-5u was obtained, the 
required amount of dry aminomethane was flash distilled into the 
storage tank (V) by codling the latter in a bath of liquid nitrogen. 

The frozen aminomethane was then allowed to thaw out and on doing so the 
potassium metal began dissolving in the liquid. This could be seen by 
observing the blue color of the solution through a sight glass (9). 

The solution of potassium in aminomethane is a slow step in the preparation 
of the PMA solution. In order to enhance the rate of solution the 


aminomethane was circulated in and out of the storage tank by a 316 SS 
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The Catalyst Make-up and Liquid Sampling Section 
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symbols for Figure 5 


PMA solution Bulk Storage Tank 

Port for Introducing Potassium Metal 
Alternate Line for Removing Feed from Tank 
Pressure Gauge 

Pressure Relief Valve 

Glass Wool Filter 

Sampling Bomb Assembly 

Polypropylene Diaphragm Pump 


Sight Glass 


- Outer Vessel for holding Liquid Nitrogen When required 


Return from Feed System 


Liquid Feed to Feed System 


- To Argon Pressure System 


To Vacuum System 


Feed Line for Flash Distilled Dry Aminomethane Vapour 
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positive displacement pump with a polypropylene diaphragm (8). Any 
solids foreign to the system were retained in a glass wool filter (6) 
installed in the circulation loop. A by-pass was installed so that 
liquid samples could be drawn from the circulating stream by introducing 
a pre-assembled sample bomb (7) in line with the system. Other by- 
passes were installed to isolate any piece of equipment from the PMA 
circulation line, should any Pea Tee anen Cir maintenance become necessary. 

Potassium metal reacted with aminomethane forming the amide and 
liberating hydrogen. The liberated hydrogen resulted in an increase 
in the total pressure of the storage vessel. This could be monitored 
on a pressure gauge (4). A pressure relief valve (5) was provided to 
bleed off any pressure in excess of a safe limiting value. A dip tube 
(3) was also installed so that the entire charge could be dumped if 
necessary. This was done as an extra safety measure. 

4.1.7 Gas Sampling System. Gas samples were drawn from the gas- 
liquid exchange system during a run, into an evacuated line of about 
5cc volume. This sample was then expanded through a set of liquid 
nitrogen chilled traps (1) (Fig. 6) into an evacuated sample bomb. 
Aminomethane vapours carried over with the gas sample were condensed 
in the chilled traps. A mercury manometer (3) was used to measure 
the pressure in the sample bomb. After the sample was obtained, the 
lines were completely evacuated in readiness for the next sample. 

The pressure in the absorption chamber was measured by a mercury 
manometer (9). In order to know the pressure at which the sample was 
drawn from the exchange system it was necessary to know the drop of 


pressure in the absorption chamber after each sample was drawn. This 
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Symbols for Figure 6 


1. Aminomethane Vapour Traps 

yop Gas Sample Bomb 

oe Sampling Pressure Manometer 

4, Ultra Pure Hydrogen Cylinder 

sic Pressure Regulator Valve 

6. 15 Litre Tank Capacity filled with Ultra Pure Hydrogen 
73 Mercury Manometer to measure pressure in (6) 

8. Differential Pressure 0i] Manometer 

om Mercury Manometer on the Gas Absorption Chamber 

10. Of] Vapour Trap 

1]. To Vent 

12. Valve isolating Sample Volume from Vacuum Pump. 

13. Valve isolating Sample Volume from the Gas-Liquid Exchange Apparatus 


14, Glass Valve isolating Sample Volume from the rest of the Sampling 
Apparatus 


15. Glass Valve to Sample Pressure Manometer 

16. Glass Valve to Sample Bomb 

17,18 Glass Valves on either leg of Sample Pressure Manometer 
19. Valve to introduce ultra pure Hydrogen in (6). 

20. Valve for Excess Pressure Bleed 

21. Valve isolating (6) from Vacuum Pump 

ie To Gas Absorption Chamber 

[1. From Gas Reservoir in the Gas-Liquid Exchange System 


l11. To Vacuum Pump 
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drop was small enough to go undetected on the mercury manometer (9) and 
was thus measured with an oil manometer (8). Octoil-S, a high vacuum 
diffusion pump oil was used as a manometric fluid. This oil has a 
very low vapour pressure which is desirable to prevent contamination 
of the exchange system. A small ice chilled trap (10) was also used 
between manometers (8) and (9) to reduce the risk of contamination. 

One leg of the oil manometer was hooked to the absorption chamber 
along with the mercury manometer (9). The other limb was attached to 
a large fifteen litre tank (6) which contained pure hydrogen. The 
pressure in this tank was measured by another mercury manometer (7). 
At the beginning of an experiment the pressure in the tank (6) was 
set to the starting pressure of the absorption system by matching 
manometers (7) and (9). This was accomplished by adjusting the pressure 
in (6) so that both limbs of the oil manometer (8) read zero. A hydrogen 
supply cylinder (4) was used to introduce gas into tank (6) if required. 

4.1.8 Analysis. 

CajmGas Samples 
The gas samples taken from the exchange system during 

an experiment were analyzed for deuterium content by a mass spectrometer. 

The instrument, designed and built by Buckley [19] is a 90° sector 
analyzer with a radius of 2.5 inches and employs dual collection of 
the isotopes D and H. The sample is introduced into the instrument 
as a gas and the signal from the collection cup for each isotope is 
amplified and read on vibrating reed electrometers. The electrical 
signal of the electrometers is also sampled on-line by an IBM 1800 
process computer which processes the data and computes the results 


in accordance with pre-programmed instructions. 
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(b) Liquid Samples 
Liquid samples of the PMA solution were taken 

occasionally from the feed tanks for analysis of total potassium. 
Atomic absorption spectrophotometry was used to estimate the concen- 
tration of the catalyst. The sample was weighed, hydrolyzed and 
diluted to known volumes. A set of potassium standards ranging from 
0.0 to 5.0 mg K/litre solution were prepared and analyzed on the Atomic 
Absorption Spectrophotometer (AAS) to generate the calibration curve. 
The sample was then analyzed and its potassium content was estimated 
from the calibration curve. The catalyst concentration was computed 
from the potassium content of the sample and the known dilution 
procedure. The AAS was a model 290-B ''Perkin-Elmer'' unit which employs 
a hollow cathode type ''Intensitron Lamp'' #566-M. The cathode of the 
lamp contains the element which has to be analyzed. The sample is 
aspirated and burnt through a flame which breaks the chemical bonds 
of the compound and the individual atoms of the element of interest 
float freely in the absorption area of the instrument. Light emitted 
from the element of hollow cathode lamp is passed through the flame 
and depending on the concentration of the sample a certain proportion 
of the light is absorbed. The instrument measures the amount of absorp- 
tion. 

4.1.9 Hydrogen Make-Up and Purification. A Milton-Roy hydrogen 
generator was used to produce a supply of hydrogen enriched with deuterium. 
The feed to the generator was enriched water which was electrolyzed 
to give hydrogen. The enriched gas was passed through a series of 
three beds to remove impurities before entering the gas absorption 


chamber. The beds consisted of activated charcoal for hydrocarbon 
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4.2 Operating Procedures 
4.2.1. The Exchange Experiment. The operating procedure used for 


conducting the hydrogen-aminomethane exchange experiment was similar 
to the procedure for the measurement of diffusivity of hydrogen in 
aminomethane as far as control of liquid flow rate, take-off level and 
system temperature of the gas-liquid exchange section was concerned. 
However, some basic differences between the exchange and diffusivity 
experiments introduced important changes in the running of the equip- 
ment. A step by step operating procedure for a typical experimental 
run is outlined in Appendix A2. 

One inch and one and one half inch diameter sphere absorbers were 
used for the exchange experiments. The large ball was found to be 
useful in effectively reducing the time taken for a given fractional 
depletion and thus improved the chances of reaching chemical equilibrium 
during the period of the experiment. A jet diameter of 0.046 cms and jet 
lengths of 0.2 to 0.265 cms were used to obtain the entire data set. 

In concurrence with the findings of the diffusivity work, a take-off 
length of 1.25 cms was always maintained and it was assumed that for 
this value of the take-off length the stagnant layer mentioned in the 
first part of this thesis was always at the junction of the sphere and 
the take-off rod. The take-off rod was 1/16 inch in diameter. 

Pre-cooled liquid feed was introduced into the refrigerated, 


controlled temperature exchange apparatus to contact with enriched 
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hydrogen which had been previously saturated with aminomethane and 
which was kept circulating at a rate of about 100 cc/sec. The total 
gas volume of the system was about 2500 cc which meant that the turn 
around time for the gas was about 25 seconds. A liquid flow rate of 
about 0.3 cc/sec was maintained by controlling the feed rate with the 
‘Whitey'' control valve. This flow rate was adequate to ensure that 
there was minimal enrichment of a liquid element during its residence 
over the sphere absorber. At the above mentioned liquid flow rate 
there was no visual evidence of any ripples being present on the 
sphere absorber. Minor disturbances of the liquid were, however, 
observed on the take-off tube level and were perhaps indicative of 
some activity on the take-off rod. Since the diameter of the sphere 
absorber was 24 times larger than the rod some activity on the latter 
could be expected. This only emphasized the importance of keeping the 
take-off level well controlled at 1.25 cms. 

A gas sample was taken before the liquid PMA solution was introduced 
into the system and subsequent gas samples were generally taken at 10 
minute intervals during the exchange experiment. This was achieved 
by expanding a 5 cc sample into an evacuated glass sample bomb. Amino- 
methane vapours in the sample were condensed in a trap immersed in 
liquid nitrogen. Before allowing a sample into the sample bomb, the 
equivalent of one sample volume (~5cc) was discarded. This was done 
to eliminate the effect of dead volume in the line up to the tank from 
which the sample was drawn. 

In order to justify the assumption of an essentially constant 


composition of the liquid with respect to the deuterated specie the 
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minimum molar ratio of liquid to gas in the system during operation was 
about 300. This could be achieved by running the experiment at a 
hydrogen partial pressure of about 1.2 atmospheres with a maximum total 
system pressure of about 2.0 atmospheres at -10°C. 

After the desired number of samples were obtained, the feed was 
stopped and the PMA solution in the system was drained into the 
receiving tanks. Any left over solution in the lines was diluted with 
pure aminomethane and the wash was separately stored. Finally the 
system was washed with methyl hydrate to destroy any remarents of catalyst 
in the system. The outgoing methyl hydrate, aminomethane and potassium 
methoxide mixture was drained. 

4.2.2 PMA Solution Make-Up. The handling of potassium or lithium 
metal requires adherance to safe working procedures so that a hazardous 
situation might be avoided. The most prominent one of these is the 
necessity of keeping the metals out of contact with air or moisture 
at all times. The work of Hayashitani [12] has outlined details of 
the methods recommended when working with potassium metal. These 
methods were followed quite rigidly. 

A high catalyst concentration in the PMA solution results in 
high rates of exchange which is a desirable feature from an industrial 
standpoint. Catalyst concentrations upto 15.8 gms K/kg Amine were 
prepared for this study. 

Aminomethane was obtained in cylinders from Matheson of Canada 
Ltd. ts analysis report indicated a maximum of 0.8% moisture content 
which was removed by contacting the required aminomethane charge 
(about 28 litres) with a slight excess of lithium metal in the absence 


of air. The reaction of lithium with water produced lithium hydroxide 
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and hydrogen. The equations that describe the reactions of lithium 


with aminomethane containing moisture are: 


Li +H,O > LiOH + 3H, (2-45) 


Li (excess) + CH,NH, > CHNH Li + 3H, (2-46) 
The hydrogen produced from these reactions increases the pressure of the 
system which in turn inhibits the reaction. The gas was bled off 
periodically and the pressure kept at about 30 psig at room temperature. 
The reaction was assumed to have gone to completion when no further 
increase of pressure was observed. The entire process of removing 
water from the aminomethane took about 48 hours. 

The PMA catalyst solution was prepared in much the same manner 
as the drying of aminomethane. The amount of potassium metal required 
for a charge of about 28 litres of PMA solution was cleaned, cut and 
contacted with dry aminomethane in a 316 SS storage vessel in absence 
of air or moisture. The aminomethane was flash distilled into the 
vessel by immersing the latter in a bath of liquid nitrogen. On 
thawing out, the aminomethane reacted with the potassium according 


to the following equation. 


1 SS 
K + CHNH, > CH NHK + 4 H, (2-47) 


This reaction is catalyzed by traces of iron which were supplied by 
the walls of the stainless steel vessel. 

Potassium metal dissolves in aminomethane giving a blue color to 
the solution. It then reacts to form the amide and on completion of 


the reaction the color changes to straw yellow. The solution step is 
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however, the slower of the two and as such the liquid of the make-up 
vessel was kept circulating so as to improve the rate of dissolution. 
As in the case of the lithium drying of aminomethane, the hydrogen 
pressure buildup was continually relieved to keep the reaction pro- 
gressing. In spite of these measures, however, the reaction may still 
take a week or more to go to completion. A detailed sequential pro- 
cedure for the make-up of PMA solution is included in Appendix A2. 

4.2.3 Analysis 

(a) Gas Samples 
The mass spectrometer used for the analysis of gas 
samples was calibrated using hydrogen standards containing 1.5, 99, 
365 and 471 ppm deuterium. 

The 1.5 ppm standard was prepared by collecting the hydrogen gas 
evolved by contacting a water standard of 1.5 ppm deuterium concentration 
with calcium metal. The water standard was supplied by the Chalk River 
Nuclear Laboratory (CRNL) of AECL. The other three standards were 
prepared by diluting a high deuterium concentration volume of hydrogen 
obtained from the hydrogen generator with ultra pure hydrogen of natural 
concentrations obtained from Matheson of Canada Ltd. These three 
standards were analyzed at CRNL relative to the 1.5 ppm water standard 
and the resulting values reported above were taken to be the absolute 
deuterium concentrations of the samples. 

The mass spectrometer was calibrated before and after the analysis 
of a set of samples belonging to an exchange experiment. The current 
ratios (1,/1,) obtained from the mass spectrometer for the standards 
were compared with the absolute D/H ratios through the following 


relationships. 
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Rx 
soi" (F be ) 1000 (2-48) 


of 
Ef osa) 
6 (2 - } 1000 (2-49) 
st 


where R* = The ratio of mass three (HD) to mass two (H.) currents 
(1,/1,) for the sample at a fixed value of i, (also 
described as a fixed value of hydrogen pressure). The 
ratio In/l, is uncorrected for H~ contribution. 


3 


oe = The ratio I,/l, for the principal standard at the same 


value of l,- 
C, = Absolute value of the sample in ppm deuterium. 


> = Absolute value of the principal standard in ppm deuterium. 


The principal standard chosen for this study was the 365 ppm deuterium 
sample but in principle any of the 99, 365 or 471 ppm standards could 
be chosen as the principal standard. I and I, are respectively, the 
current signals generated by the HD and Ho beams bombarding the collector 
caps of the mass spectrometer. A plot of | and coe was considered to 
be the calibration curve of the mass spectrometer. The plot was 
essentially a straight line whose equation was determined by standard 
least square techniques. A typical example of this curve with the 
plotted points representing the standards is shown in Figure 7. The 
values 'R' obtained from the mass spectrometer for each of the samples 
taken during a run, were converted to the corresponding 1S using 


equation (2-48). The least square equation was then used to determine 


the 64s from which the absolute values yee were obtained by equation 
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(b) Liquid Samples 


Liquid samples were taken from the system periodically 
and analyzed for potassium methyl amide using the atomic absorption 
spectrophotometer (AAS). 

Two 316 SS sample bombs consisting of a sample volume and an 
expansion volume were thoroughly evacuated and then used to sample the 
PMA in line. An average sample was in the range of 4 to 6 gms and was 
taken in duplicate. The sampled PMA was then hydrolzed in distilled 
water and the resulting solution was heated in an oven at 200°F to 
drive off the aminomethane. On most occasions the solution was 
evaporated to dryness but this is not essential. Usually heating for 
about four hours is sufficient to drive off the aminomethane. 

Six liquid standards containing 0.0 (distilled water), 1.0, 2.0, 
3.0, 4.0 and 5.0 mgms K/litre solution were prepared by dissolving 
appropriate amounts of KCI in distilled water. The dried sample was 
re-dissolved in distilled water to make up a volume such that its 
estimated concentration of potassium was between 0.0 and 
5.0 mgms K/litre solution. The standards and the samples were then 
analyzed on the AAS. A calibration curve of potassium concentration 
vs AAS reading was prepared by using the standards. The AAS reading 
obtained from the samples was then converted to potassium concentration 
by using the calibration curve. The concentrations thus obtained were 
interpreted to give PMA concentrations in gms K/kg amine by using the 
dilution procedure and the weight of the samples. A typical calibration 


curve for the AAS is shown in Figure 8. 
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4.2.4 Method of Interpretation of Data. The number of moles of 


HD in the circulating gas phase deplete cumulatively with time. This 


depletion rate may be expressed as 


d ae ae 
-n oh i ae BALD Ph 2 = 
ake Rance wae a) 
g 
where a = Total number of moles of hydrogen gas. 
2 


P = The partial of hydrogen (atm) 


V = Total volume of the system occupied by the gas (cc) 
R = Universal gas constant = 82 .06(cc) (atm) 
g gm mole K 


T = Absolute temperature of the circulating hydrogen (°K) 
y = Mole fraction of hydrogen deuteride (HD) in hydrogen (H.) 


t = time (sec) 


The negative sign indicates a decreasing concentration y of HD with time 
in the gas phase. 
A mass balance of the hydrogen deuteride across the gas-liquid 


interface yields 


: Rt e =k AS (Cc - Cc.) (2-51) 
where A. = The surface area of contact for the dry sphere (cm?) 
KP = The overall mass transfer coefficient (cm/sec) 
Gs = The concentration of HD in the liquid phase in physical 


equilibrium with the HD concentration in the bulk gas phase 


gm moles 
Ce 
Ce = The concentration of HD in the liquid in chemical equilibrium 
with the concentration of CH,NHD in the bulk liquid 97 moles 
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Since the molar ratio of the liquid to gas used was always very high 
the value of CHNHD in the bulk liquid was essentially constant which 
implies that C. is a constant for a given experiment. Assuming that 


the total resistance to mass transfer can be represented by a two 


film theory concept one can say that 





ge 
re ieee os aoe 
K ke Hees 
where oe molar density of the liquid an moles 
7 ais atm 
Ho The Henry's Law coefficient eT SG 
k, = Liquid phase mass transfer coefficient sal 
L sec 
%e = Gas phase mass transfer coefficient SESS 


cm sec atm 
Since the solubility of hydrogen in aminomethane is very low (i.e. H, 
is very high) the resistance to mass transfer is essentially in the 


liquid phase and therefore 


~ (2-52a) 


4 
and C€ corresponds to the interfacial concentration as well and may be 
replaced by C.. Representing the solubility of HD in aminomethane by 


the Henry's Law as 


yP =P, =H, x, (2-53) 
pie toa ae (2-54) 
oT 
2 
where P. = Partial pressure of hydrogen deuteride 
x, = Mole fraction of HD in the liquid phase. 
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Noting that 
Se ea 
(2-55) 
and Fo ool ttc 


Equation (2-51) may be expressed as 
(2=56) 


Assuming that the concentration of HD in the gas at time t=o is given 


by Y.,2 and replacing y, by y, integration of equation (2-56) between 


0 and t leads to 
Riel 
(2-57) 
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lf the interfacial area of contact for the sphere is corrected for the 


thickness of the film then A, in equation (2-57) may be replaced by 


An where 
A + 
“AL 2.58 A : 
Aaah LE aoa R (2-58) 


Be being the surface area of the dry sphere and AY the liquid film 


thickness at the sphere equator. 
Yale 
A plot of In ace vs t yields a slope of 


in ‘e 
K An Py Rol 
Ya from which KP the mass transfer coefficient may 
2 


For high values of the chemical rate constant ky (i.e. 


be obtained. 


in the fast reaction regime) equation (2-6) may then be used to determine k 


l 








.0¢1 “Tess ic ¥ ; : mt 
sail hia 
7* jo" po as | > 


(22-¢) | ma, 
* 5” ‘bas 


26 bezeorgée od yam (12-S) noldsup3. 
eae | _ xb a 
(a2~s) (Vy - sy) gga A Se gee 


navip zi o=3 smid te ese sit ni GH to nolsexdnsonad efd PeAd enimuezA 








neewied (d92@-S$) noitguns To nolgexnpssni .¥ Vd ¥, oniséiqes bas watt yd ; 
oi ebse!l 3 bas 0 ~ 
T § 46 A Y"V ; 
a 1° 2 4 3 
({a=2) 5 —ts aule ( = : Jar ;, - 
% a s' ni ; 
si 207 dst0297705 2) sisiqe say 10% so6%nos To sete’ Islostiezai snd Th 
yd bes6igs1 ad vam (Ye-8) doigeups ni oA morls mf i? aid Yo 2eendsing 
> wl 
onerw 4A 
> 
+ ; g 
(82-5) ase hs x 
2. 
ve 
mbit biupil sz *) bas eisige yb siz to sone esstiue si nied . A 


= 


-Tateups srerge ef Js <a ' 


i 
Io sqole 6 eblsly tf 2¥ * (sia) Yo ee ' 





1Ba% 


Urey K, = vk.D (2-6) 


For lower values of k, equations (2-35) to (2-38) may be suitably 


rearranged to isolate K, as a function of k 


L and known system properties. 
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CHAPTER 5 


RESULTS 

Experiments were conducted to measure the rate of exchange of 
hydrogen with PMA solution at four temperatures in the range of -10°C 
to -30°C. Both the 1 inch and the 1.5 inch diameter sphere absorbers 
were used to obtain the data. Three catalyst concentrations were 
investigated in the range of 5.0 to 16.0 gms K/kg amine. The rate 
constants and the mass transfer coefficients calculated from the data 
are given in Table 2. Equation (2-6) was used to relate K, with k, 
and the surface area exposed by the sphere absorber and the jet was 


considered to be the interfacial area of contact. A detailed sample 


calculation showing the procedure used is given in Appendix B2. 
Yaya 


Vine 


Fi e e f 

REIS SP UES SULA §eheKe te In ( )vs time and indicates 
the pseudo-first order nature of the exchange reaction. The ratio 
YoYe 


vinge 


is indicative of the cumulative depletion of hydrogen deuteride 


( 





in the gas phase. The equilibrium value y, was estimated either by 
direct gas analysis or from the liquid concentration of deuterium and 
the separation factor defined by equation (2-17). Equation (2-21) was 
used to represent the separation factor as a function of temperature. 
For cases where the experiment was allowed to reach equilibrium, Yo. was 
obtained directly by gas analysis and the corresponding liquid phase 
equilibrium concentration was estimated using equation (2-17). Once 
known, the liquid phase concentration was suitably adjusted to account 
for liquid enrichment at the end of each experiment. This was achieved 


by a simple mass balance of deuterium between the gas and the liquid 
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phases. Details of this calculation are given in Appendix B2. Owing 
to the very large molar ratio of liquid to gas (~300) the enrichment 

of the liquid for each experiment was seldom above | ppm which meant 

a change in the gas phase equilibrium value of about 0.2 ppm. 


The logarithm of the rate constant k, is plotted vs the reciprocal 


I 
absolute temperature in Figure 10. The relation of Ky with temperature 
may be expressed by (a) the Arrhenius Equation, or (b) the Transition 


State or Absolute Reaction Rate Theory. The temperature dependency 


for the two cases jis expressed as 





a me 
RG 
a) Kine ikee (2-59) 
and J 
' R_T 
b) k, =k Te? (2-60) 
fe) 
where s 
Ko = pre-exponential constant (sec ) 
i at ra 
eae pre-exponential constant (sec ok 4 
= 4 cal 
ES = Activation energy ie 
T = absolute temperature (°K) 
” cal a. 
RG =| gaseconstant ——-=s7 = 1.987 


Taking natural logs of equations (2-59) and (2-60) one obtains 


E 
a) In k, = Ink, - er (9-61) 
g 
| ES 
and b) In k, = In k, + InT - aT 
g 
; kK, ; ES 
eye ha k, = In <a In ko -o (2-62) 
g 


A linear model fit to equations (2-61) and (2-62) suggests the following 


parameters based on a least square estimation 
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ae ln ko an bc 
Bo = 3.554 kcals/gm.mole (2-63) 
and b) Ink. = 6.641 
26 = 3.053 kcals/gm.mole (2-64) 


The activation energy of 3.6 k cals/mole reported by Bar-Eli and Klein 
[1] was apparently determined using the model represented by equation 
(2.59). Rochard [2] on the other hand used equation (2-60) to calculate 
the activation energy of 6.8 kcals/mole reported in her work. She 
attributed the discrepancy between the values obtained by her and by 
Bar-Eli and Klein [1] to different values used for the solubility of 
hydrogen in aminomethane. This is in contradiction to the fact that 


the Henry's Law coefficient at -62.2°C calculated from equation (2-13) 


and the value reported by Rochard [2] are 1.24 x 10" eee ee rend 
mole fraction 
Waghele 5:2 107 ny respectively and indicate a reasonable 


mole fraction 


agreement. 

The variation of the rate constant Ky vs the catalyst concentration 
in the PMA solution is shown in Figure 1]. The catalyst concentration 
was analyzed in duplicate and a mean of two such duplicates was taken 
to be the correct concentration. The maximum deviation of the samples 
about the mean was within +1%. Figure 11 indicates that the rate con- 
stant apparently increases almost linearly with catalyst concentration 


"c' less than about 5.0 gms K/kg amine and subsequently 


for values of 
levels off and appears to become independent of 'c' at values larger 


than about 20.0 gmsK/kg amine. 
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CHAPTER 6 


DISCUSSION 
6.1 Gas Phase Resistance 


[It may be seen from Table 2 that the values of k ~ for each 


i'ma 

experimental determination are larger than 25.0 which, in light of the 

theoretical development of chapter 3 indicates that the exchange 

reaction always took place in the fast reaction regime for both the 

] inch and the 1.5 inch diameter absorbers. When the rate of chemical 

reaction increases and becomes high, the resulting increase in the 

chemical absorption coefficient K implies that the resistance to mass 

transfer in the liquid phase has been reduced. Under such conditions 

it is not uncommon for the gas phase resistance to become important, 

in which case the absorption process would cease to be liquid phase 

controlling. Due consideration must therefore be given to ascertain 

whether the effect of the gas-phase resistance has to be accounted for 

when estimating the rate constant of the liquid phase reaction. 
Froessling's equation [20] for forced convection mass transfer 

in gases around spheres was used to estimate the mass transfer coefficient 

in the gas phase. The liquid and gas phase resistances were compared on 

the basis of the two-film theory. Froessling's equation may be repre- 


sented by the following expression. 





! 
kned 
Gaans a Wz omnes z 
5 el ee S. (2-65) 
Vv 
: ey, cms 
where ke = the gas phase mass transfer coefficient aa 
d. = diameter of the sphere absorber (cm) 
Do w= difTusiVvity Ole HD Inthe gas phase ome 
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>.) 
It 


Reynolds number based on sphere diameter 





n 
It 


‘ v 
% Schmidt number = D~ 

Vv 
The diffusivity (D.) of hydrogen-deuteride in hydrogen saturated with 
aminomethane vapour was estimated using the equation recommended by 


Hirschfelder, Bird and Spotz [21] which is given as 





=3 | l l 
D = 1.492 x 10 ec —— + —— (2-66) 
. Sed reel Mer ota 
ie leZ ane 
where D) = the diffusivity expressed in (ft7/hr) 
T, = absolute temperature in (°R) 
P. = absolute gas pressure in (atm) 
rio = mean collision radius of species | and 2 (A) 
sea iayae te A SEY 
(ro), = collision radius of species 1 (A) 
(ize = collision radius of species 2(A) 
| = collision integral for diffusion 
P at 
ly Ussexpressed aswdm(unctionvof ja—— 
2 
where 


E re E 
2 NE) on 
and k = Boltzmann's constant (erg/°K) 


€,,= energy of molecular interaction (ergs) 


ez 


Further details and tabulated values of Ly vse kT, for the use of 


ele 
equations (2-66) and (2-67) are given in Perry [22] where a modified form 


of equation (2-66) is reported. 
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Based on the two film theory, the overall liquid phase resistance 


may be expressed as 





l l l 
Va ertia & i (2-68) 
K. Kk, k aH 
where K = the overall liquid phase mass transfer coefficient <a 
rs P| wa cms 
Ke the liquid phase mass transfer coefficient 2a 
Ke = the gas-phase mass transfer coefficient expressed per unit 
of partial pressure driving force li) Moyes 
cm sec atm 
H' = The Henry's Law constant expressed as cose 
gm moles 


L} 
The mass transfer coefficients ke and ke may be related by the following 


expression 
t 
ee (2-69) 
G 
g 


Al x 


For a typical experiment the following values were obtained: 


eal diChigae 


oO 
It 


Vv 

So =20 5050 

R 9 felt) 0) 

e 

K = 0.245 cm/sec 

ke = 3.08 cm/sec 
Leese) 2cCe=e2olelonk 


Using the above values it may be seen that 
= 0.996 (2-70) 


This implied that the resistance to mass transfer was predominantly in 


the liquid phase under conditions where the highest reaction rates were 
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obtained and the assumption of a liquid phase controlled absorption was 
justified. 
6.2 Effect of a Change in y on k, 

The manner in which an error in the estimated equilibrium value 
would affect the slope of the plot given in Figure 9 and consequently 
would alter the value of the rate constant ky determined from that 
Slope needs to be investigated. It is important to realize that for 
a given value of temperature and catalyst concentration, the rate 
constant ky is fixed but an error in the estimate of ve could lead to 
yey, 


| VS metedealtal. 
Ye 


a misinterpretation of in( 
in 


Re-writing equation (2-57) as 


Ye) . 
f= ak ae (2-71) 
Y: ae, 


and differentiating w.r.t. Y, one obtains 
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substituting for from (2-71) 
Yes 
du v,wer' 3k” (2-71b) 
— come one = 2 = t Aa = 
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At. it 0 va ie and est + |] which leads to the conclusion that 
ou peel ;. - Poe ae “ . 
ave > 0. This implies that the ordinate u is invariant with ve at 

e 


small times. Since the ordinate u is invariant with Ds the plot of 


uvs. t is fixed and hence k* is invariant with patel — et UR MN S 
e 
may be shown rigorously by considering the second part of equation (2-71b) 


i.e. 
kk t 
as f oy 
era) e 
" (y = Vileeac 
bee escalate eat Rett (2-71c) 
eT a Je 
in ’e 


Taking the expression (2-7lc) to the limit t-+0; VieaY Ge aONe obtains 


(y. -y) , k*t (y. Ve )k* 
t-0 ed eer pam) Tees ee ee: = i" ="0 


2a 
OT nae) at va AMERY nay) OP YS 


In general, as long as (k*t) is a small number, the derivatives a. 
and — and also the resulting error due tons change in Y, are ops iin 
The ene of k* is usually smal] (the maximum value obtained was 

of the order of 3.3 x lo eee) and as a result low errors in NAS would 
not affect the results appreciably if the data obtained in the earlier 
part of the experiment are used. Generally, the error in estimating As, 
is not expected to be larger than +1 ppm and this was found to have a 
negligible effect on the final value of the rate constant Ky: An 


analysis was conducted wherein the equilibrium value was changed from 


22.0 to 25.0 ppm for a set of experimental data obtained at -11.7°C. 
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The corresponding values of the rate constant obtained are given in 


Table 3. 


With the above arguments in mind, if a plot of in ase t 


for an experiment indicated a deviation from linearity after a time 
't', the data prior to that time was used to estimate the rate constant k}- 


6.3 Modelling the Data of kK, vs. t and k, vs. ¢ 


Two models were tested for the data of k, vs. T and three models 
were proposed and tested for the kK, vs. c data. The best model for each 
set of data was selected on the basis of the Bayesian expected likelihood 
criteria which utilizes the Bayesian interpretation Olepnopsebila ty. s tnis 
approach to model discrimination has been used with considerable success 
by Singh [23]. The selected models were then combined to give a joint 
model incorporating both temperature and concentration and the best 
parameters for the model were determined by a non-linear estimation 
technique suggested by Marquardt [24]. 

6.3.1 The Bayesian Approach to Model Discrimination. Consider that 
a set of m models which contain the true model describing the observed 
data is selected after preliminary screening. Let the set of m models 


be represented as 


Y] a o,(, ay) + oy 
yy = ox, ap) + ey (2-72) 
i 
Vana (Saat em 
where 
x = the vector of independent variables 
a, = the vector of unknown parameters for the kth model 
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TABLE 3 


Effect of a Change in y 
ec —) 


on the Rate Constant Ky 
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yy the dependent variable 


ey the error term in the kth model. 


The truth (or conversely the uncertainty) of any model form k may be 
expressed quantitatively by a discrete probability PL: The set of 


P's then satisfy the following relationships 


m 
OR FP ca) Osan) ppasale0 (2-73) 


The Bayesian probabilities Pp. are indicative of the current state of 
information about a certain quantity and are subject to updating with 
the acquisition of additional information by the yse of the Baye's 

Theorem which states that given n sets of observations (yyr> Xi5 i=l, 


a] 
2---n) theposterior probability for the jth model P; is expressed as 


‘ +L (yet..) 
= aaa jgea lie em) (2-74) 
j Py LO 1e) 


P; is the probability prior to observing the data Y| and L(y /o;) is the 
expected likelihood of observing the data, given that the model form be 
is the correct one. 

In absence of any conclusive prior information it is generally 
assumed that all models are equally likely which implies that the prior 
probabilities are taken to be all equal. Similarly, the prior parameter 
values are assumed to be represented by an improper uniform distribution 
implying that they are distributed over the range -~ to +m. For the 
case where the parameter prior is represented by an improper uniform 


distribution the expected likelihood for the jth model may be represented 


by the following expression derived by Singh [23]. 
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2 =n/2 1 r*/2 SA 
L(y, |¢.) = (210~.) () exp (2-75) 
5 q 2 2 
20 « 
J 
where r* = the number of parameters in the jth model. 
5 = the minimum sum of errors squared. 
n 
a ) eb di ‘calcul d Ne 
i=] erved i calculated i 
minimum 
a. = the parameter vector which minimizes the sum of errors 
J Squared. 
fete reg? 
fo} ‘ = (Caer = the unbiased estimate of the error variance. 


The 95% confidence interval for the estimated parameters may be expressed 


as A 
= + - 
“ite as ie wen 
where a, ae the kth element of the least square estimated vector a; 
for the jth model. 
2 aver the kth diagonal element of the variance covariance matrix 
’ 


for the jth model. 
The variance co-variance matrix may be computed from the following relation- 
ships 


oe ES ea (2-77) 


nok (aia @ (2-78) 
and x; is the (nxr) design matrix given by 
x; Sets Claw e2sn ence Kem 2onecr) (2-79) 


The elements of the design matrix are simply the first derivatives of 
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the model with respect to the parameters evaluated at the least square 


estimate Spe They are represented as 


a (x., a;) 
7} ik = la BSUS ot eat ce ” (2-80) 
if oa. a.=a, 
jk Ss. zi 


6.3.2 The Models 


k, vs. Temperature. The models attempted for a fit to 
the data represented as In ky VS. i. were 
mei 
NE ie oe (2259) 
fa 
or In kK = In Ke BY ae (2-61) 


Equation (2-59) is the Arrhenius equation 


See / Rem 
2)iekie= KET Te g (2-60) 
; ES 
or In ky = In ko =a) af "RT 
g 
ky . ES 
Lo In ae = In KO RT (2-62) 
g 


which corresponds to the Transition State Theory with an exponent of | 
on T. The data was fit to equations (2-61) and (2-62). The values of 


i} 
Koy ake aha ES obtained for the two models are reported as equations 


op Xa: 
(2-63) and (2-64). The posterior probabilities of the two models, 

the 95% confidence limits for the parameters and the standard deviation 
of the fit are given in Table 4. The parameter estimates have also 


been included. 


The posterior probabilities of the models indicate that in the 
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TABLE 4 


Model Estimates for Rate 
Constant vs Temperature Data 
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range of the data the transition state theory equation is only marginally 
superior to the Arrhenius equation. This implies that the In T term in 
the second model contributes very little to the overall representation 
of the data. A variance equal to the variance of the model fit was 
used to estimate the posterior probabilities. 
ky vs. Concentration. The models selected to represent the data 
of ky vs. concentration were limited to those that satisfied some of 
the basic criteria observed by previous investigators and during the 
course of this investigation. These criteria may be stated as follows: 
a) The curve of k) vs. ¢ should pass through the point (0,0) since 
no reaction was observed in absence of the catalyst. 
b) At very low catalyst concentrations, the dependency of kK VS.¢ 
should become almost linear [1,2]. 
c) As the concentration of the catalyst is increased the curve 
of kK) vs. C should level off, eventually tending to become independent 
Oe concentration vlls2\- 
The data obtained during this investigation also suggests that all three 
of the above criteria hold. There is evidence however, (Fig. 11) that 


at low concentrations the plot of Kk) vs. c is not quite linear. Based 


on the above criteria the following models were proposed. 


(i) k, wh ce (2-81) 
-a Cc 

ii Ts Send = a (2-82) 
-a ek 

(iii) k, 5 a,3(1 -e a ) (2-83) 


It may be noted that model (i) satisfies only criteria (a) and (b) but 


not (c); it was however selected so that a polynomial fit of the data 
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over the concentration range of this investigation may be analyzed. 
Model (ii) satisfies all the three criteria while model (iii) satisfies 
criteria (a) and (c) but suggests that at very low 'c' the rate constant 


would be approximated by 


3¢ (2-83a) 


This was introduced to ascertain whether the linear relationship was 
valid for the data obtained during the course of this work. If the 
linearity did indeed hold the constant 233 was expected to be 1.0. 

The kK, vs. c data was all obtained in the temperature range of 
-11.0°C to -12.5°C. It was all reduced to -11.0°C using the transition 
State theory model. The resulting values of the least square parameter 
estimates, the 95% confidence intervals, the standard deviations of the 
fits and the posterfor probabilities of the models are given in Table 5. 

In order to compute the posterior probabilities of the models the 
variance of the data experimentally observed should be used. This is 
in general smaller than the variance predicted by the model fit because 
the latter incorporates uncertainties due to the model form. Based on 
all the data taken for ky VSeal and ky vs. C, a variance of 56.25 was 
used which corresponds to a 20 value of 15.0 or a +20 limit (95% 
confidence) of 30.0. 

lt is interesting to note that as far as models (i) and (ii) are 
concerned the closeness of the fit reflected by the standard deviations 
is approximately the same, with model (i) faring somewhat better. The 
main deficiency of model (i) apparently is reflected in its inability 
to predict a levelling off with increase of concentration. Model Uris 


however, is designed to level off but is apparently quite poor because 





al 
-besylsre od yam noisepit2eval alas tO epnet noize tsNSoNSS nid) sae a 
. ay 
egiieizee (111) lsbom slldw sivei iso aoxdd od fev eeptetgee (18) Tecen 
tratenog sie. fz '>' wo! yrov 16 96d atesppue dud (9) baa (e)>oeinesia9 
vd basamtnovqqe:edublube | 
Ff 6 Ke oe ke ia : 
(e€S-5) . ee 36 96 (8 = 7 
tow qirenoiseles rsenil sd iatderw niesseoes od bsoubotsn! dew e18T 
sit +1 .advow efdt Jo se 1uo9 sf7 paltub bentletdo si#sb sd 307 bi tay 
0.1 ad of be3oSqxe 2ow pps Jnetenos ory bio beebai bib ytinsentl 
iG Sons. SIwIsI4sqned ond of} benisico Lis 26w 636b°O' sev he sit 
noitizenss3 eft prieu 3°09. 1)- oF beaubey bhe eew of -.9°S.8)= eseoro. ice 
TsTsmsexeq Sisupe tesel oid to goufev anit{pest oAT —. lebom yaeerns7asaze 
ent to sAohinivas bsabnece oft ,2lavieteal esnebitaos $2 oft .estemités . 
2 aideT ni savip oss elsbom od? to zelzilidsedorg soinereed ead bas 2317 
onj se) sbom Sis te ihilidsdeo'g Yolsedeog 47s Stuqpon od ratho nl 
2i 2idy .been sd bluode bevisedo yi lesnemiveqxs steb afd te esone}iaV 
2u62ed 3)1 Tobom ej. yd bessibeag sansiisv eid neniy veiieme: te1sngp a) ies 
7 


no bsesd .mrot [sbom ant o3 sub estigletdiesnu esses 0gtooni 183381 snd 

26w 25,02 to evininaey 6.2 .e¥ -A bee T.ev )¢ 10% nedted steb orld ie 

$28) timit wStis-10.0:2f Yo eulev BS» oF abnegearnas daidwibeal 
0-08 Fo, (satebi tno ; 

e768 (17) bas (1) elobom es ae es tert mal os pene eid : 


peepee: 


ae 





a 


152. 


TABLE 5 


Model Estimates for Rate 
Constant vs. Catalyst Concentration Data 


Oo 
+20 aoe (Model 
20h P98 ie nestor 











MODEL | 

EQUATION | 56.269 | -1.274 0.0188 | 3.846 | 0.261 250411 
(2-81) 

MODEL 2 

EQUATION | 863.61 | 0.0689 0.0004 | 143.84 | 0.0189 15.550 
(2-82) 







MODEL 3 
EQUATION 
(2-83) 
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it does not level off as fast as is required by the data. The posterior 
probabilities indicate that model (iii) is 98.08% probable while models 
(i) and (ii) have a 1.88% and 0.04% chance respectively, of being the 
correct models to represent the data. Model (iii) is therefore clearly 
the most likely one from the set of models chosen to represent the data. 
Its standard deviation also is the least of three with a 95% confidence 
interval for the data of about £3%. It may also be noticed that the 
value of a3, = 1.357 indicates a slight departure of the model from 
linearity at low catalyst concentrations. 

An investigation of equations (2-60) and (2-83) suggests that if 
a joint model form was selected to penredert the data so that the effect 
of temperature and concentration could be predicted from one equation, 


it becomes apparent that the parameter a,., of equation (2-83) could be 


13 
the only temperature dependent term. With this in mind the following 


model form may be proposed. 


E 


a a 
- <= , 34 
R_T a54¢ 


= g = -84 
ky ah Te (l -e ) (2-84) 
The above equation incorporates four parameters that need to be estimated 
namely ary: Es ao, and 234° A non-linear least square estimation based 
on the method of Marquardt [24] led to the following values of the 


parameters and their 95% confidence limits. 


Shih 883.09 + 593.34 
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Standard deviation of the fit = 11.943. 
The above given standard deviation of the fit corresponds to a 95% 
confidence interval of about +47 

lt may be noticed that the Senne anh has a very high band 
width whereas its counterpart in the linearized equation (2-62) 
namely k. can be estimated with much better confidence. This may be 
attributed to the fact that the amount of data available for estimation 
of a four parameter non-linear model is limited. A larger number of 
points would improve the confidence level of the parameters. 

Equation (2-84) was used to calculate Ky vs. T data at a catalyst 
concentration of 15.80 gms K/kg amine and also ky vs. c data at a 
temperature of -11.0°C. The experimental data is compared with 
calculated values in Tables 6 and 7 respectively. Figure 12 is a 
series of plots of k) vs. ¢ as predicted by equation (2-84) over a 
concentration range of 0 to 30 gms K/kg amine and a temperature range 
Gums i0s GC to -30°G. 

6.4 Comparison with Existing Data 

The investigations of Rochard [2] and Bar-Eli and Klein [1] are the 
only ones to date dealing with hydrogen amine systems in the A her: 
of metal alkali amide catalysts. Bar-Eli and Klein [1] used lithium 
methyl amide to measure the rates of exchange of the deuterium-amino- 
methane reaction. Since the effect of the cation on the rate of exchange 
of hydrogen-amine systems is not known it is not possible to compare 
the results of Bar-Eli and Klein [1] with the data of this investigation. 
The above authors [1], however, quote an Arrhenius activation energy 


of 3.6 k cal/gm mole obtained for the D,-CH.NH LiNHCH, system. This 
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TABLE 6 


Comparison of Calculated vs. Observed Data 
of Rate Constant vs. Temperature using the Joint Model 


Catalyst Concentration c = 15.80 (gms K/kg Amine) 


-| 
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TAG Ese), 


Comparison of Calculated vs. Observed Data 
of Rate Constant vs. Catalyst Concentration using the Joint Model 


Temperature = -11.0°C 
% Error 
=| -] ena 
c kjp(sec ) cE eae Ex 100 
Date (gms K/kg Amine) Eqn. (2-84) (Experimental) 1E 
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value compares quite closely with the value of 3.5 k cal/gm mole 
obtained from this investigation but the values of the rate constant 
ky are substantially higher than the corresponding values obtained by 
Bar-Eli and Klein [1]. 

The data of Rochard [2] for the hydrogen deuteride-aminomethane 
exchange was obtained in a stirred reaction vessel at temperatures 
of -62.2°C, -77.2°C and -90.0°C and potassium methyl amide concentra- 
tions ranging from 0 to 70 mm/1 (0 to 3.5 gms K/kg amine). The low 
temperature and the low catalyst concentrations were apparently 
necessary because of difficulties encountered in eliminating liquid 
phase mass transfer resistance. This difficulty was particularly felt 
for the HD-CHNH, system wherein stirrer speeds had to be increased 
beyond 3000 rpm as a result of which temperature instabilities were 
also encountered. Rochard's data suggests that the rate constant Kk) 
varies linearly with catalyst concentration only up to about 20 mm/1 
at -90.0°C and about 10 mm/1 at -77.2°C. Subsequently, the curve of 
kK, vs. c levels off suggesting an invariance of ky with c for higher 
catalyst concentrations. This effect was attributed to the variation 
of the degree of dissociation of the PMA catalyst with total catalyst 
concentration. The values of the rate constant Ky reported by Rochard 
are up to 25 times lower than the corresponding values obtained by 
extrapolating the data of this investigation in temperature and concen- 
tration. In light of the fact that the present data indicates a 
linear dependence of k, with c up to about 5.0 gms K/kg amine the curves 
of ky vs. c obtained by Rochard [2] were linearly extrapolated from 0.45 
and 1.0 gms K/kg amine at -77.2°C and -90.0°C respectively to 5.0 


k 


gms K/kg amine and the resulting values of ky were plotted as In ~ 
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VS. 7 The slope of this plot suggests an activation energy of 6.0 
k cal/gm mole as compared to the value of 6.8 k cal/gm mole obtained 
by the above authors [2] through analysis of the data without extra- 
polation. Table 8 shows a comparison of the values of ky obtained 

by extending the above plot to a temperature range of -10°C to -30°C 
with the values calculated by the joint model represented as equation 
(2-84). It may be noted that the discrepancy in the two sets of data 
has been reduced to a factor ranging from 1.3 to 2.2. This could be 
attributed to the different values of the activation energy used and 
insufficient knowledge of the variation of the degree of dissociation 
of the amide catalyst with total catalyst concentration and with 
temperature. This analysis, however, suggests that in addition to 

the above mentioned uncertainties the liquid phase mass transfer 
resistance might not have been eliminated from the work of Rochard 
[2]. Additional work is required to relate the dissociation of the 
catalyst with temperature and total catalyst concentration in order to 
interpret the data on a more sound basis. 

The rates of exchange of the hydrogen-aminomethane system have 
been reported to be an order of magntiude higher than those for the 
analogous hydrogen-ammonia exchange [1,2]. Bar-Eli and Klein [1] 
used an indirect method to ke the rates of the two exchange 

] 


processes and indicate that CH.NH 


BZ 
[kj] 


NR 


are in the range of 10 to 35 at temperatures ranging from -48°C to -66°C 
respectively. The above authors however, could not adequately characterize 


the catalyst for the aminomethane exchange process but suggest that the 


Ral 






1 : 
0.3 to yo1sne nélisvidse me 2sésppue tolg 2ids to agota ent a, «eV 


bonisido sicm m@\iso + 6.6 Io sulev andy oF beteqma> 26 stom mpe\lso a 
“sixes 1uonyiw ei6b art Yo eleylens deuerdy [8] 210mtve svods ery yd 
baniaido ,4 to aeuléyv ory to noelieqmos 6 ewode 6 sider snolisiog 


J°OE- 0) 3°Ol- 30 spnay er uTe e0med & of sorq avods af? pnibnaixe ya 


ad 
nok EUS e6 baingestas 1 isbom Jniol of} yd beteiuslso a2ouley oa} noiw 
6isb lo 2242 awd oft ni yanseqs roe|b sdt tend badgon ed yam Zi , (88S) ; 
s¢ blooa eidT .&.8 of €& | mov? enloney tofse? se 63 beapbet nesd ean 
big bees voles asitevitos ofd jo eouvlev iné7e7718 edz 02 bezudi3376 
. fei so2e16 70 satosb sit To noistcivay Sf? Yo epbelword tnsioitivent 4 
njiw 0n& noite Nedno> Jeyiszoo ‘eter Ablw Jeviaties abime a3 to ) 
OJ noldibbe ni seco elaseppuve , isvewor ,eleyviens e2iql Ss iuJeieane? | 
Jaan zzum ee6rfe biupi!l srl ceitnisdveony banoijaem avods sz q 
biafiooh To Arow edd movt bateninits sesd aved Jon Jnpim sonsseiesy 
ait to mnolisisoseib ad% sisher o3 betiuper 2% drow fanoizibbA LIS) : 
93 yabyo ni noitsxtnssnos J2yisiso [8903 bee e1useseoms? Ai iw geviesas ; 
.vized brude ssom @.no sisb ani Jaiqieini - ( 
SV5n MaJzye snedsomonims-odporbyn aid to epnataxe To 25367 eit ; Z 
SfAy 1Ot seo neria isipi obuidnpem to sebI0 ne ed of betteger nsed | a 


Lil) wield bos ifs+v6ee . [8,1] ie | nonma-nagarbya zuopolsns 

SenEHOxS Ow) 93 Yo keIe7 ois a of ‘portzgm josatond gti . 
A) * a 

a — eure ot 





160. 


TABLE 8 


Comparison of the Rate Constant Kk predicted by the Joint 
Model with kK obtained by a Linear Extrapolation of Rochard's [2] Data 


Catalyst Concentration = 5.0 gms K/kg Amine 
k 


1 This Work k ky This work 


Eqn(2-84) (sec) Ky Rochard 





reported ratios are perhaps the lower limits of the expected values. 
At -59°C in the linear range of catalyst concentration Rochard [2] 

indicates that the ratio of the rate constants for the two exchange 
processes is about 25. The catalysts used were KNHCH, and KNH., for 


the aminomethane and ammonia reactions respectively. The data of 


Bourke and Lee [3] corrected for diffusivity and corresponding to 5.2 
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161. 
gms K/kg ammonia (i.e. 4.5 gms KNH,/1it.ammonia) and 20°C was extra- 
polated to -10°C using an activation energy of 5.8 kcal/gm mole. 

The rate constant thus obtained was compared to the value obtained 
from this investigation as calculated by equation (2-84) at similar 
conditions of temperature and catalyst concentration. A similar 
comparison was also made at a catalyst concentration of 20.9 gms K/kg 


ammonia (18.0 gms KNH,/lit. ammonia) and -10°C. The ratio Tkyloy NH 


3NH, 
Tk yyy 
3 


was found to vary from 56 at the lower to 37 at the higher concentrations. 
This decline was due to the fact that the linear increase of ky with c 
extends to higher values of c for ammonia than for aminomethane. It 
is estimated that these ratios will increase if similar comparisons 
were made at lower temperatures because of the difference between the 
activation energies of the two exchange processes. 
6.5 Accuracy of the Data 

An error analysis of the data based on the method of interpretation 
used for this study is given in Appendix C2. The above analysis 
indicated that the accuracy of the values of the rate constant is 
dependent largely on the accuracy with which the diffusivity of HD 
in aminomethane, the hydrogen solubility and the gas analysis are 
known. Aminomethane density, total gas volume and the time of sampling 
are also significant factors. Based on the accuracies of the above 
mentioned data it is estimated that the data obtained for the rate 


constant ky is reliable to within +102. 


, 
- ro} 
*619x6 u Of 61 
' 
hn 
é : % 
4 wae * 
‘ € 
r 
4 : 
. «wits 
: r{ j 
4 
. 
= 
NO17es Ss" ~ : 
rs ai : S 
iia 
‘H ) ¥ ie, 
"c ai ' 
: Sie 2A24'2viE6n c5p 
2 —— * : 
 gnifgqnmse to amiz aris 
i 


a 


S$vods sz 


mn efi yot benis 


7M 





we 


b is 


3 D te te 
‘OO 221967055 














\Sinomné.371\ HUA emp 2.8% 18.1) einomm 
‘ Ty 
a t ) VE ’ > q) 7 / 6 ie oni ew 3 
nD ve “'iS3dO evry Ineie 
Os usis> 26 nolsephiez 
6.7) hes te he | tev! 5'65 bs 33036571 =] t sibs Bt 
, ' F $ : eae be 
ent & JG SBM O2165 26W NOziyEqno> 
bi _ 
wisaeg Wb — + vo 2 7 
s G 'T\ MMM emp 0.51) sinomme 
a 
, 
c O77 YISV OF bivot eaw : - 
i 2 | 
H a de 4 } — wsh aif. 
Ge sub 26w snilfseb etAT 
é 'o eaulév isApid of ebnatxe 8 8 8] 
'‘«@ 
yun A sqr ~ ’ hater: 4 2 >| t 
4 : 25 67 OS 7EMmMises ei : 
wy ee kt es als) ce 
JIS 1TSQGMNS 2WOl J6 SDE SSW - 
7 
> > d a 4 * , ; ; 
ar © esioyens nolssviitos 
,e2 
) a3 to yoswosA 2.0. 
a ' _ 
é > bse : 2 ‘fsn6 017s AA af 
; j - 
. - . ~ " 
£ fi (buse 2idt r boeu 
. 10598 sf3 Jeni betseasibna 
i . iw yoslvsos edt no ylepie! tnsbneqsb 
19« 
+ j P| (arn wh art 
i \iiguioe neporbyd, addi start oma siate ae A! 


worst a 
ikea 






xa lov l6302 


é€&> 





soenrony ik 
* ' 


1 -tnB9 Vingizoe 


6 afd Ao & 








ido sieb 






£ 
jani3 basemites 21 4 


162. 


CHAPTER 7 


CONCLUSIONS 

The use of a single sphere absorber for measuring diffusivities and 
reaction rate data for the hydrogen-aminomethane and the hydrogen deu- 
teride-aminomethane-potassium methyl amide systems respectively, has 
been successfully demonstrated. 

Experience has been gained in the handling of the above systems 
and suitable operating techniques have been developed after overcoming 
difficult experimental problems. Information has been obtained on the 
acceptable materials of construction for use with the aminomethane- 
potassium methyl amide solution. 

The sphere absorber was operated in the fast reaction regime to 
obtain pseudo-first order rate constants for the hydrogen deuteride- 
aminomethane exchange between -10°C and -30°C, in the presence of pot- 
assium methyl amide as catalyst. The accuracy of the reported data is 
expected to be within + 10% and is dependent upon the accuracy of the 
diffusivity, the hydrogen solubility and the gas analysis of the data. 

The data obtained was found to be an order of magnitude higher than 
similar data obtained in a stirred type reactor [2] where mass transfer 
effects are expected to be limiting. In addition, the rate constants 
obtained were found to be at least 35 times the corresponding values 
reported for the hydrogen-ammonia exchange reaction. 

Additional work is required to relate the degree of dissociation 
of the PMA catalyst with temperature and total catalyst concentration 


before the data can be interpreted on a firm footing. 
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163. 
The experimental rate constants were correlated with temperature 


and catalyst concentration through a model which may be expressed as 


RT 7 1.3588 
k, = 883.09Te 9 (1-e eet’ (2-85) 


This model can be used to calculate the experimental rate data to 
within + 4% and is recommended for the estimation of rate constants 
where data is not yet available. 

The information obtained on the diffusivity and the rate constants 
can be employed to predict transfer coefficients and thus tray 
efficiencies that may be expected in low temperature contacting towers 


used in the heavy water industry. 
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NOMENCLATURE 


Constants in equation (2-83a) 

Constants in equation (2-81) 

Constants in equation (2-82) 

Constants in equation (2-83) 

Constants in equation (2-84) 

Vector of unknown parameters for models | to m. 

The parameter vector which minimizes the sum of errors 
square for the jth model 


Activity of the amide ion 


Activity of the potassium ion 

Interfacial area (ena 

Matrix defined by equation (2-78) 

Interfacial area of the sphere absorber corrected for 
the film thickness (cary 

Surface area of the dry sphere (ena) 

The catalyst concentration in PMA-aminomethane solution 
gms K ) 

kg Amine 


gm moles ) 
Concentration of HD in aminomethane liquid cc 


Concentration of HD in liquid aminomethane in equilibrium 


with CH.NHD ( senores ) 
5 cc 


Interfacial concentration of HD in liquid aminomethane 


gm mores) 
cc 


gm uel | 


Bulk concentration of HD in the liquid phase ( a 
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The concentration of HD in the liquid in physical 


equilibrium with the HD concentration in the bulk gas 


gm moles ) 
phase ( s.noles 


Absolute concentration of (D/H) in the gas sample (ppm) 
Absolute concentration of (D/H) in the principal 
standard (ppm) 

The variance-covariance matrix 


The concentration of catalyst in solution (eames) 


Diameter of the sphere absorber (CM) 

Molecular diffusion coefficient of HD in liquid amino- 
zi) 

sec 

Molecular diffusion coefficient of HD in a gas mixture 


2 
of hydrogen and aminomethane (=) 
Sec 


Activation energy of the reaction (2h) 
; gm mole 


gm voles 
sec’ 


The total gas absorption rate ( 


e 1 2 
The instantaneous gas absorption rate (2 rales 


: AL atm 
Henry's Law coefficient expressed as ace 


AGGracm ) 


; eae 
Henry's Law coefficient expressed as (= IGE 


Current signal corresponding to the Ho concentration in 
the gas sample (amperes) 
Current signal corresponding to the HD concentration in 


the gas sample (amperes) 


Collision integral for diffusion 
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Pseudo-First Order rate constant for the forward 
reaction eee 
Pseudo-First Order rate constant for the reverse 
reaction (secua) 

gs : -n l-m 
Specific reaction rate constant ([cat] (C.-C.) 
sect 
Pre-exponential factor for equation (2-59) (sec!) 


Specific rate constant corresponding to the undissociated 


[KNH.] specie gm moles rh boo 

Pre-exponential factor for equation (2-60) (scene) 
Specific rate constant based on the activity NHS (sec 
Specific rate constant based on the activity product 
NH; ayt (sec!) 


Boltzmann's Constant (erg/°K) 

Defined by equation (B2-5) 

The gas phase mass transfer coefficient (cm/sec) 

The gas phase coefficient expressed per unit of partial 


sek m moles 
pressure driving force eas |i slbed ln 


cm sec atm 
The overall chemical absorption coefficient (cm/sec) 
Chemical equilibrium constant defined by equation (2-16) 


cc 


Liquid flow rate over the sphere absorber ie 


The likelihood function expressed as equation (2-75) 
Number of exchangeable hydrogen atoms per molecule of 


amine 
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Number of exchangeable hydrogen atoms per molecule 
of hydrogen 


Number of moles of hydrogen in the gas phase 


gm moles 
; Gc 


Concentration of the amide ion in solution 
Discrete probabilities for models 1 to m 
Prior probability of the jth model 
Posterior probability of the jth model 
Partial pressure of hydrogen (atm) 


Partial pressure of hydrogen-deuteride (atm) 


Total absolute gas pressure (atm) 


gm moles 
cc sec 


The rate of chemical reaction 


The number of parameters in the jth model 
Collision radius of species 1 (A) 
Collision radius of species 2 (A) 

Mean collision radius of species 1 and 2 (A) = Ltr a * 
(r),/2 

Radius of the sphere absorber (cm) 

Ratio of currents (1,/1,) corresponding to HD and H, 
concentrations uncorrected for H3 contribution (i.e. 
measured at a value of I, other than zero) 

Ratio of (13/15) for the principal standard measured 


at the same value of |, as R- 


Universal gas constant cae or la 
ett ac 
Moe? tS gm mole °K gm mole °K 
Paved: 


Reynolds number based on the sphere diameter 7 
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Schmidt number > 
Vv 


The minimum sum of errors squared 

Time (sec) 

Maximum contact time (sec) 

Absolute temperature (°K) 

Absolute temperature (°R) 

Expression defined by equation (2-69) 
cm 


e f e e penned 
Velocity of the circulating gas ee 


The velocity vector representing velocities of fluid 
elements flowing in the three coordinate directions over 
cm 


the sphere absorber — 
sec 


Velocity of a fluid element along a stream line a 


Total volume of the system oecupied by gas (cc) 


Interfacial velocity of the fluid sl 
sec 


Interfacial velocity of the fluid at the sphere equator 


cm 
sec 


Coordinate representing distance into the liquid film 
flowing over the sphere absorber (cm) 

Vector of independent variables. 

Mole fraction of hydrogen deuteride in aminomethane 
Mole fraction of hydrogen deuteride in aminomethane at 
the gas-liquid interface 

Mole fraction of hydrogen deuteride in aminomethane in 
equilibrium with CH,NHD 
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The design matrix 
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Mole fraction of hydrogen deuteride in hydrogen 
Dependent variable for equation set (2-72) 

Mole fraction representing the equilibrium concentration 
of hydrogen deuteride in the hydrogen gas 

Mole fraction of hydrogen deuteride in hydrogen at time 


t=0 


The separation factor of D to H between aminomethane 
liquid and hydrogen gas at equilibrium 

Coordinate defined in Fig. 1 (radians) 

Coordinate defined in Fig. 1 (radians) 

Dimensionless expression defined by equation (2-43) 
Representation of model forms for models 1 to m 

Total thickness of the liquid film (cm) 

Liquid film thickness measured at the sphere equator (cm) 
Dimensionless ratio of gas phase concentrations of HD/H., 
defined by equation (2-48) 

Dimensionless ratio for gas standards defined by equation 
(2-49) 

Dimensionless variable defined by equation (2-32) 

Maximum value of w~ for the sphere 


Dimensionless expression defined by equation (2-41) 


gms moles 


Molar density of liquid aminomethane = 


Gas density ~~ 


Gas viscosity (poise) 
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Kinematic viscosity of the gas = 


Energy of molecular interaction (erg/°K) 
Error terms in models 1 to m 
The unbiased estimate of the error variance of the 


jth model 
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APPENDIX A2 
Detailed Operating Procedures 
A2.1 The Exchange Experiment 
The detailed operating procedure for conducting a typical hydrogen 
aminomethane exchange experiment using the single sphere absorber is 
outlined below: 

(1) Evacuate the entire gas-liquid exchange section and its connecting. 
lines to a vacuum of better than 5n. 

(2) Pressure the exchange system between the 'Whitey'' feed and liquid 
return valves with enriched hydrogen to a pressure of about 25 
cms Hg (gauge) 

(3) Close the lid of the refrigerated air bath and begin refrigerating 
the system. As the temperature of the system falls the pressure 
in the absorption chamber will also fall as noticed on the 
mercury manometer attached to the chamber. 

(4) Also start refrigerating the ethylene glycol-water baths around the 
receiving tanks, the intermediate tank and the catalyst wash 
return tank. Steps (3) and (4) should be done on the day prior 
to the experiment because it takes about 12 to 16 hours to attain 
steady temperatures. 

(5) About one hour before the experiment start pumping the methy] 
hydrate coolant through the shell side of the concentric pipe 
heat exchangers. This ensures that the lines which carry the 
liquid are cold when the experiment is started. 

(6) Check the pressure on the absorption chamber manometer, if the 
pressure reads less than 12 cms Hg, add some more enriched 


hydrogen; if in excess of 12 cms Hg, bleed off some gas so as 
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to set the pressure to the above mentioned value. 

(7) Start the hydrogen circulator, the gas rotameter float should indicate 
a flow rate. The hydrogen circulator is driven by an air motor. 
The gas circulation rate can be controlled by adjusting the 
speed of the air motor. This in turn may be accomplished by 
varying the air pressure to the motor. Adjust the hydrogen 
circulation rate till the rotameter reads between 15 to 20% 
of full scale which corresponds to an approximate flow rate of 
about 100 cc/sec at the conditions of the system. 

(8) Introduce some pure dry aminomethane into the saturator from the 
dry aminomethane tank. The liquid would come in under a vapour 
pressure differential when the saturator solenoid valve is 
opened. Hydrogen gas will be seen bubbling through the amino- 
methane column in the saturator. 

(9) Open the distributor solenoid valve. The pressure of the system 
will be seen rising due to an increase of aminomethane partial 
pressure in the circulating hydrogen and due to a slight temperature 
rise due to the heat of the distributor solenoid valve. 

(10) Open all the valves in the return line downstream of the ''Whitey'' 
liquid return valve upto and including the valve isolating the 
receiving tanks. This ensures that the return line downstream 
of the 'Whitey'' valve is at the pressure of the receiving tanks 
which should be the pressure corresponding to the vapour pressure 
of the aminomethane at the temperature of the receiving tanks. 

(11) Start the vacuum pump aaa the gas-sampling section and evacuate 


its lines including the traps and the sample manometer (3) (Fig. 6). 
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Surround the traps of the glass section (1) (Fig. 6) with liquid 
nitrogen and immerse the trap between the oil manometer and the 
gas absorption chamber manometer in a mixture of ice and water. 

Turn on the power to the electric timer and the potentiometer 
which is used to measure temperatures of the exchange section. 
The temperatures are measured by copper-constantan thermocouples 
located as shown in Figure 2. The cold junction of these thermo- 
couples is immersed in a mixture of ice and water. 

Attach a gas sampling bomb at the sampling location (16) (Fig. 6) 
and evacuate it both on the bomb and on the valve side. 

Approximately an hour after the hydrogen circulator is started, 
the manometer serving the gas absorption chamber should reach 
a steady pressure, indicating that the hydrogen is saturated 
with aminomethane. 

Take the initial gas sample in the following manner. All numbers 
in parenthesis refer to equipment locations on Fig. 6. 

(a) ltsolate the vacuum pump of the gas-sampling section by closing 
the vacuum valve (12). Also isolate the sample bomb by closing 
the glass valve (14) between the steel sample line and the glass 
expansion lines. 

(b) Open the sample valve (13) to introduce a gas sample into the 
steel line. Then close this valve. 

(c) Open the vacuum pump valve (12). This would evacuate the gas 
sample trapped as in (b). One gas sample is discarded to 
eliminate gas in any dead volumes and thus obtain a represent- 


ative sample. Then close vacuum valve (12) "after about’ 30) sec. 
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(d) Open the glass valve (14) and gradually crack open the sample 
valve (13) till the sample pressure manometer (3) records 
about 3 inches Hg pressure (absolute) in the sample bomb. 
Close sample valve (13). 

(e) Close the sample bomb glass valve (2) and then close valve (16). 
(f) Remove the sample bomb and open vacuum valve (12) to evacuate 
gas pressure in the glass line. The mercury limbs of the 
sample manometer (3) should come back to zero. The sampling 

section is now ready for the next sample. Attach a new 
sample bomb and open valves (2) and (16). 

Open the liquid feed valve and introduce PMA solution into the 
exchange system. Maintain the flow rate at a level higher than 
the discharge rate of the nozzle. This will cause the PMA 
solution to accumulate in the distributor. 

Step (17) will cause the system pressure to increase substantially 
since the entry of the liquid results in an effective compression 
of the vapour space. 

When the jet (nozzle) starts flowing and the sphere absorber is 
completely wetted by the PMA solution, start the timer. 

When the liquid level appears in the take-off tube, open the valves 
in the liquid return line and adjust the 'Whitey'! control valve 
‘$0 as to maintain the liquid level in the take-off tube. 

When the level of the PMA solution in the liquid distributor 
reaches the glass tube above the distributor, close the distri- 
butor solenoid valve and adjust the feed control valve to main- 
tain a flow rate of 7.0 to 7.5 on the liquid rotameter. This 


approximately corresponds to a liquid volumetric flow rate of 
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0.275 cc/sec. 

(22) Carry out the final adjustments with the liquid return control 
valve to set a take-off height of 1.25 cms. with the aid of a 
cathetometer. 

(23) Set the pressure in the reservoir tank of the gas sampling section 
equal to the system pressure. This may be done by first adjust- 
ing the pressure in the reservoir to within 1/2 cms Hg of the 
system pressure by visual inspection and then opening the valves 
connecting the oil manometer to the absorption chamber on one 
side and the reservoir on the other. The final adjustment may 
be done by balancing both the limbs of the oi] manometer to Zero. 

(24) Record the time at which all the adjustments for setting the 
conditions of the experiment have been completed. Since the 
system operation during this time did not conform to a fixed 
set of conditions, it may become necessary to shift the time 
axis by this amount during analysis of data. 

(25) Maintain the liquid flow rate and the take-off height at the 
desired values for the entire duration of the run. 

(26) After every 500 or 600 seconds take a gas sample in the following 
manner: Refer to Fig. 6. 

(a) A sample bomb attached to the gas-sampling system should 
already have been evacuated. 

(b) When the timer reads about 50 seconds less than the time at 
which a sample is desired, close valve (14) and the vacuum 
valve (12). 

(c) Open sample valve (13) and note the oil manometer indicate a 


depletion of pressure in the absorption chamber. Close sample 
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valve (13). 

(d) Discard the sample by evacuating through vacuum valve (12). 
Close this valve about five seconds before the desired 
sampling time. 

(e) Open the sample valve (13) at the desired sampling time and 
introduce sample in the steel line. Close the sample valve 
A 3))e 

(f) Open glass valve (14) to expand the sample into the sample 
bomb. The sample manometer will indicate the sample pressure, 

(g) Close sample bomb valve (2) and valve (16). 

(h) Remove sample bomb, open vacuum valve (12) and record the sample 
bomb number, the time at which sample was taken, the reading 
on the oil manometer and the temperatures of the gas and 
liquid leaving the absorption chamber. 

(i) Attach a new sample bomb and open valves (2) and (16) so that 
the bomb is adequately evacuated at the next sampling time. 

(27) Continue taking gas samples according to step (26) till a sufficient 
number of samples are available. About 10 to 15 samples were 
taken for each experiment conducted during this investigation. 

(28) If an experimental measurement of the gas phase equilibrium value 
is desired, the time interval between samples after the first 
10 samples should be increased to about 30 minutes and the 
experiment should be conducted for a total time of about 5 1/2 to 
6mhou rs: at) 10 2G. 

(29) Take the last gas sample in the same manner as the initial sample 
procedure outlined in step (16). 


(30) Close liquid feed valve and drain the PMA solution in the distributor 
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into the receiving tanks. 

When the distributor is drained, open the valves leading to the 
wash return tank and isolate the receiving tanks. 

Introduce aminomethane into the feed line and run the system 
in exactly the same fashion as for the PMA solution (steps 
(17), (19) and (20)). Continue feeding aminomethane til] 
distributor is full. 

Stop feeding aminomethane and drain the distributor into the wash 
return tank. 

Repeat steps (32) and (33) till the liquid in the distributor and 
the take-off tube is colour less. 

Stop the hydrogen circulator and record the level of the top of the 
saturator and the level of the aminomethane in the saturator using 
a cathetometer. These levels are required to estimate the vapour 
space in the saturator which must be included in the total gas 
volume of the system. 

When the colour of the liquid in the distributor has faded, isolate 
the wash return tank and the aminomethane feeding tanks. 

Introduce methyl hydrate into the liquid feed line and run the system 
with methyl hydrate in the usual fashion; the return alcohol should 
be sent to drain. 

When alcohol is introduced into the system the chamber pressure wil] 
fall rapidly due to the reaction of the remaining catalyst in 
solution and the alcohol. The liquid in the system will turn milky 
initially but will start clearing up as the product of reaction are 
taken out of the system. 

When the alcohol in the distributor and the take-off tube appears 


clear, the feed may be stopped and the system should be completely 
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drained till the absorption chamber pressure is about | or 2 cms 
Hg (gauge) . 

(40) Remove the liquid nitrogen baths from around the traps (1) (Fig. 6) 
of the gas-sampling section and let the frozen aminomethane be 
evacuated into the trap of the vacuum pump. It can be discarded 
from there. 

(41) Shut down the coolant flow to the concentric pipe heat exchangers. 

(42) Transfer the cooling load to the feeding tanks if they are to serve 
as the receiving tanks for the next experiment and allow the 
receiving tanks to come to room temperature. 

(43) The experimental run is comptete at this point and the gas-liquid 
exchange system may be cleaned and completely evacuated in 
preparation for the next experiment. 

A2.2 PMA Solution Make-Up 

The following sequential procedure may be used to make up a charge 
of PMA catalyst solution in ami nome thane. 

(1) Estimate the quantity of aminomethane and potassium metal required 
for a desired charge of PMA solution by using stoichiometric 
relationships. 

(2) Assuming a 0.8% moisture content by weight in the aminomethane deter- 
mine the amount of lithium metal required to dry the amine using 
equation (2-45). | 

(3) Lithium is stored under toluene and is available in lump or stick 
form from Fischer Scientific Co. 

(4) Cut the required amount of lithium metal using a sharp knife and 
shave off the top skin thus exposing the shining surface of the 


metal. This operation should be performed under toluene. Cut 
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the cleaned lithium into small pieces. 

(5) Weigh the required amount of cleaned lithium pieces and immerse 
them in n-pentane to remove the film of toluene. 

(6) Introduce the lithium in a dry, clean steel vessel against a 
stream of argon. This is fone to ensure minimal contact with 
air and moisture. 

(7) Seal the port through which the lithium was introduced and evacuate 
the vessel to pump out the remanents of argon, and n-pentane. 

The lines connecting the vessel to the aminomethane cylinder, 
should also be evacuated at this stage. 

(8) When a vacuum of 1 yw has been achieved, immerse the vessel ina 
bath of liquid nitrogen and open the feed valve of the amino- 
methane cylinder. The aminomethane will distill into the steel 
vessel. 

(9) The level of liquid nitrogen should be maintained close to the 
top of the vessel making sure that the fittings through which 
the aminomethane is entering do not dip in the liquid nitrogen. 
This is done to prevent freezing up a line with aminomethane thus 
blocking entry into the tank. 

(10) The aminomethane cylinder should be placed on a balance so that the 
weight of aminomethane distilling can be monitored. 

(11) When the desired weight of aminomethane has been distilled into the 
vessel, close the amine feed valve and lift the vessel out of the 
liquid nitrogen bath. | 

(12) In order to limit the pressure buildup in the vessel due to hydrogen 


evolution a relief valve set to about 80 psig should be provided. 
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Allow the vessel to come to room temperature; the aminomethane 
will Regia to thaw out and this can be seen on a pressure gauge 
which will start showing an increase of pressure (when cold the 
pressure gauge will indicate a vacuum in the vessel). 

When the pressure gauge indicates a positive pressure, connect 
the relief valve into the system. 

The aminomethane will take about 12 hours to thaw out. The 
reaction of lithium with moisture in the amine begins as soon 
as first contact with liquid aminomethane is made. Since the 
reaction results in heat evolution, the thaw out process is 
accelerated and pressure buildup will occur fairly rapidly. 

Bleed the hydrogen pressure into the vent system periodically 
to maintain a pressure of about 30 psig, till the pressure 
does not rise anymore after bleeding. 

The aminomethane is now dry and ready for use. The process may 
take upto 2 days to reach this stage. 

Cut the estimated amount of potassium metal under mineral oil in 
the same manner as indicated in step (4). 

Follow steps (5) to (7) using the 316 SS storage vessel for the 
PMA solution. 

When a vacuum of 1 uv has been reached immerse the vessel in a bath 
of liquid nitrogen and flash distill the aminomethane dried with 
lithium. 

Follow steps (9) to (16) keeping the solution in the vessel 
circulating by the use of a polypropylene diaphragm pump. 


When the colour of the solution changes from blue to straw yellow 










ae 
7 
a ~ a 
UWIGIOGNSS Moor + » oF Ise2ov edz wollA (€1) 


£) 
* 


>? nined if iw 


a « 
: ft. J 
a i s22510 
{ 
a : 
if Sr q 
: - avi 
ie] 
. : of : 
5 
- 
- “a *) 
be jO13565" 
: 
j is (> 
? ' , 5S) 2 ,¢c 
' 
* ’ 
: ic 
. '> 
+ 5 
74 rn 1) 
4 * * 
' sUU 
: > 
we 
, , oo fal 
3 V < {© wey a7 
. . 
nolzuloe AMS 
* ! 
i ane eat ne ~AoaA as 2 > P - a - i ij 
jJec & TT! 22 m3 2 ism: O8N3601 As sa a | 'O MUNOBY 6 nonwW (OS. 
= Tas 
Ss 






tj er’ Saal a 


, : aS 4 
: ia waa bo le ie oy or. wv 
ie e25V any nm nous ql 2 ar 


183. 


the reaction may be considered to have gone to completion. 
A sight glass should be provided with the vessel so that the 


colour may be conveniently seen. 
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APPENDIX B2 
Sample Calculation 

The following procedure was followed to estimate the value of the 
reaction rate constant k (sec |) from the analysis of the gas samples 
taken during an experiment. The experiment conducted on 11th October 
1972 is used to illustrate the calculation procedure. 

The raw data sheet as obtained from the experiment is given below. 
During the experiment columns 1-5 and column 7 were recorded. The 
analysis of the gas samples was conducted subsequently and the results 
were transferred to column 6 of the data sheet. The numbers in column 
6 represented as R* are the ratios of the currents 13/1, uncorrected for 
the HS contribution and obtained at a constant value of I, = 24.0 x 10n 
amperes. 

For the closed system under consideration, the rate of depletion 
of D from the gas phase is essentially independent of system pressure 
and as such the readings in column 4 namely the oi] manometer differential 
pressure, were not needed but were recorded in case some correction had 
to be applied. 

The gas analysis was done on the mass spectrometer and the results 
were interpreted to give the absolute values of D/H. This was done with 
the aid of four known standards which were analyzed before and after 
the analysis of the gas samples obtained during the experiment. The 


values obtained from the mass spectrometer for the standards along 


with their absolute values are given in Table B2-2. 
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TABLE B2<2 


Mass Spectrometer Analysis for the Gas Standards 


R% 10° D/H Absolute(ppm) 
-51.6 1.5 
Before 61.3 99 
346.3 365 
ig ich A 47] 
-54.] lee5 
58.3 99 
After 
349.0 365 
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The analysis of the standards were converted in the following 


manner: 
R* 
Sy plows (a = ) 1000 (2-48) 
st 
cs | 
ae “ge ) 1000 (2-49) 
st 


where the terms of equations (2-48) and (2-49) have been defined in 
chapter 4. For the above experiment 


es = 346.3, 349.0 Ce = 365 
From above equations one‘obtains 


TABLE B2-3 
Converted Form of Gas Standards Analysis 


6 6 


DI st 
-1149.0 =e es] 
- 823.0 -728.8 
Before 0.0 0.0 
366.5 290.4 
-1155.0 sli 5) 
- 833.0 -728 .8 
After 0.0 O20 
SA 290.4 


The above data was fit to a linear and a quadratic equation using 
the linear least square criteria. The equations are of the form; 
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The results indicate that the linear equation (1) is quite adequate 


for the analysis. The following parameters were obtained: 


For Model (1) 8B 17.136 B 


WW 
—— 
. 
oa 
—S 
——S 


h 


For Model (2) By Tone Boe Wels 44 = Ons 37 econ 


The magntiude of B, indicates that the quadratic contribution is 


3 


negligible; in addition, the variance-covariance matrix of the parameters 


indicated that the 95% confidence limit of B, was +0.7 x 107! implying 


3 
that it is not reliable. A linear model was thus considered adequate. 


The gas analysis data was converted to the corresponding Sy 


values with equation (2-48) using an average value of Rey given by 


_ 346.3 + 349.0 _ 
Solent) ? 2 = RY) 


The linear model was used to get the corresponding values of boy? and 
using equation (2-49) the absolute values of the gas samples in ppm 
were obtained. The final values obtained are given below: 

TABLE B2-4 


Data of Absolute D/H Concentrations vs. Time 


Time y = D/H 
(sec) (ppm) 
0 527 
500 479 
1000 435 
1500 394 
2000 353 
2500 | ac 
3000 295 
3500 263 
4000 238 


4500 216 
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TABLE B2-4 (Cont.) 189. 


Time(sec) Y=D/H (PPM) 

5000 200 
5500 181 
6000 164 
6500 152 End of Run; Created Spill 
8000 62 
9000 4] 

10,000 ? 

11,000 27 

12,000 28 


Since the earliest part of the run (the period between the initial sample 
and t = 500 sec) includes some stabilizing effects like setting the 
take-off level, adjusting the flow rates, setting the oil manometer 
etc., the first point was usually neglected. The time axis was thus 
moved to T' = t - 500 (sec). At t = 6500 sec the liquid was allowed 
to spill out of the take-off tube into the gas absorption chamber. This 
was done to increase the area of contact and thus reach equilibrium 
faster. The data points to be used for the estimation of the rate 
constant are 

TABLE B2-5 


Data of Absolute D/H Concentration vs. Time 
Used for Rate Constant Estimation 


li y 
(sec) (ppm) 

ao 479 
500 435 
1000 394 
1500 353 
2000 322 
2500 293 
3000 263 
3500 238 


4000 216 
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TABLE B2-5 (Cont.) 190. 


TaeGsec) y (PPM) 
4500 200 
5000 181 
5500 164 
6000 152 


Deliberate spilling of the PMA solution to obtain the equilibrium value 
was not done as a routine procedure for every experiment. This was 
done a few times during the entire course of this investigation so 

that the equilibrium value obtained in this manner could be compared 
with the value obtained by a mass balance of the deuterated specie in 
the liquid phase. The separation factor for hydrogen-aminomethane 

was assumed to be the same as that for the hydrogen-ammonia system. 


The initial value of (D/H) was determined from an experiment 


aminomethane 


where the gas phase was taken to equilibrium. The concentration of 


was then corrected after every experiment by esti- 


SOW inonethane 
mating the liquid phase enrichment as calculated by the following 


material balance equation 


ae er OR oe (B2-3) 
l 
where Xe = D/H ratio. in liquid after the experiment (ppm) 
ae D/H ratio in liquid before the experiment (ppm) 
Y; = D/H ratio in the gas phase at the beginning of the experi- 
ment (ppm) 


Mi te D/H ratio in the gas phase at the end of the experiment (ppm) 
n= Total number of moles of gas in the system 
n. = Total number of moles of liquid that passed over the sphere 


absorber. 
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The value x for the next experiment was estimated by another material 


1 
balance which involves Xoo Nis X; 


i (previous) and the liquid remaining 


back in the feeding tanks (not used for the present run). 

Equation (B2-4) was used to calculate Y, as 
vic =p ee (B2-4) 
where See is the separation factor. 

The value of Ye obtained from this particular experiment by direct 
gas analysis was equal to 25 ppm. The PMA solution which was allowed 
to spill into the gas chamber was a stagnant pool of liquid which 
aided in depleting the gas phase of its deuterium content. The 
stagnant liquid was enriched in the process and the transfer of 
deuterium was eventually expected to take place from the liquid to the 
gas phase. A plot of y vs T would thus go through a minimum. The 
value of y corresponding to the minimum was taken as the equilibrium 
value. Fig. B2-1 shows a typical plot of this nature. The value of 
hs obtained with equation (B2-4) was 24.9 ppm. It may be mentioned 
here that an average value of x; and Xe or, for runs where the initial 
gas concentration Yi, was low, the use of xX, may be better justified 
to calculate Wess However, the very small differences in x and x 


(~1 ppm) make it immaterial whether equation (B2-4) is used with xs 
1 


Xe 


or an average of the two. 
US 2 Ae 
Yin Ye 


be estimated. These values are given below. Yin is the D/H concentra- 


Using a value of ee 24.9 ppm, the ordinate In ( 


1 
tion in the gas at T = Q. 
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Plot of y vs. t to Determine the Gas Phase Equilibrium Value 
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TABLE B2-6 


Data of Cumulative Depletion vs. Time 





In ( Ye 7 
Yin Ye (sec) 
0.0 0 
-0.1024 500 
-0.2080 1000 
-0.3256 1500 
-0.4250 2000 
0 CYA] 2500 
-0.6440 3000 
-0.7540 3500 
-0.8668 4000 
-0.9500 4500 
-1.0669 5000 
-1.1856 5500 


193. 


£2 




















3-£8 JISAT 


- 
omit. Vv Rol jai ged avite! 3 to al f @ = iat 
1 . av¥ ; = - 
(ose) (one ym ul 
a* nit. 3 


y 0.0 
= | Asoflo- 
0901 oROe . O- 
aes 28.0 
6005 025. 0- 
Obes rt82 0» 
— od. 08 
OGct ovat o- 
— 8388. O« 3 
| rae o02e. o- 
— eaao. f- | 
we 8281. i- 


2ets. 





194. 


The slope (-k ) of the least square line through the above data of 





Me \ 
In (———-) vs. T was found to be 
Yin Ye 
Saree OP hax. 105° 
eyes EAE Ripe alone rss) 


Referring to equations (2-57) and (2-58) 


& Ken On hee 
k = os (B2-5) 
Z 
where A, = Ay + Ay (B2-6) 
An = Area of the sphere absorber corrected for liquid film thick- 
ness gn) 
A, = Area of the liquid jet (cm*) 


Interfacial Area of Contact 
Sphere diameter = 3.81 cms 
Jet diameter = 0.0460 cms 
Length of Jet = 0.2650 cms 
Liquid film thickness 


at the sphere equator = 0.0066 cms 


“, Area of the dry sphere = nde Sl Gp 3.81° = 40.6 war 


Correcting for liquid film thickness by using equation (2-58) 


> 
Il 


2.58 x 0.0066 
Ne (i+ aeene0S 


45.6 (1.00894) = 46.0 cm” 


Area of Jet = md iL where di; = diameter of the jet (cm) 


Ay = 1 x 0.0460 x 0.265 = 0.0383 cm” 
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A, = 46.0 + 0.0383 = 46.04 om 
The same calculation on the computer resulted in A, = 46.05 cm? 


System Temperature 
The average temperature of the liquid was taken as being represen- 
tative. From the data sheet, average thermocouple output from t = 500 sec 
to t = 6500 sec = -0.4509 mv. 
Use linear interpolation of the standard E.M.F. data of copper-constantan 
thermocouples (Ref. 0°C) 
eee 1.0 mv 


l 


T, 


-0.417 


-12.0°C mv 


-0.455 


The equation of the mean curve representing the calibration of 


the thermocouples with respect to temperature is given by 


e, + 0.0064 
= B a 
Ss 0.9962 | (B2-7) 
where 

ss the standard E.M.F. generated by the thermocouple as listed 

in the ''Leeds and Northrup'' thermocouple tables (copper- 

constantan) (mv) 
en © MMIC actual E.M.F. generated by the thermocouple (mv) 


The actual calibration data for the above thermocouples is given in 


Moore's thesis [9]. 
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-0.4509 


_ -0.4509 + 0.0064 
5 0.9962 


-0.4463 mv 
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-0.4170 - (-0.4463) (_ 
TOS eo eEoy NO) 


Factor 


ee 0295 ee ° 
ime 


femppot Systemes: rl) .OSChep O0cJ/e Cems hdc 


Density of Aminomethane 


The following equation was used [9]. 


peeia tsa, + ae (B2-8) 
where a, = 0.905798 

a, = ~4.009472 x 10" 

a, = -1.457690 x 10°? 

T = Temperature °K 

0 = aminomethane density a 

oe 
T =273010 = 11.77 = 261.41 -K 


Using equation (B2-8) the density p may be found to be 


P gms 
6 =—0)./014 a 


Henry's Law Coefficient 


Using equations (2-13) and (2-14) one obtains 


In H, = 4.1403 + 1.6612 ser ir \. 0.1160 Gas a 


atm 
mole fraction 


or H, = 6625.7 


System Gas Volume 
The total volume occupied by the gas in the gas-liquid exchange 


system was estimated by expanding a known test volume of a gas at a 
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known pressure and temperature into the system. The initial pressure 


of the system was also known. The final pressure of the test volume 


and the system was recorded after the expansion. The following 


equation 


Since no 


equation 


Equation 


may then be used to estimate the unknown volume. 





a) Vy (B2-9) 


temperature change was observed during the test the above 
assumes a constant temperature. 


The unknown system volume (cc) 


Volume of the test tank and its connecting lines (cc) 


= Initial pressure of the test tank (cm, Hg) (absolute) 


Initial pressure of the system (cms Hg) (absolute) 


Final pressure of test tank and system after expansion (cm Hg) 
(absolute) 


(B2-9) was derived by a simple molar balance of the gas before 


and after expansion. The procedure involves the following: 


No. 


No. 


ee) 
of moles of gas in test tank + lines = at aN 
l 
ills 
of moles of gas in system before expansion = aT ana 
l 





Total number of moles of gas = (n, + ns) 


Total volume after expansion = Vy + Ue 
P_(V,+V_) 
Tain lie cS i 
Sap + re (Assumes no temperature change on expansion) 
l 
or PeVi tee we = Pi, ts Ares 


(pH mz) noiansqxe i93%6 metzeye bns 4eed fe59 To siueeetg feali 
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(P,-P,) 
or Mig iGxap Vy (B2-9) 


The system volume was estimated two or three times during the course of 
the investigation. Whenever a change in the apparatus constituting 

the gas-liquid exchange system was incorporated, a new test for gas 
volume was conducted. The last volume test was conducted on 18th May 


1972. The following results were obtained. 


Vy + connecting lines = 1796.3 cc 

Barometric Pressure = 703.8 mm Hg 
P, = 703.8 mm Hg (Barometric Pressure) 
1 = 171.37 mm Hg 
Pe = 388.53 mm Hg 


So loca, c 
Me = Qe: (1796.3) = 2607.8 cc 


A cathetometer was used to take the pressure readings from a 
mercury manometer. The above volume includes the volume of the liquid 
lines in the system that were not isolated and the volume of the liquid 
aminomethane that would be introduced into the saturator before each 


experiment. The volume of the liquid lines was esitmated as TO ATS 


Volume of gas in the system + volume 


of liquid aminomethane in saturator = 2607.8 - 43.0 = 2564.8 cc. 


The volume of liquid aminomethane in the saturator was estimated at the 
end of each experiment by measuring the height of the saturator not 
occupied by the liquid. A cathetometer was used for the purpose. 

For the experiment of the 11th October 1972 the aminomethane volume in 


the saturator was 49.2 cc. 
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Volume of gas in the system = 2564.8 - 49.2 = 2515.6 cc. 


Diffusivity of HD in Aminomethane 

The following method was used to obtain the value of the diffusivity 
of hydrogen-deuteride in aminomethane. 

Over a temperature range of -5°C to -35°C the diffusivity of 
HD-CHNH, and Ho7CHANH, was estimated at 5°C Ds Ieye using equation 
(1-49) given in the first part of this thesis. The mean ratio of 
Dip Oy over the entire temperature range was then estimated as being 
02 9:102: 

The experimental data of the diffusivity of Ho-CH NH, was taken 


from the mean line of Fig. 6 in part | of this thesis and a least 


square fit to the data of log, 4D ViS.6 ae was generated. The 
equation is given as 
log, ,D = -2.25564 - 0.432524 8 (B2-10) 
10 H, 
where 8 = a (B2-11) 


The 95% confidence interval of the parameters is better than 0.5%. 


The value of Dip was then taken to be given by 


= B = 
Dip 0.9102 af (B2-12) 


For the specific experiment under investigation 


1000 
109104, = -2.25564 - 0.432524 x se, 
= =3..9102 
(-3.9102) Cuchi a 
D =" 10.0 . 2123.0 °x% 10 —— 


Res 
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r -6 cm 
and Dua = 111.9 x 10 Cec 


The overall liquid phase mass transfer coefficient K 


estimated from equation (B2-5) 





é ‘ aT ihg! 
V H, 
_ 0.7014 gm moles 
aan 31.07 cc 
a cc atm 
Me =9O2 506 Gikmole °K 
k= 20 bleh) 
V0 F=225 yoeor(ce) 
= atm atm 
Ho = 6625.7 mole fraction ©" moles 
moles 
Z 
Ar = 46.05 (cm) 
= 20.2 bax 107? (eecm) 
a su Ho oes ue atm moles | 
L ied moles cm< 


_ 0.214 x 102 x 2515.6 x 6625.7 x 31.07 
46.05 x 0.7014 x 82.06 x 261.41 


2.14 x 2.5156 x 6.6257 x 3.107 
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may now be 


Cc 
moles 





cms 


sec 


(B2-5) 


moles °K 
cc atm 


Now for the fast reaction regime neglecting gas phase resistance 
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ie BL SERIO 
py APE: 


= 229.0 (ence. 
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Details of the calculation of bs by estimation of Liquid Enrichment 
Vapour Pressure of Aminomethane 

The following equation was used to express the aminomethane vapour 
pressure as a function of system temperature. 


C 
BS as e 
Mi SEN, es Cat + C, InT (B2-14) 


where C, = 100.90554 
C, = -5716.32506 
Gee 801022400 
C,, = -15.3000 
P = vapour pressure of aminomethane (atm.) 
T = temperature (°K) 
Using T = 261.41 °K one obtains 
P = 0.7806 atm. 
The partial pressure Pa of hydrogen was assumed to be given by 


P,=P_-P (B2-15) 


where PL = total system pressure at the beginning of the experiment. 
One could use an average total pressure of the system given by the 
mean of the initial. and final pressures but since this is used only 
to estimate the total number of moles of hydrogen in the system (of 
the order of 0.15) use of the initial pressure only is sufficient for 


a reasonably close estimate of the liquid enrichment 


ae. (265s Ace = 143.69 _ 
Pe eerste ne Bry |S ae t 
Py = 1.11 atm 


Number of moles of hydrogen in the system are 
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He eye ORT. x. 255156 : } 
i = RT nines sc DARI wee 0.1302 moles 


Consider the time of the run before the spill i.e. at 6500 sec. 
Average Liquid flow rate = 0.292 cc/sec 

The details of converting the liquid rotameter readings to the corres- 
ponding aminomethane flow rates are given in Moore's thesis [9]. 


Number of Moles of liquid passed over ball absorber is given 


_ 0.292 x 6500 x 0.7014 


by ny 31.07 


42.9 moles 


The final D/H concentration of the gas at time t = 6500 sec is 151.6 ppm. 
Y initial at time t =o is = 526.9 ppm 


(Y.-Y¢) Ms = (526.9 - 151.6) 0.1302 = 48.9 


Initial liquid concentration of D/H from analysis of previous 
experiment is = 114.3 ppm 


Using equation (B2-3) one obtains 


PARA SD SRN e tse Utshe 


fh oe 2.9 
= 114.3 + 2 
= 115.4 ppm. 


The equation of the separation factor is represented as 


Fae 237.59. | “ 
1091 9% aq = -0.2422 + * (2-21) 


i. 237.59 
Bede = 0.2422 oF OGLE 


-0.2422 + 0.9089 = 0.6667 
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= 199-6667 _ 4.6419 
eq 
x 
~The Pa ea 
te = Bog i T6419 = 24.9 ppm. 


The overall liquid concentration % used for the next experiment was 


estimated by using a linear mixing rule essentially of the same form 


as equation (B2-3). 


-£os 
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APPENDIX C2 


Error Analysis 


The standard error in a measured quantity may be defined by 


o} 
=o (C2-1) 
oe 


where o is the standard deviation of the measured quantity and n is 


the total number of observations. 
For any function 


Gr= sf (ASB tC ,D) (C2-2) 


The standard error in G is given by 
Z Ne 2 2 
Ziel of. 2, of 2 Of Zz of Z. : 
“aie (35) > a (22) Coes (26) Ce (+) ay (2-3) 


where G is the mean of G and Ayr Ups A and dy are the standard errors 


of the means of A,B,C,D. 


The equation used to express the depeltion of the hydrogen-deuter ide 


concentration with time is given as 


y-y KoAT Oak. 
In = -- +t ¢ (2-57) 
Yin Ye 2 


where AL includes the area of the liquid jet. 


For a fast first order reaction, the relationship of K with the 


rate constant kK, is given as 


(2-6) 


Expressing KP as a function of the other quantities 
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V H aye | 4) 
E re Mbal € Yin Ye 
OF A i F(V,H,AL,0) oT, tsYs¥Qo¥i 2) (C2-5) 


The universal gas constant R, is not included in (C2-5) because it 


aK, 
OR 
g 


is assumed invariant and as such = 0 (C2-6) 


Expressing equation (C2-5) in terms of standard errors one obtains 


2 2 2 2 
Demme te \ ated (cok \uo 282 He 2 i) 2 
ae - (25) 02 + (2) + (3 a + (a a 
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2 . 
+ (se oe (C2-7) 
oy. : 
in in 
The derivatives may be expressed as 
wt tometer in( pas) - - ae 
TFL g Yin “e 
oK K 
Lea: u oy 
a = (C2-9) 
pn UGA Hamer 78/3 *\ a (c2-10) 
Ay a ee "1 Y.-y 4 A, 
T AL Poe it inpee 
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iL Ab 3 
dK K 
Bee ol i 
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K 
an VH 1 : 
LE 5 y-y 
¢) AxtouReTt ( - IF a e (C2-14) 
Y Tees ee, (y ¥_) 0(- = 
in e 
SSE K (y-y;,,) 
ae y-Y ale 
In y - I ey 
in Ye 
aK K 
L L 
a z (C2-16) 
aN Tae 





eae. tec 


Substituting equations (C2-8) to (C2-16) in (C2-7) 
a2 
Oy ‘ 2 a 2 On 2 o 2 7 2 ae Z 
=~ - (*} + (-4) aoe ees e 
L T PL : : 
2 
SUN RAR Wo esc 
‘AMGEN (2 ACEI hone) LAG REVS! C2-17 
Y Ye R aN) Ye Yin Ye Yin Ye B 
Y-Ye 
where 8 = Jh [———— (C2-18) 
Yin Ye 


It may be noted that the terms involving 


ae Y,) 


are the ones 


which could contribute significantly to the total error if the data 


att > *@ (i.e. y > Y,) are used. 


It is advisable therefore, to shift 


the time scale so as to include data upto a time 't' where the 'y' 


value is still significantly different from Yer 


In general 
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most of the data taken during this investigation was obtained from 
experiments which were conducted for about two hours. 


Re-arranging equation (2-6) one obtains 


Me 
k, er (C2-19) 
2 ale ‘ 2 a8 ‘ 2 
sic “ky c= 3K as + 3D On (C2-20) 


The derivatives are given as 


ok, 2K) kK | 
——— = eetesinasl — iy —_ C2-21 
3K, D KP 

ak, ee k, 

ee ary (C2-22) 
oD D 

D 
a 72 2 


,{ aL (C 2-23) 
oe =) 
I K D 


Error in Gas Volume Estimation 


Re-writing equation (B2-9) 


P.-P 
ee ee v, (B2-9) 
Rey 








where 9V__ " = 5 Ss (C2-25) 
Lae eal jae 
WACPSEP ae 
SE ee teal (C2-26) 
oP, P| =P) (P_-P 
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av V : 

p> = OPP (C2-27) 
1 Poe 

aV Vv 

= (C2-28) 


Substituting the derivatives in (C2-24) 


2 2 
ty We " Preudllaals "Py a 
— =(5—-] + - +{=——] + l (C2-29) 
: meee PynP ed MPEP Cea e 
Representative estimates of the standard and fractional errors of the 


pertinent variables are: 


Op = + 0.2 mm 
] 
a =t O.l—mm = a 
Be A 
yy = am WO) Reker 
] 
a =a =a = + 1 ppm 
Mf Ye an 


A Fortin's barometer was used to measure P, and a cathetometer to measure 
l 


the heights of the two limbs of a mercury manometer to obtain Pe and AL 


The fractional errors are of the following order of magnitude. 
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Using the general equations derived above along with typical values 
Oley; Yin and Yor the fractional error in the overall mass transfer 
coefficient may be represented by 

ett iecs o( 10>") 


AE 


(C2-31) 


Thissredds to an estimated tractional error tn the rate constant of 


“Te 


eal 6x10— 


Ky 


‘ (C2-32) 


lt may be safely stated therefore that the rate constant values obtained 
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APPENDIX D2 
Detailed Experimental Data 
The detailed experimental data taken for the depletion of deuterium 
in the gas phase with time is given in Table D2-1. The gas phase 
compositions 'y' shown in the table represent the absolute ppm 
concentrations of D/H obtained after conversion from the raw data 
in the manner shown in Appendix B2. Table D2-1 may be used to obtain 


ky vs. T and kK) vs C data reported in Chapter 5 of this thesis. 
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TABLE 02-1 Zi 


Detailed Experimental Data 


Run No./Date 1/( 23/3/72) 2/(5 /4/72) 3/ (17/4/72) _ W(5/7/72) 
0 0 a bj a b 
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650 600 499 
967 1200 470 206 
Time vs. Gas 1240 1800 438 
1500 2400 410 
Phase Deuter- 2057 3000 379 
2550 4200 334 
ium Concentra- 3150 292 
3650 256 g 
tion Data 4250 214 Mean of ky values 
ape er obtained from both sets 
6850 was taken as representative. 
7850 
8850 
9850 


Total Interfacial 
















Area 20.67 (cm) 20.65 46.03 
Catalyst Concen- 
tration 15.80 (gmsK/kg . Amine) 15.80 15.80 


Henry's Law Coeff!8204.7(atm/mole frac.) 
Gas Phase Equili- 

brium Value 21.1 (PPM) 
System Tempera- 


6643.2 
24.8 


7006.0 
24.4 


















ture -30.11 (°C) -12.00 -16.66 
Total Initial 

Press. 1.476 (atm.Absolute) 1.866 1.68) 
Gas Volume in 

System 3108 cc 3091 2538 
HD-Amine Diffus- ef 5) * ie 

ivity 83.9x10 ~(cm*/sec) 111.5x10 104 .0x10 












Liquid Flow Rate |0.400 (cc/sec) 0.409 0.261 


Ami nome thane 








Density 0.7222 (gms/cc) 0.7017 0.7071 
SS Sl 
Run No./Date 5/(7/7/72) 6/ (11/7/72) 7/ (26/7/72) 8/(31/7/72) 


T (sec) | Y = D/H(PPM) T (sec) Y = D/H(PPM) 


Time vs. Gas 
Phase Deuter- 
jum Concentration 


Data 


Total Interfacial 


Area 46.03 (cm) 46.06 46 .06 46 .06 
Catalyst Concen- 

: tration 15.80 (gmsK/kg .Amine) - 15.80 15.80 15.80 
Henry's Law Coeff }6917.5(atm/mole frac.) 7644.9 7581.9 7571.5 
Gas Phase Equili- 

brium Value 24.5 (PPM) 22.9 23.1 : 23.1 
System Tempera- 

ture -15.55(°C) -24.15 -23.45 -23 .33 
Total Initial 

Press. 1.507 (atm Abs.) 1.435 1.481 1.514 
Gas Volume in 

System 2540 (cc) 2526 2523 2519 
HD-Amine Diffus- ees a a3 #3 

ivity 105.8x10 ~(cm°/sec) 92.6x10 93 .6x10 93.8x10 
Liquid Flow Rate |0.263 (cc/sec) 0.280 0.281 0.282 


Ami nome thane 
Density 0.7058 (gms/cc) 0.7156 0.7148 0.7147 
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Run No./Date 








Time vs. Gas 
Phase Deuter~ 
lum Concentra- 


tion Data 


Total Interfacial 
Area 

Catalyst Concen- 
tration 


Gas Phase Equili- 
brium Value 

System Tempera- 
ture 

Total Initial 
Press. 

Gas Volume in 
System 

HD-Amine Diffus- 
ivity 

Liquid Flow Rate 

Ami nomethane 
Density 


Run No./Date 


Time vs. Gas 
Phase Deuter- 
ium Concentra- 


tion Data 


Total Interfacial 
Area 

Catalyst Concen- 
tration 


Gas Phase Equil- 
ibrium Value 

System Tempera- 
ture 

Total Initial 
Press. 

Gas Volume in 
Sys tem 

HO-Amine Dif fus- 
ivity 

Liquid Flow Rate 

Aminomethane 
Density 


9/ (2/8/72) 10/ (28/8/72) ~-11/(04/9/72) | 12/ (22/9/72) 
T (sec) Y = D/H(PPH) 


46.04 (cm?) 


15.80 (gmsK/kg .Amine) 
Henry's Law Coeff}6704.6( atm/mole frac.) 


25.3 (PPM) 
-12.81 (°C) 
1.815 (atm.Abs.) 
2526 (cc) 


110.2107 (cm2/sec) 
0.277 (cc/sec) 


0.7026 (gms/cc) 


13/ (30/9/72) 


46.05 (em?) 


10.42 (gmsK/kg.Amine) 
Henry's Law Coeff}6597.4(atm/mole frac.) 


25.0 (PPM) 
-11.39 (°C) 
1.909 (atm.Abs.) 
2534 (cc) 


112.5x107°(cm2/see) 
0.298 (cc/sec) 


0.7010 (gms/cc) 


TABLE 02-1 Cont. 


46.04 


15.80 
6677 .6 


25.6 
-12.46 
1.442 
2531 


110.8x10- 
0.291 


46.05 


4.95 
6575.4 


25.0 
-11.09 

1.85) 

2525 


1130x1076 


0.296 
0.7006 


46.05 


15.80 
6618.8 


24.8 
-11.68 
1.890 
2521 


112.0x107° 
0.292 


T (sec) | Y = D/H(PPM) 


46.05 











4.95 
6625.7 


24.9 
“VW.77 
1.891 
2516 


111.9x1076 
0.292 


0.7014 


eM IRS 


46.05 


10.42 
6604.9 


24.9 
“11.49 
1.900 
2524 


112.3x107° 


0.293 
0.7011 
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